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ABSTRACT 


This  technical  report  covers  Task  II  of  a  two-task  investi¬ 
gation  into  the  application  and  implementation  of  advanced  discrete 
guidance  and  control  theory  and  analysis  to  a  United  States  Amor 
missile  system  of  the  future.  Task  I  was  completed  during  the 
period  7  August  through  30  September  1979  and  is  described  in 
Applied  Control  Systems  Dynamics  Company' Report  ACSD-78-2.  dated 
September  1978.  £ 

During  Task  H  a  stability  analysis,  summary,  and  implementa¬ 
tion  associated  with  the  discrete  control  for  a  missile  system  was 
performed  via  the  parameter  space  technique.  Emphasis  was  placed 
upon  the  development  of  analytical  tools  to  perform  the  design  and 
analysis  of  the  system.  The  three  design  tools  developed  under 
this  contract  are  described  r 


SAM*  ,  an  Alternative  to  Sampled -Data  Signal  Flow  Graphs 


Determination  of  Digital  Control  System  Response  by 
Cross-Multiplication ;  (f 


<%  Sampled -Data  Analysis  in  Parameter  Space. 

The  three  techniques  have  been  cast  in  a  form  amenable  to  imple¬ 
mentation  within  digital  computers  or  desktop  calculators,  f - 

Also,  during  Task  II,  a  technical  seminar  was  conducted  for 
selected  members  of  the  Guidance  and  Control  Directorate.  In 
addition,  summarized  within  the  report  are  the  presentations  pre¬ 
pared  by  Dr.  Seltzer  and  trips  undertaken  in  its  behalf. 


Thi»  document  ha*  beoa  approved 
for  public  reloaeo  and  sale;  its 
distribution  is  unlimited. 
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SECTION  I.  INTRODUCTION 


A,  Study  0b.1ecti.T8 

Hie  purpose  of  this  study  Is  to  investigate  the 
application  and  implementation  of  advanced  discrete  guidance  and 
control  theory  and  analysis  to  a  United  States  Any  Missile  System 
of  the  future*  An  aim  of  this  investigation  is  to  achieve  modularity 
of  the  guidance  and  control  system  functions.  The  end  result  of 
this  study  is  the  analysis  and  preliminary  design  of  a  discrete 
(or  digital)  guidance  and  control  system  (or  systems)  that  can  be 
used  for  a  U.  S.  Amy  Guided  Missile  System* 

The  study  effort  is  composed  of  two  tasks  described 
in  the  Scope  of  Work,  contained  in  Technical  Requirement  T-0107, 
entitled  "Application  of  Discrete  Guidance  and  Control  Theory. 

The  requirements  imposed  by  these  two  tasks  are  contained  in  the 
referenced  contract  document  and  re-eta ted  (for  convenience  of  the 
reader)  below. 

Task  I  requirements t 

2*1  Obtain  or  derive  several  candidate  guidance 
and  control  laws  applicable  to  boost,  coast,  and  terminal  flight 
phases  of  a  long  range  ground-to-ground  missile  system. 

2.2  Obtain  or  derive  mathematical  models  of  the 
missile  system  applicable  to  the  three  flight  phases  identified  in 
Subsection  1*1  above.  The  models  shall  include  significant  dynamics 
of  missile  plant,  sensors,  effectors,  on-board  computer  or  processor, 
and  external  disturbances  (such  as  aerodynamics)*  In  particular, 
identify  any  significant  nonlinearities  in  the  model.  Since  a  major 
portion  of  this  modeling  already  exists  in-house,  this  task  be 
oriented  to  extending  that  currently  available. 

2.3  Mathematically  cast  the  several  candidate 
guidance  and  control  laws.  Subsection  1.1  above,  in  a  form  amenable 
to  their  implementation  in  an  on-board  conputer  or  processor. 

Task  II  requirements : 

2  .4  Extend  the  guidance  and  control  laws  of 
Subsection  2.1  and  Subsection  2.3  via  the  parameter  space  technique 
far  the  analysis  and  design  of  digital  control  systems.  Investigate 
the  capabilities  and  limitations  of  this  technique  to  the  discrete 
control  system  of  the  missile  developed  in  Subsection  2.3  and  com¬ 
pare  the  results  obtained  through  the  use  of  extant  sampled -data 
techniques. 
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2.5  Confine  the  outputs  resulting  from  Sub¬ 
sections  2*2,  L.3,  and  i.f  above  so  the  dynamics  of  each  system 
may  be  determined  analytically,  using  both  existing  and  newly- 
developed  digital  and  sampled -data  techniques. 

2.6  Deleted  by  Contract  Amendment,  Modification 
No.  P00002,  dated  31  January  1979. 2 

?. 7.  Document  In  complete  detail  the  modeling, 
analysis,  and  simulation  results  described  in  Subsections  2.1 
through  2*C  above* 

Conduct  a  technical  seminar  for  selected 
members  of  the  Guidance  and  Control  Analysis  Group,  Guidance  and 
Control  Directorate  (DRDMI-TGN).2  Hie  purpose  of  the  seminar  is 
to  develop  a  design  expertise  among  members  of  the  Group  to  enable 
them  to  design  and  evaluate  discrete  guidance  and  control  systems* 
Both  existing  and  newly-developed  sampled -data  techniques  will  be 
described  and  discussed  in  detail  with  practical  applications  con¬ 
ducted  during  the  course  of  the  seminar  conference* 

B,  Study  Schedule 

The  requirements  of  Task  I  were  completed  during  the 
period,  called  Phase  I,  7  August  1978  (date  of  contract  initiation) 
through  30  September  1978.  They  are  -reported  upon  in  Applied  Control 
Systems  Dynamics  Company  Report  ACDS-78-2,  dated  September  1978, 
entitled:  "Application  of  Discrete  Guidance  and  Control  Theory. 
Technical  Report:  Phase  I* "3 

The  requirements  of  Task  II  are  to  be  completed  during 
the  period,  called  Phase  II,  that  culminates  in  twelve  (12)  months 
duration*  TM«  is  construed  to  mean  the  period  ending  6  August  1979. 
It  began  with  Amendment  1,  Modification  P00001,  effective  20  October 
1978.5 


C.  Contents  of  Report 

TMa  technical  report  covers  the  effort  performed  in 
completing  the  requirements  of  Task  II  of  the  study,  i.e».  Phase  U* 
As  stated  in  the  referenced  contract  requirements,  the  documentation 
of  this  technical  report  shall  Include,  but  not  be  limited  bo,* 

a*  A  stability  analysis,  summary  and  implementation 
of  the  discrete  control  for  the  missile  system  applicable  via  the 
parameter  space  technique.  The  results  and  comparisons  made  for  the 
extant  sampled -data  technique  shall  be  documented. 
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b.  Any  simulation  initializing  operations. 

c.  Any  simulation  input  parameters  required 

and  their  formats. 

d.  Any  analytic  results  produced. 

e.  Operating  instructions  far  a  user  of  the 
computer  simulation(s)  generated. 

f .  (deleted ) 

g.  Summary  of  computer  simulation  runs,  results, 
and  analyses  performed. 

Hie  documentation  shall  cover  work  done  to  perform 
Task  II  from  20  October  1978  (beginning  of  Modification  P00001  of 
the  contract)  until  completion  of  contract.  The  required  documenta¬ 
tion  shall  be  a  technical  report  to  be  delivered  to  the  Guidance  and 
Control  Analysis  Group  (DRDMI-TGN). 
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SECTION  H.  DEVELOPMENT  PROGRAM  FOR  ADVANCE)  GUIDANCE  &  CONTROL  SYSTEM 

The  definition  of  the  development  program  for  the  Advanced 
Guidance  and  Control  System  was  developed  during  Task  I  and  may  be 
found  in  the  referenced  technical  report :  Riase  1,3  The  original 
development  program,  then  entitled  "Development  Program  for  an  Optimal 
Digital  Guidance  and  Control  System  for  U„  S.  Amy  Future  Missiles" 
is  included  as  Appendix  A  of  that  report.  That  development  program 
c emprises  the  framework  within  which  the  study  is  being  performed. 

A  summary  follows  of  that  program. 

The  U.  S.  Army  -lissile  Research  and  Development  Cormand 
(M3RADC0M)  recently  embarked  on  a  task  to  develop  an  Advanced  Guidance 
and  Control  (G&C)  System  far  future  U,  S,  Army  Modular  Missiles.  The 
intent  is  to  "leapfrog"  systems  currently  under  development.  The  reason 
far  embarking  on  this  task  is  that  existing  weapon  systems  may  be  ser¬ 
iously  degraded  in  engagements  against  targets  with  predicted  character¬ 
istics  of  the  1990's  and  beyond  and  in  the  battlefield  environments  of 
the  time  frame.  It  has  been  established  that  the  guidance  laws  currently 
in  use  may  not  be  adequate  to  control  these  dynamically  elusive  threats. 
Thus,  it  is  projected  that  fundamental  advances  in  G&C  systems  theory 
are  required  to  enhance  the  effectiveness  of  future  weapon  systems. 
Additionally,  missile  airframe  and  propulsion  systems  may  require 
advances  cannons  urate  with  predicted  target  scenarios.  In  particular, 
air  defense  weapons  currently  in  R&D  may  be  seriously  hampered  in  the 
combat  scenarios  envisioned,  from  an  overall  system  viewpoint,  this 
task  shall  address  the  issue  of  creating  new  theory  in  the  G&C  area 
to  meet  the  high  performance  threats  of  the  future,  from  the  systems 
requirements,  the  objective  clearly  emerges  to  develop  and  validate 
research  that  will  enable  the  digital  optimal  guidance  and  control  of 
future  U.  S.  Army  Modular  Missiles .  The  work  that  has  taken  place  to 
meet  this  objective  began  in  fiscal  year  1979.  It  began  with  a  con¬ 
tinuing  definition  of  the  dynamics  of  predicted  future  targets.  This 
was  followed  by  performing  and  documenting  an  extensive  search  of 
literature  pertaining  to  the  theory  and  implementation  of  missile 
guidance  laws,  A  trade  study  was  conducted  to  compare  capabilities 
of  the  guidance  law  classes,  leading  to  the  selection  of  an  optimal 
formulation  of  the  advanced  guidance  law.  Development  began  of  analyti¬ 
cal  models  of  candidate  future  missile  and  subsystems  (to  enable  realistic 
G&C  system  research).  New  digital  control  theory  also  was  developed 
with  an  eye  toward  simplifying  elusive ly  difficult  mathematical  tools 
far  designing  digital  systems.  For  example,  an  alternative  to  the 
usually  difficult  and  cumbersome  digital  signal  flow  graph  method  of 
design  already  has  been  developed.  In  fiscal  year  1980.  it  is  planned 
to  further  extend  and  develop  new  digital  (sampled -data)  theory. 

Computer  software  will  be  developed  to  implement  these  (and  other) 
theoretical  advanced  analysis  tools.  Guidance  law  performance  indices 
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will  be  evaluated,  and  the  need  for  on-board  state  estimation  techniques 
will  be  determined*  Parallel  to  this  effort,  fluidics  technology  will 
continue  to  be  developed.  Beyond  fiscal  year  1980,  candidate  advanced 
G&C  systems  configurations  will  be  defined  and  a  figure  of  merit  developed* 
Advanced  analysis  and  simulation  techniques  will  be  used  to  evaluate 
the  candidate  G&C  systems.  Additionally,  the  merits  of  advanced  elect¬ 
ronics,  fluidics,  and  distributed  controllers  will  be  evaluated* 

The  foregoing  G&C  efforts  are  being  conducted  within  the  Technology 
Laboratory  of  MIRADCOM.  The  Advanced  G&C  System  program  is  defined 
within  the  G&C  Directorate.  It  is  being  implemented  by  both  in-house 
engineers  and  by  contracted  effort  at  universities  and  local  industry* 

It  is  to  support  this  Advanced  G&C  System  program  that  the  contract 
described  by  this  technical  report  was  made. 

The  participation  of  the  G&C  Analysis  Group  in  the 
Advanced  G&C  System  development  has  been  through  the  assignment  of 
particular  areas  of  research  and  development.  Hie  assignments  are 
indicated  in  Appendix  A  of  Technical  Report:  Phase  I  and  in  Appendices  I 
and  J  of  this  report.  Progress  has  been  indicated  by  individual  oral 
quarterly  reports  to  Hr*  Russell  T.  Gambill.  It  is  planned  to  formalize 
a  progress  report  at  the  end  of  fiscal  year  1979  by  the  Guidance  and 
Control  Analysis  Group. 
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SECTION  m.  SURVEY  OF  PERTMENT  STATE-OF-THE-ART 

As  reported  in  the  Technical  Report:  Phase  I,  the  U.  S. 
Air  Force  presently  is  investigating  and  developing  advanced  digital 
guidance  and  control  concepts  for  air-to-air  missiles.  Or.  Seltzer 
has  maintained  close  liaison  with  the  principal  investigators  at  the 
TJ.  S.  Air  Force  Armament  Laboratory,  Eglin  Air  Force  Base.  Of  parti¬ 
cular  interest  has  been  their  development  of  an  algorithm  to  predict 
"time-to-go"  (remaining  time  of  flight).  This  important  parameter  is 
needed  by  both  optimal  guidance  law  formulations  and  Disturbance 
Accommodating  Control  formulations.  It  presently  is  being  investigated 
for  applicability  and  ease  of  implementation. 

An  awareness  of  the  rapidly  changing  state-of-the-art 
has  been  maintained  by  Dr.  Seltzer,  partially  through  his  participation 
in  American  Institute  of  Aeronautics  and  Astronautics  (AIAA)  conferences 
pertinent  to  the  subject  at  hand  (see  Section  VII).  This  is  augmented 
by  extensive  selected  reading  of  pertinent  literature  and  technical 
publications. 


SECTION  17.  PERFORMANCE  REQUIREMENTS 

As  a  result  of  extensive  studies  of  predicted  future 
target  characteristics  conducted  by  personnel  of  the  Guidance  and 
Control  Analysis  Group,  Dr.  Seltzer  conducted  extensive  discussions 
and  analyses  of  the  gathered  information  with  these  personnel  and 
personnel  of  the  Department  of  Defense  (see  Section  XI).  The  results 
were  the  following: 

1.  The  surface  targets  of  the  1990’s  era  pro¬ 
bably  will  be  capable  of  being  engaged  by  surface-to-surface  missiles 
now  in  the  field  or  presently  under  development. 

2.  The  aerial  targets  of  the  1990’s  era  that 

may  be  too  elusive  dynamically  for  surface-to-air  missiles  presently 
under  development  are: 

a.  Tactical  ballistic  missiles; 

b.  Aircraft  (manned  or  unmanned); 

c.  Cruise  missiles;  and 

d.  Remote  pilotless  vehicles  (RPV's). 

The  predicted  large  number  of  RPV’s  probably  dictates  their  engage¬ 
ment  by  less  expensive  weapons  than  a  guided  missile.  Hence,  the 
first  three  predicted  future  target  classes  are  considered  to  be 
those  that  must  be  engaged  by  an  Advanced  Guidance  and  Control  System. 
Hence,  it  is  their  characteristics  that  are  used  to  define  the 
necessary  performance  characteristics  of  the  Advanced  G&C  System 
and  its  missile.  To  this  end,  a  design  point  for  the  proposed 
future  missile  is  developed  as  follows: 

Mach  No. :  li 
lateral  g's:  12 
Altitude:  70,000  feet. 

It  is  expected  that  these  values  will  be  exceeded  by  a  significant 
margin,  e.g.,  perhaps  20  g’s. 

In  an  attempt  to  determine  analytically  if  such  per¬ 
formance  requirements  are  feasible,  the  Aeroballistics  Directorate 
was  contacted.  Members  of  that  Directorate  are  now  investigating 
three  possible  candidate  missile  configurations  with  the  hope  of 
determining  one  (or  more)  that  can  meet  the  desired  characteristics. 
They  are  shown  in  Appendix  K. 
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SECTION  V.  GUIDANCE  AND  CONTROL  LAW  EXTENSION 


A.  Introduction 

Mathematical  models  of  candidate  guidance  laws 
were  categorized  and  described  in  Section  17  of  the  Technical 
Reports  Phase  1.3  A  means  of  extending  guidance  and  control 
laws  is  developed  in  Part  B  of  this  section.  It  is  then  applied  - 
via  the  parameter  space  technique  in  particular  -  to  a  selected 
model  of  a  future  missile  system,  using  as  a  basis  Section  17 
of  the  Technical  Report:  Phase  I.  This  is  performed  in  Part  C 
of  this  section.  The  results  are  then  compared  with  those  that 
would  be  obtained  by  extant  sampled -data  techniques. 

B.  Design  Tool  Development 

As  a  means  of  meeting  Section  2,h  of  the  referenced 
Scope  of  Work  in  Technical  Requirement  T-0107,  three  new  design  and 
analysis  tools  were  developed  by  Dr.  Seltzer,  They  are  described 
summarily  in  the  following  three  subparagraphs  and  In  detail  in 
Appendices  A,  B,  and  C  of  this  report. 

1.  SAM:  An  Alternative  to  Sampled -Data  Signal 

Flow  Graphs. 


Th  Appendix  A,  an  alternative  to  the  use  of 
Signal  Flow  Graphs  (SFG*s)  and  Mason's  Gain  Rule  (Formula)  is 
presented  to  use  in  the  analysis  of  complicated  sampled -data 
digital  control  systems.  Usually  in  such  systems,  block  diagram 
algebraic  manipulation  nay  become  unwield ly,  particularly  when 
such  systems  include  multiple  loops  and  samplers.  The  Systematic 
Analysis  Method  (SAM)  may  be  applied  to  such  systems  (as  well  as 
to  simple  single-loop  feedback  systems,  of  course).  This  is  shown 
in  Section  II  of  Appendix  A.  Also  shown  is  how  to  apply  SAM  to 
make  use  of  Modified  Z  -  transforms  (Section  7  of  Appendix  A). 

The  advantages  of  using  SAM  are  that  the 
cumbersome  application  of  Mason's  Gain  Formula  can  be  avoided. 
Further,  the  entire  method  of  drawing  Signal  Flow  Graphs  may 
be  circumvented.  Since  only  the  equations  describing  the  system 
are  needed  for  SAM,  even  the  customary  block  diagram  is  not  needed. 

If  the  analyst  prefers  to  use  one  of  the  Signal 
Flow  Graph  (SFG)  methods,  a  modified  SFG  technique  is  also  described 
(Section  in  of  Appendix  A).  It  is  simpler  and  less  cumbersome  to 
apply  than  the  conventional  SFG  method  which,  far  purposes  of  com¬ 
parison,  is  described  in  Section  17  of  Appendix  A. 


All  three  techniques  are  applied  to  two  examples 
in  Appendix  A,  This  Is  done  to  better  describe  the  application  of  SAH 
and  the  modified  SFG  and  to  help  provide  a  basis  for  the  following 
comparison  of  the  three  methods. 

The  Systematic  Analysis  Method  (SAH)  can  be  used 
to  determine  the  states  of  a  digital  system  in  terms  of  that  system’s 
transfer  functions  and  the  inputs  to  the  system.  This ,  of  course,  can 
also  be  done  by  the  application  of  other  methods,  such  as  signal  flow 
graph  (SFG)  techniques.  The  usual  advertised  advantages  of  the  latter, 
when  they  are  compared  to  block  diagram  or  algebraic  manipulation,  are 
that  they  are  particularly  amenable  to  the  analysis  of  complicated 
systems . 

It  is  demonstrated  in  Appendix  A  of  this  report 
that  SAH  can  handle  digital  system  analysis  as  capably  as  can  SFG 
methods.  It  has  the  advantage  of  not  requiring  the  cumbersome  Mason's 
Gain  Rule.  It  hence  avoids  the  oft-committed  errors  associated  with 
SFG  analysis;  such  as  over-looked  closed -loops,  non -obvious  forward 
paths,  etc..  As  with  SFG's,  a  block  diagram  is  not  needed  -  the 
system  equations  are  sufficient. 

Finally,  it  is  demonstrated  that  SAM  Is  easy  to 
implement.  It  appears  to  take  less  lengthy  analytical  manipulation. 

If  the  analyst  prefers  using  SFG’s  to  either  block 
diagrams  or  algebraic  manipulations,  a  modified  SFG  technique  is  pre¬ 
sented.  It  is  based  on  SAM  techniques.  As  such,  it  is  systematic. 

■While  the  SFG’s  produced  by  this  technique  are  usually  different  from 
those  produced  by  standard  SFG  techniques,  they  yield  the  same  results. 
Mason's  Gain  Rule  is  only  applied  at  one  stage  of  the  analysis  in  the 
modified  SFG  technique  -  as  opposed  to  the  standacd  SFG  technique  which 
requires  several  applications  of  Mason’s  Gain  Rule. 

In  conclusion  then,  an  alternative  to  the  Signal 
Flow  Graph  (SFG)  technique  has  been  presented  and  compared  to  a  standard 
and  a  modified  SFG.  Hie  alternative,  termed  Systematic  Analytical 
Method  or  SAM,  is  claimed  herein  to  be  simpler  and  more  straightforward 
to  implement  than  the  SFG  methods.  Not  only  does  it  appear  to  be  quicker 
to  use  In  system  analysis,  but  it  obviates  the  cumbersome  use  of  Mason's 
Gain  Rule.  Far  the  analyst  who  desires  to  use  SFG  techniques,  a  modified 
SFG  technique  is  presented.  It  is  systematic  and  reduces  the  number  of 
required  applications  of  Mason's  Gain  Rule. 


2.  Determination  of  Digital  Control  System  Response 
by  Cross-Multiplication. 

3h  Appendix  B,  a  means  is  described  for  obtaining 
the  response  of  a  digital  control  system  from  its  closed-loop  transfer 
function.  It  is  assumed  that  the  closed -loop  transfer  function  is 
available  in  the  z-  or  modified  z-transform  domain.  The  numerator 
and  denominator  of  each  side  of  the  transfer  function  are  cross -multi¬ 
plied.  The  Real  Translation  Theorem  is  then  applied  to  the  result, 
yielding  a  difference  equation  in  the  time-domain.  This  may  be 
solved  for  the  system  response  in  terms  of  the  reference  (or  other) 
input  (s)  to  the  system  as  well  as  in  terms  of  system  state  initial 
conditions. 


Two  different  modifications  to  the  basic  technique 
are  described:  one  using  the  submultiple  method  and  one  using  the 
modified  z -transform  technique.  These  are  applied  when  it  is  desired 
to  determine  the  intra-eampling  responses  of  the  system. 

All  three  techniques  are  applied  to  a  single  example. 

A  summary  of  the  techniques  and  their  application  is  provided  at  the 
conclusion  of  the  report. 

It  has  been  assumed  that  a  given  digital  or  sampled - 
data  system  can  be  described  by  a  closed -loop  transfer  function  that 
relates  a  single  controlled  output  of  the  system  to  a  single  reference 
Input.  If  there  is  more  than  one  input,  the  technique  can  also  be 
applied  to  the  resulting  sum  of  closed -loop  transfer  functions  relating 
the  controlled  output  to  each  of  the  inputs.  Although  the  report  refers 
only  to  a  single  controlled  output,  the  technique  can  be  applied  to  find 
any  system  state  if  it  is  related  to  the  Inputs  to  the  system  in  the  z- 
doskdn.  These  relationships  may  be  derived  by  using  any  of  the  standard 
techniques  (such  as  signal  flow  graphs)  or  by  the  newly  developed 
Systematic  Analytical  Method  (SAM)  technique. 

Several  analytical  techniques  for  obtaining  the  response 
of  a  digital  control  system  have  been  described.  They  are  based  on  a 
single  principle:  cross-multiplication  followed  by  applications  of  the 
real  translation  theorems.  Each  is  applied  to  a  single  exasple.  As  a 
starting  point  for  application  of  each  of  the  techniques,  it  is  required 
that  the  dynamics  of  the  digital  control  system  be  described  in  the  z- 
or  modified  z -domain. 

The  advantages  of  the  three  techniques  over  extant 
classical  methods  are: 

a.  The  response  may  be  obtained  for  any  deterministic 
reference  input  into  the  system  as  long  as  its  value  is  known  at  the 
sampling  instants.  It  need  not  be  described  by  a  differential  equation, 
and  the  z-transform  for  a  specific  reference  Input  need  not  be  determined 
before  Obtaining  an  expression  for  the  response. 
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b.  Initial  conditions  and  Instantaneous  changes 
can  be  accommodated  readily  by  the  equations  obtained  through  use  of 
the  cross-multiplication  methods. 

c.  The  difference  equations  obtained  are  parti¬ 
cularly  amenable  to  programming  on  a  desktop  calculator  or  digital 
computer. 

d.  A  detailed  knowledge  of  the  theory  underlying 
digital  or  sampled -data  control  systems  is  not  required  (althougi  it 
certainly  is  helpful)  by  the  analyst  In  order  to  apply  the  recipes 
described  herein. 

3.  Sampled -Data  Analysis  In  Parameter  Space 

In  essence,  the  parameter  space  technique  provides  the 
digital  guidance  and  control  system  designer  with  a  means  for  determining 
the  stability  and  dynamic  characteristics  of  a  digital  control  system 
In  terms  of  several  selected  system  parameters.  The  method  requires 
that  the  system  characteristic  equation  be  available  In  the  complex 
z -domain.  It  is  this  capability  that  makes  the  method  mare  powerful 
than  most  design  techniques  (which  describe  stability  In  terms  of  only 
one  variable  parameter  or  gain). 

The  parameter  space  method  provides  an  analytical  tool 
developed  for  use  In  control  system  analysis  and  synthesis.  Although 
not  necessary,  its  application  is  facilitated  by  augmenting  the  analyti¬ 
cal  results  with  graphical  portrayals  in  a  selected  multi -parameter  space. 
The  method  requires  that  the  control  system  be  described  by  a  character¬ 
istic  equation  which,  far  sampled -data  or  digital  systems,  may  be  expressed 
In  the  z -domain.  The  technique  is  based  on  the  analysis  and  synthesis 
methods  far  linear  and  nonlinear  control  system  design  which  are  amply 
described  in  Siljak’s  excellent  monograph  on  the  subject.6  Reference  7 
describes  the  application  of  the  technique  to  the  analysis  and  synthesis 
of  linear  sampled -data  control  systems. 

Ohce  the  system  characteristic  equation  has  been  obtained, 
the  parameter  plane  method  enables  the  designer  to  evaluate  graphically 
the  locations  of  roots  of  the  equation.  Hence,  he  may  design  the  con¬ 
trol  system  in  terms  of  the  chosen  performance  criteria;  e.g.,  absolute 
stability,  damping  ratio,  and  settling  time.  He  is  able  to  see  the  effect 
an  the  characteristic  equation  roots  of  changing  two  adjustable 
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parane tars.  Siljak  farther  simplified  the  design  procedure  by  intro¬ 
ducing  Chebyshev  functions  into  the  equations,  thereby  putting  them 
in  a  form  particularly  amenable  to  their  solution  by  a  digital  computer* 
The  method  has  been  extended  to  portray  the  effect  of  varying  the 
aampUng  period.  The  extended  method  permits  one  to  see  the  effect 
of  the  choice  of  values  assigned  to  the  sampling  period  on  absolute 
and  relative  stability.  Also,  the  recursive  formulas  shown  therein 
are  simpler  in  farm  than  the  Chebyshev  functions  of  Reference  7. 

The  resulting  formulation  is  deliberately  cast  in  a  form  that  makes 
it  particularly  amenable  to  solution  by  a  digital  computer  or  a  desk 
calculator,  again  emphasizing  the  interplay  between  analysis  and  com¬ 
puting  machines.  When  portrayed  graphically,  the  results  show  the 
dynamic  relation  between  the  selected  parameters  and  the  character¬ 
istic  equation  roots,  as  a  function  of  the  nondimensianal  independent 
argument,  Hence  one  readily  can  deduce  the  dynamic  effect  upon 

the  system  of  various  combinations  of  values  of  the  selected  para¬ 
meters  defining  the  parameter  space. 

The  history  of  the  continuous-time  domain 
version  of  the  parameter  plane  technique  is  well— i escribed  with  suit¬ 
able  references  in  Reference  6.  Briefly  summarizing  that  history,  in 
1876  I.  A.  Vishnegradsky  of  the  Leningrad  School  of  Theoretical  and 
Applied  Mechanics  developed  and  used  the  first  version  of  the  para¬ 
meter  plane  techniques  to  portray  system  stability  and  transient 
characteristics  of  a  third  order  system  on  a  two -parameter  plane. 

Ih  19U9  Professor  Yu  I.  Neimark  of  the  Russian  School  of  Automatic 
Control  generalized  Vishne grad sky’s  approach  to  permit  the  decom¬ 
position  of  a  two-parameter  domain  (0;  describing  an  nib  order  system 
into  stable  and  unstable  regions.  The  technique  was  called  D -decom¬ 
position.  During  the  period  1959  through  1966,  Professor  D.  Mitrovic, 
founder  of  a  Belgrade  group  of  automatic  control,  extended  the  method 
to  enable  the  analyst  to  relate  the  system's  variable  parameters  to 
the  system  response,  using  the  last  two  coefficients  of  an  nth  order 
characteristic  equation.  Beginning  in  196U,  Professor  D.  D.  Siljak, 
then  a  student  of  'H-trovic's  at  the  University  of  Belgrade,  generalized 
the  method  and  called  it  the  Parameter  Plane  method.  His  method  per¬ 
mitted  the  analyst  to  select  an  arbitrary  pair  of  characteristic  equa¬ 
tion  coefficients  (or  parameters  appearing  within  the  coefficients) 
and  portray  both  graphically  and  analytically  the  dependence  of  the 
system  response  upon  the  selected  parameters.  The  method  was  extended 
subsequently  by  Siljak  and  others  to  encompass  a  host  of  related  pro¬ 
blems.  In  1967  Professor  J.  George  modified  the  D-decomposition  method 
to  enable  the  portrayal  of  the  absolute  stability  region  In  a  multi¬ 
parameter  space  (George  also  showed  how  to  portray  contours  of  relative 
stability,  as  did  Siljak).  All  of  the  foregoing  work  is  carefully  and 
completely  referenced  within  Reference  6.  In  1966  and  subsequently. 


Seltzer  has  applied  the  parameter  space  method  to  the  design  of 
missile,  aircraft,  and  satellite  controllers  including  systems  con¬ 
taining  one  or  two  nonlinearities;  the  analysis  of  the  dynamic  effects 
of  the  nonlinear  "Solid  Friction"  (Dahl)  model  for  systems  with  ball 
bearings,  such  as  control  moment  gyroscopes  and  reaction  wheels;  and 
the  specification  by  the  system  designer  of  the  dynamic  structural 
flexibility  constraints  to  the  structural  designer.  Most  of  this  work 
has  appeared  in  the  technical  journals  of  the  IEEE,  the  AIAA,  the 
International  Journal  of  Control,  and  the  journal.  Computers  and 
Electrical  Engineering.  A  portion  of  the  history  that  has  not  been 
reported  upon  previously  (with  one  exception  to  be  noted)  is  the  con¬ 
trol  system  work  conducted  by  the  German  rocket  scientists  in  the  early 
19U0's  at  Peenemunde.  There,  Dr.  W.  Haeussermann  and  others  applied 
the  D -decomposition  technique  to  the  design  of  the  7-2  Rocket,  following 
Dr.  Haeussermann’s  earlier  (pre-World  War  U)  application  to  the  control 
of  an  underwater  torpedo.  This  work  was  not  published  in  the  open  liter¬ 
ature  because  of  national  security  constraints.  When  the  group  came  to 
the  United  States  to  work  with  the  Army  Ballistic  Missile  Agency  (first 
at  Fort  Bliss,  Texas,  then  at  Redstone  Arsenal,  Huntsville,  Alabama), 

Dr.  Haeussermann  and  his  associates  continued  to  apply  the  method  to 
U.  S.  Army  Missiles  (and  later  to  NASA  space  vehicles).  Again,  national 
security  (this  time,  another  nation)  precluded  publication  in  the  open 
literature  until  1957 0 

In  196U  Professor  Siljak  published  the  first 
application  of  the  parameter  plane  technique  to  sampled -data  systems.? 

As  mentioned  above,  this  was  extended  in  References  8  through  9.  In 
1971  Seltzer  presented  an  algorithm  for  systematically  solving  the  Popov 
Criterion  applied  to  sampled -data  systems.  Applications  of  these  sampled - 
data  parameter  space  techniques  are  found  in  References  8,  9,  and  11. 

In  summary,  an  analytical  method  far  portraying 
stability  regions  in  a  selected  parameter  space  is  described  for  a  digital 
system.  The  method  requires  that  the  system  characteristic  equation  be 
available  and  expressed  in  the  complex  2 -domain.  It  also  is  possible  to 
apply  pole  placement  to  obtain  desired  dynamic  characteristics  using  this 
modified  parameter  space  technique.  The  advantage  of  the  technique  over 
existing  classical  sampled -data  methods  is  that  the  stability  and  dynamic 
response  characteristics  are  expressed  in  terms  of  several  (rather  than 
merely  one)  selected  parameters.  Also,  the  sampling  period,  T,  need  not 
be  expressed  numerically  before  the  design  technique  begins,  giving  the 
system  designer  one  mare  degree  of  freedom. 


C.  Surrey  and  Implementation  of  Discrete  Control  for 
Missile  Systems  via  the  Parameter  Space, 

1.  Introduction 

In  this  section,  the  results  of  mathematical 
models  of  Sections  IV  and  V  of  the  Technical  Report:  Phase  I  are 
used  to  demonstrate  the  application  of  the  newly-developed  digital 
and  sampled -data  techniques*  A  simplified  planar  model  of  a  missile 
system  is  used  to  demonstrate  the  application  of  the  three  design 
and  analysis  techniques  described  in  Part  B.  above.  It  is  described 
by  the  block  diagram  of  Figure  1,  where  G(s;  is  used  to  portray  the 
transfer  function  describing  the  characteristics  of  the  plant.  For 
lucidity,  aerodynamics  are  neglected  and  the  coefficients  of  the 
equations  describing  the  dynamics  of  the  system  are  assumed  to  be 
time-invariant.  The  symbol  Gjj0(s)  is  used  to  represent  the  transfer 
function  associated  with  a  zero  order  hold  in  a  digital  computer. 
Trapezoidal  integration  is  used  to  approximate  the  mathematics  assoc¬ 
iated  with  taking  an  integral  of  a  state  and  is  represented  by  the 
transfer  function  (in  the  complex  z -domain) ,  B(z).  These  three  transfer 
functions  may  be  written  as 


G(s)  -  3.  , 
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The  onboard  digital  computer  sampling  period  is  represented  by  the 
syafcol  T,  and  symbols  s  and  s  represent  the  complex  Laplace  (s)  and 
z  variables,  where 

d-  sT 

»  *  .  (I7-U) 

(to  Figure  1,  control  gains  are  represented  by  symbols  a-| ,  a 0,  and  a  ,  • 
The  attitude  of  the  missile  is  represented  by  9,  and  the  overdot  repre¬ 
sents  the  time  derivative  of  the  symbol  beneath  it. 


The  approach  taken  in  this  section  is  to: 

a*  First  determine  the  system  closed -loop 

transfer  function. 

b.  Then  determine  conditions  to  ensure  system 

stability. 

c.  Finally,  determine  representative  system 
responses  to  given  Inputs  to  the  system. 

In  each  of  the  three  portions  to  the  approach,  the 
methods  developed  under  this  contract  'Hill  be  applied.  It  is  at  this 
point  that  one  sees  the  logical  need  far  the  development  of  three  tools. 
"SAIT"  was  developed  to  provide  a  better  means  far  obtaining  the  closed- 
loop  transfer  function.  The  parameter  space  technique  was  developed  to 
provide  an  efficient  means  for  determining  conditions  that  insure  sta¬ 
bility  (and  that  also  can  ensure  a  desired  system  response).  The  cross¬ 
multiplication  technique  was  developed  to  provide  a  relatively  simple 
means  far  determining  the  system  response.  It  is  stated  that  these 
three  techniques  appear  to  be  superior  to  known  extant  methods.  An 
indication  of  this  superiority  Hill  be  provided  summarily  within  this 
section  at  the  conclusion  of  each  of  the  above  three  steps. 

2.  System  Closed -Loop  Transfer  Function. 

The  system  closed -loop  transfer  function  will  now 
be  derived  by  applying  the  "SAIT"  technique  described  in  Appendix  A. 

The  System  Equations  (see  Step  No.  1  of  the  SAM  procedure)  may  be 
written  directly  upon  inspection  of  Figure  1. 


E  -  X,  -X  +  Xq 
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0  *  Xg  ■  Xy/S 
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The  crverbar  over  a  state  Indicates  that  the  state's  value  at  the 
beginning  of  arty  sampling  period  is  held  (at  that  value)  until  the 
beginning  of  the  sampling  period.  The  asterisk  following  a  state 
indicates  the  instantaneous  value  of  the  state  at  the  instant  of 
the  assumed  sampling  impulse,  i.e. ,  it  Indicates  the  pulsed  trans¬ 
form  of  the  variable.  All  states  are  considered  to  be  written  in 
the  Ieplace  domain,  l.e.,  the  usual  notation  Xj(s)  has  been  short¬ 
ened  to  Xj  for  convenience. 

The  Modified  System  Equations  (see  Step  No.  2  of 
the  SAM  procedure)  corresponding  to  the  above  System  Equations  may 
be  written  directly  from  the  System  Equations. 
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The  Poised  System  Equations  (see  Step  No.  3  of  the 
|  SAM  procedure)  corresponding  to  the  above  equations  nay  be  written 

|  directly  from  the  System  Equations. 
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For  convenience,  one  nay  define  G1  (s)  as 
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which  leads  to 

1 

G^U)  -  ^  |gi  (s^  -  (1-z-1 )  | - - 

1 

-  -isJ2T2(z4l)/(z-l)2 

(I7-15Z) 

m 

One  nay  define  G2(z)  as 
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The  trapezoidal  approximation  far  the  integration  operation  (see 
Section  7  of  the  Technical  Report:  Phase  I)  nay  be  written  in  the 


T 

i 

z -domain  as: 

D(z)  -  %T(z+1  )/(z-l ) 

(17-17) 
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Making  use  of  these  definitions,  one  nay  substitute  then  into  the 
Pulsed  System  Equations  to  obtain 
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and,  after  some  algebraic  manipulation,  one  obtains  the  closed-loop 
transfer  function. 
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If  one  defines  the  dimensionless  gains,  kp,  k^,  kd,  as 

1 

kp  -  ^cgT2  , 
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and 
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one  may  rewrite  Equation  (17-19)  as  2 

9(z)  (kp^z2  +  21^  -  (kp-ki)  cfj* 

X(z )  C.E,(a)  C.E.Cz) 

(17-23) 

where  the  denominator  C.E.(z)  of  the  closed-loop  transfer  function 
is  the  characteristic  equation  for  the  system,  i.e. , 

C.E.(z)  -  z3  +  (kp-Hc1-Hcd-3)z2  +  (21^-21^  +3)a  +  (-kp 

•1) 

(I7-2L) 
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The  characteristic  equation  will  be  used  In  the  sequel  to  determine 
conditions  to  ensure  system  stability.  The  closed-loop  transfer 
function  will  be  used  to  determine  the  response  of  the  system. 

The  question  is  now  asked :  "How  well  does  SAM 
compare  to  other  methods  for  determining  the  closed -loop  transfer 
function?"  Several  means  will  be  examined  briefly  far  comparative 
purposes.  They  are  (1)  the  well-known  and  oft-used  Signal  Flow 
Graph  (SFG)  method:  and  (2)  block  diagram  manipulation. 

Even  if  the  modifications  put  forth  in  Section  III 
of  Appendix  A  are  implemented,  the  SFG  method  is  tedious  and  amenable 
to  mistakes  (such  as  not  recognizing  the  existence  of  loops  on  the 
graph).  Using  the  System  Equations  (17-5)  through  (I7-lh)  and  the 
Pulsed  System  Equations  (I7-5P)  through  (I7-1UP),  one  may  construct 
the  modified  SFG  found  in  Figure  2.  Mason's  Gain  Rule  (also  found 
in  Appendix  A)  is  then  applied  as  follows. 

There  are  three  loops  on  the  SFG. 

Loop  No.  1:  A  closed-loop  beginning  and  ending 
at  Xy*  with  a  gain  of  ^  ■  a-j  G2*. 

Loop  No.  2:  A  loop  passing  through  Xj*,  Xy*,  Xg*, 
and  back  to  X-|*  with  a  gain  of  K2  ■  a|G|*G2*(ao  +  a_iD*). 

Loop  No.  3*  A  loop  passing  through  X-j*,  Xg*, 
and  back  to  X1*  with  a  gain  of  ■  G1*(a0  +  D*) . 

It  is  observed  that  Loops  No.  1  and  No.  3  do  not  touch  each  other. 

The  value  of  A  (actually,  the  characteristic  equation)  is  found  to 
be 

A  .  1  _  k,  -  X2  -  k3  +  k1k3 

-  1  -  a1G2*  -  (sq  +  a^D*)  G,*  (a^*  +  l) 

+  a1G1^KJ2*  (sq  +  a^D*) 

-  1  -  a,G2*  -  (a0  +  a^D*)  G.,*  (17-25) 

This  may  be  shown  to  be  Identical  to  Equation  (I7-2lt)  which  is  the 
system  characteristic  equation  found  by  applying  SAM. 


I 

t 

I 

Returning  to  the  SFG,  one  sees  that  there  are  two 

I  forward  paths  from  "X*  to  9*.  Along  the  first  forward  path  (the 

upper  of  the  two), 

M,  *  a1G1^2*  (ag  +  a^D*)  (17-26) 

^  and 

A1  -  1  .  (17-27) 

I  Along  the  second  forward  path  (the  lower  one). 


&,*  (a0  +  a_-,D*) 


(17-28) 


and 


*2  “  1  - 


(17-29) 


The  expression  far  the  closed -loop  gain  between  1*  and  0*  is 

«i  A  *  ''A 

A 


9* 

TC* 


a1G1^-G2*  (aQ  +  a_.,I>*)  +  0,*  (aQ  +  a^G*)  (1  -  a^*) 


G^*  (ag  +  a_^D*) 


1  -  a,G2*  -  (aQ  +  «”■,&*)"  Q,* 


(17-30) 


It  is  seen  that  Equation  {17-30)  may  be  manipulated  so  that  it  is 
identical  to  Equation  (17-19)  obtained  by  using  SAM. 

The  same  result  may  be  obtained  by  block  diagram 
and/or  algebraic  manipulation.  Since  this  technique  is  well  known, 
it  will  not  be  described  here.  Suffice  it  to  say,  it  can  become 
tedious. 

The  advantages  that  SAM  appears  to  possess  are  that 
it  can  be  applied  to  a  system  having  any  level  of  complexity  or  order. 
It  is  not  subject  to  not  recognizing  a  path  or  closed -loop  (not  always 
obvious)  on  a  SFG.  3ecause  it  is  systematic,  it  is  less  amenable  to 
errors  by  the  user. 
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3.  System  Stabilit 


Hie  conditions  for  stability  of  the  system  are  now 
determined  by  using  the  parameter  space  technique  described  in 
Append  ix  C.  It  will  be  applied  to  the  characteristic  equation  ( 17-21*) • 
Since  it  is  a  third  order  polynomial,  it  has  three  roots  -  either  all 
three  are  real  or  one  is  real  and  two  are  complex  conjugates. 


The  z«1  boundary  is  found  by  substituting  the  value  of 
unity  for  z  in  Equation  (17-21*),  yielding 


k*  »  0 


(17-31 


The  z*-1  boundary  is  found  by  letting  z=-l  in 
Equation  (17-21*),  yielding 


kd  -2 


(17-32) 


The  complex  conjugate  root  boundary  (associated  with 
the  unit  circle  in  the  z -domain)  is  found  by  substituting  the  value. 


z^  =  R1  +  i  1^  (1  -  B2)^ 


(17-33) 


for  each  value  of  z^  (where  j  is  an  integer)  appearing  in  Equation  (17-21*), 
where  values  of  Rj  and  I.  corresponding  to  the  unit  circle  boundary  may 
be  found  from  Table  1  (see  Appendix  C  for  how  these  values  may  be  obtained) 
and  3  is  defined  as  coa/T.  (Recall  s  may  be  defined  as  s=*  o'  +  i«). 


TABLE  1 . 

7alues  of  Rj,  Ij.,  and  ^  j 


R3 

T3 

1 

0 

-kp  +  ki  +  -1 

3 

1 

2ki  — 2kjj  +  3 

2B2-1 

2B 

kp  +  ki  +  kd  -3 

1lB3-33 

l*B2-1 

1 

■i 


The  result  is  an  equation  containing  both  real  and  imaginary  quan¬ 
tities.  If  the  reals  and  the  imaginaries  are  separated  into  two 
equations,  one  obtains  the  following  two  simultaneous  algebraic 
equations : 


R  «  f  ^2  +  f  2^2  +  ^  1  ^  0^0  *  ® 


(I7-3J4H) 


i  -  r3i3  +  +  xi  +  a  0 


(I7-3UI) 


Substituting  the  values  from  Table  1  into  Equations  (I7-3U),  one  may 
obtain 

R  =  (B+1  )(B-1  )kp  +  B(B+1  )k±  +  B(B-1  )kd  +  (B-1  )2(2B+1  )  =  0  (I7-35R) 


I  ■  Bkp  +  (3+1  )k£  +  (B-I)k^  +  (B-1)(2B-1)  -  0 


(I7-35D 


Since  there  are  two  equations,  they  may  be  solved  for  any  two  para¬ 
meters  contained  within.  The  two  control  gains,  kn  and  Iqj,  may  be 
the  selected  two.  If  one  uses  Cramer’s  Rule  to  solve  for  kp  and  ky, 
one  may  obtain  the  Jacobian, 


(B+1 )  (B-1 )  B(B-1 ) 

B  (B-1 ) 


1  -B?0YRO 


(17-36) 


The  parameters  (control  gains)  may  be  solved  for,  obtaining 
-3(B+1  )ki  -(B-1  )2  (2B+1 )  B(B— 1  ) 

-.(B+lJk-  -  (B-1  )(2B-1 )  (B-1) 


(B+1 ) (B-1  )  -B(B+1  -  (B-1  )2(2B+1 ) 

B  -  (B+Olq  -  (B-1 )  (2B-1 ) 


A  summary  of  the  three  stability  boundaries  in  the  k_  ,  kj ,  kd 
parameter  space  is:  p 


2*+1  :  k^  ■  0  , 
z«-1 :  kd  ■  2  , 


z-e 


i«T. 


1-B  , 

k±  +  (1  -B)  . 


Although  these  may  readily  be  shown  in  three-space  with  proper 
graphics,  they  are  shown  in  this  report  in  a  k_,  kd  two-space  (or 
plane)  on  Figure  3  for  several  values  of  k±  (i.e.,  k^  -  0,  0.1,  and 
0.5).  The  shading  rules  described  In  Appendix  C  are  applied  to 
determine  the  stable  and  unstable  regions.  The  number  of  stable 
roots  in  a  particular  region  are  denoted  by  an  encircled  number. 
The  program  for  obtaining  numerical  values  for  generating  these 
curves,  using  a  Hewlett-Packard  9100B  calculator,  is  found  in 
Appendix  D. 


If  the  root  locus  technique  were  used  to  plot 
the  dynamic  characteristics  of  the  system,  this  could  be  accom¬ 
plished  only  after  two  of  the  three  parameters  were  assigned 
numerical  values,  thus  losing  a  considerable  amount  of  potential 
information.  Also,  the  numerator  and  denominator  of  the  root  locus 
equation  would  have  to  be  factored  into  first  and  second  order  terms 
to  determine  locations  of  poles  and  zeros.  Far  example,  if  k.  were 
selected  as  the  variable  parameter  (gain),  the  root  locus  equation 
(obtained  tram  the  characteristic  equation)  would  be: 


G(z)  - 


(z2  -2a  +  1 )  kd 

z3  +  (kp  *  k±  -3)z2  +  (2k±  +3)z  +(k^-Hc1-l) 


(17-39) 


The  corresponding  root  locus  equation  far  the  selection  of  kp  is 
equally  cumbersome: 


G(z) 


(z2  -2z  -H  )  kp _ 

z3  ♦  (Vkd-3)z2  +  (2ki-2kd+3)z  +  (kj^-l) 


(I7-U0) 
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of  several  generic  (rather  than  numerical)  values  of  gains.  In  this 
study,  the  dynamic  characteristics  of  the  system  have  not  been  speci¬ 
fied,  although  Appendix  C  tells  how  this  can  be  done  (essentially 
using  pole  placement  techniQue).  In  that  case,  the  parameter  space 
technique  provides  additional  strength,  since  a  numerical  value  of 
the  sampling  period  (T)  need  not  be  specified  before  design  commences. 

If  one  wishes  to  resort  to  array  techniques  (i.e,, 
the  sampled -data  analogue  of  the  Routh-Hurwitz  Criterion),  one  may 
use  the  Schur-Cohn  Method,  the  more  versatile  Jury's  Test,  or  the 
considerably  mare  versatile  Raible's  Test.  Even  in  the  case  of 
the  latter,  the  use  of  the  table  (when  implemented  with  a  digital 
calculator  or  by  hand  analysis)  is  cumbersome,  as  shown  in  Table  2. 

TABLE  2. 

Raible's  Test 


a3 

*2 

a1 

a0 

^a  "  ac/a3 

aoKa 

a1ka 

a2ka 

- 

b0 

n 

”5 

- 

kjj  -  b2Ao 

Vb 

D 

c0 

ci 

■B 

kg  -  ct/cq 

Cl^c 

mm 

d0 

- 

The  values  of  the  coefficients  in  Table  2  are  found  in  Table  3. 


!  “ 


J 


I 


TABLE  3. 

Ralble^  Test  Coeffieiants 


LU 

aj 

bJ 

ci 

d3 

B 

-kp+l^d-1 

a3“a0^a 

b0"b2kb 

Cg-C-jkc 

B 

21ci-2lcd+3 

^"Va 

bTb1kb 

- 

2 

VW3 

- 

- 

3 

1 

- 

- 

- 

To  use  the  Raible's  Array  of  Table  2  for  determining  the  conditions 
that  ensure  stability,  one  determines  the  values  of  parameters  that 
cause  bQ,  Cq,  and  dg  to  always  be  greater  than  unity.  While  an  the 
surface  this  appears  to  be  simple,  the  actual  implementation  requires 
great  algebraic  effort  (unless  numerical  values  -  with  their 
attendant  decrease  in  generality  -  are  used). 

U.  System  response 

The  system  response  to  an  input  is  now  determined, 
using  the  cross-multiplication  method  described  in  Appendix  B.  It 
will  be  applied  to  the  closed-loop  transfer  function  of  Equations  (17-23) 
through  (I7-2U).  Cross-multiplication  of  the  terms  in  Equation  (17-23) 
leads  directly  to  the  equation, 

z^0(z)  +  y2z20(z)  +i'1ze(z)  +fg0(z)  ** 


-52z2^(z)  -<£|  zX(z)  -<TgX(z)  .  (17-ljl) 

If  one  solves  Equation  (17-Ul)  for  the  value  of  0(z)  that  is  the 
coefficient  of  the  highest  power  of  z  (in  this  case,  three),  one 
obtains 

0(z)  -  -  [f2z-1X(z)+i1z-2X(z)+50z-3x(z)] 

-  [j2z"Vz)+f1t“Ze(zW0s“30(z)]  .  (I7-U2) 


25 


Equation  (I7-Uk)  nay  be  put  in  the  form  of  a  difference  equation, 
where  the  symbol  9v  denotes  the  value  of  0(t)  at  the  instant,  t*kT, 
k  an  integer,  and  denotes  the  value  of  X(t)  at  time  t-kT: 

®k  *  "[<$2^-1  +<5l*k-2  "^0^-3] 

"  2^-1  +*iek-1  +&09k-3"|  .  (I7-U3) 

The  power  of  the  difference  equation  is  respectable  in  the  era  of 
digital  machinery.  For  example,  in  the  case  of  Equation  (I7-li3), 
any  value  of  6(kT)  may  be  found  if  only  the  next  earlier  values  of 
9  are  known.  This  may  be  done  for  any  reference  input  signal,  "X(t), 
whose  value  is  known  at  the  sampling  instants  t-kT.  In  other  words, 
contrary  to  many  existing  methods  for  determining  system  response, 

"X(t)  need  not  be  specified  in  mathematically  closed-form.  Let  us 
dwell  on  this  power  briefly.  Suppose  the  initial  condition  for  the 
missile  system  under  study  herein  is  denoted  as  0Q,  i.e.,  the  value 
of  the  missile  attitude  at  the  instant  of  time  t-o.  Using  Equation  (I7-U3) 
and  assuming  that  t-0  was  selected  (as  it  surely  can  be)  at  an  instant 
before  which  the  state  8  is  at  rest,  or  zero  ,  one  obtains  the  value 
of  6  at  the  next  sampling  instant,  t*T: 

®1  *^2^0  "^2®0  •  (IF-UU) 

Knowing  9.,  one  may  again  use  Equation  (I7-U2),  this  time  to  find  the 
value  of  6  at  time  t*2T,  i.e., 

02  -<f2X,  +<^*0  -  (I r2a,  •  (I7-li5) 

Similarly,  values  of  8  may  be  determined  at  later  sampling  instants 
of  time  in  this  recursive  manner.  This  farm  obviously  is  amenable 
to  implementation  in  a  digital  computer  or  calculator. 

As  a  simple  example  of  the  application  of  the  technique, 
assume  the  system  is  commanded  to  follow  a  unit  step  input  applied  at 
time  t®0,  i.e.. 


T(t)  -  1  HiO 

-  0  V  t  <  0  .  (I7-U6) 
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The  response  to  the  step  input  may  be  determined  readily  from 


Equation  (iV-ljj)  as: 

(k-0)  0O  -  0  (I7-U7) 

(k»1 )  0-|  *  6  2*0  ■  kp  +  ki  (I7-U8) 

(k=2)  02  =  2+  (f ,  -r29i  “  kp  +  3ki  -  r26 1  (I V-U9) 

(k»3)  0^  *  (S2  +<f-|  *6  q)  “  (i*202  +  0-|)  ■  Ijkj  -  OC202  +)Ti9i)  (17-50) 

(k&j)  Gfc  -  1%  -  (fgOk.-i  +  JT 1 0k-2  +y09k-3)  •  (17-51) 

Equations  (I7-U7)  through  (17-51)  may  be  programmed  on  any  digital 


machine.  A  program  for  calculating  values  of  0(kT)  is  found  in 
Appendix  E. 


As  an  independent  check  on  the  response  values 
obtained,  the  initial  value  and  steady-state  value  of  0(kT)  may  be 
determined  by  applying  the  sampled -data  Initial  Value  and  Final  Value 
Theorems  to  the  closed-loop  transfer  function  of  Equation  ( 17-23). 


Initial  Value  Theorem: 

lim  0*(t)  »  lim  0(z)  ■  lim  (  6  2z2  +  <So^z) 

t-^0  z z+~  z3  +y2z2  +  JC-j2  +  8*0  # 


(17-52) 


The  answer  is  dependent  on  the  farm  of  VC(z).  For  a  unit  step  function. 
Equation  (I7-U6)  may  first  be  transformed  into  the  z -domain,  yielding 


Vt(z)  ■ 


(17-53) 


Because  of  the  nature  of  the  feedback  signal  in  the  missile  case  (see 
Figure  1),  the  value  of  TC(z)  should  be  assigned  a  negative  sign.  If 
the  negative  of  Equation  ( 17-53)  is  inserted  in  Equation  (17-52),  the 
initial,  value  of  0*(t)  is  seen  to  be  zero. 


lim  9*(t) 
t  <» 


lim  (1-z-1)  0(z) 
z->  1 


Im  -  f(kp-Hc1)z2  +  2k1z  -  (kp-kj) 

s3  +y2z2  +  l*i  a  +J*o 


lim 
z 


X(z)  (I7-5H) 
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Again,  if  the  theorem  is  applied  to  the  case  of  a  unit  step  input 
to  the  system  and  the  negative  of  Equation  (IV-53)  is  substituted 
into  Equation  (IV-5U),  the  resulting  steady-state  value  of  0*(t) 
becomes 

lim  6*(t)  -  *  *1  *  -  fcp  +Jfi  u  #  (17-55) 

t -roo  i  +]f2  *^0 

The  program  found  in  Appendix  E  has  been  applied  to  the  design 
example,  assuming  a  unit  step  input*  The  result  has  been  plotted 
as  Figure  U. 

The  brief  discussion  found  in  Appendix  B  of  other 
techniques  that  may  be  used  to  determine  the  system  response  is 
sufficient  to  convince  the  reader  of  the  advantages  of  the  cross¬ 
multiplication  technique*  Renee,  they  sill  not  be  applied  to  the 
design  example. 

i 

4a» 

1 

l 
1 
1 

i 
l 

1 
1 
1 
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SECTION  71.  TECHNICAL  SEMINAR 

A  technical  seminar  was  conducted  by  Hr.  Seltzer  for 
selected  members  of  the  Guidance  and  Control  Analysis  Group,  Guidance 
and  Control  Directorate  (DRDMI-TGN).  The  purpose  of  the  seminar  was 
to  develop  a  design  expertise  among  members  of  the  Group  to  enable  then 
to  design  and  evaluate  discrete  (i.e.,  digital)  guidance  and  control 
systems.  Both  existing  and  newly-developed  sampled -data  techniques 
were  described  in  detail  with  practical  applications  conducted  during 
the  course  of  the  seminar  conferences.  In  particular,  the  new  tech¬ 
niques  described  in  Appendices  A,  B,  and  C  of  this  report  were  fully 
developed  for  the  participants.  The  duration  of  the  seminar  was  123 
contact  hours. 

The  contents  of  the  seminar  are  described  in  the  following 
outline  format. 

A.  Introduction. 

B.  The  sampling  process. 

1.  Mathematical  description  of  the  sampling  process 
in  both  the  time  and  the  frequency  domains. 

2.  The  Sampling  ("Shannon's")  Theorem. 

C.  The  ideal  sampler:  the  impulse  sampling  approximation. 

1.  Motivation, 

2 .  Mathematical  description. 

a.  Time  domain. 

b.  Frequency  domain. 

(1)  Fourier  analysis, 

(2)  Laplace  analysis. 

(3)  Complex  convolution, 

(a)  Closed  form, 

(b)  Series  form, 

D.  Comparison  with  carrier-modulated  systems. 


z -plane. 


E,  Z -transform  analysis  of  linear  sampled-data  systems. 

1.  Introduction  to,  and  derivation  of,  the  z -transform, 

2.  Evaluation  of  z -transforms. 

3.  Happing  from  the  complex  s -plane  to  the  complex 

U.  Theorems  of  the  z -transform. 

£.  The  concept  of  pulse  transfer  functions* 


cascaded  elements* 


6.  Z -domain  transfer  function  relationships  between 


7.  limitations  of  the  z -transform  method. 


8.  Determination  of  the  control  system  response 
between  sampled  instants* 

a*  Use  of  the  submultiple  method* 

b.  Use  of  the  Modified  and  the  Delayed  z -transforms. 

9*  Modified  z -transform  theorems. 

E.  Data  reconstruction,  emphasizing  zero-order  and  first- 

order  holds* 

G.  Sampled-data  systems:  Determination  of  outputs  and 
other  states  of  digital  control  systems* 


pulation* 


1.  Block  diagram  manipulation  and  algebraic  rnani- 


2*  Signal  Flow  Graphs  and  'las on's  Gain  Formula, 
a*  Conventional  (B.C.  Kuo)  technique, 
b.  Seltzer's  technique. 

3.  A  newly-developed  systematic  algebraic  mani¬ 
pulation  technique  termed  "SAM"  (Systematic  Algebraic  Manipulation) . 


H.  System  Response. 

1.  Comparison  of  second-order  response  for: 

a.  Continuous  system. 

b.  Sampled -data  system  without  data  holds. 

c.  Sampled -data  system  with  zero-order  hold. 

2.  Stability  and  other  dynamic  determinations. 

a.  Jury’s  Test, 

b.  Raible’s  Test, 

c.  lapping  from  the  complex  s -plane  to  the 

complex  v-  and  r-planes. 

d.  Extension  of  the  Routh-Hurwitz  Test  to  the 

w-  and  r -domains. 

e.  Effect  of  pole-zero  locations  an  system 
dynamic  response  (in  the  z -domain). 

f.  The  root  locus  method  (in  the  z-domain) 
for  sampled -data  systems. 

g.  The  Nyqtdst  Criterion. 

h.  The  Bode  diagram. 

i.  The  gain-phase  plot  (Nichols’  Chart). 

j.  The  parameter  space. 

3.  Response  between  sampling  instants. 

a.  Submultiple  sampling  method. 

b.  Modified  z -transform  method, 

c.  The  newly-developed  cross-multiplication 

method. 


d.  Hidden  Instabilities 


SECTION  VH.  AMERICAN  INSTITDTS  CF  AERONAUTICS  AND  ASTRONAUTICS 

(AIAA)  ACTIVITIES 

During  the  contract  period.  Dr.  Seltzer  has  supported 
certain  pertinent  AIAA  activities.  This  has  taken  the  following  form. 

A.  AIAA  Aerospace  Sciences  Meeting 

Dr.  Seltzer  was  selected  to  organize  the  Guidance 
and  Control  (G&C)  sessions  at  the  AIAA  Aerospace  Sciences  Meeting 
held  in  New  Orleans,  Louisiana,  15-17  January  1979.  Also,  at  that 
meeting.  Dr.  Seltzer  presented  a  paper  co-authored  by  Dr.  Harold  L. 

Pas trick  (Guidance  and  Control  Directorate,  U.  S.  Army  Missile 
Research  and  Development  Command)  &  Professor  Michael  Warren 
(University  of  Florida  in  Gainesville)..  The  paper  was  entitled 
"Guidance  Laws  for  Tactical  Missiles".'3  It  includes  a  comparison 
of  guidance  laws  applicable  to  short  range  tactical  missiles. 

These  lavs  are  segmented  into  several  classes  and  the  principles 
underlying  each  class  are  discussed.  Specific  attention  is  given 
to  the  structure  of  the  guidance  technique  and  the  requirements  for 
its  implementation.  Evaluation  and  comparison  of  the  performance 
of  each  guidance  law  versus  the  cost  of  implementing  it  are  con¬ 
sidered.  An  extensive  bibliography  of  relevant  literature  is 
included.  A  copy  of  the  paper  is  included  as  Appendix  F. 

B.  AIAA  Guidance  and  Control  Conference 

A  paper  co-authored  by  Dr.  Pas  trick  and  Dr.  Seltzer 
has  been  accepted  for  delivery  at  the  AIAA  Guidance  and  Control  Con¬ 
ference  to  be  held  6-6  August  1979.  The  paper  is  entitled  "Future 
U.  S.  Amy  Missile  Guidance  and  Control  Systems".  It  described  the 
task  that  was  recently  begun  by  MIRADCOM  to  develop  an  advanced  G&C 
system  for  Mature  Amy  Modular  Missiles  (see  Section  II  of  this 
report).  The  paper  defines  the  problem  and  described  the  approach 
being  taken  as  part  of  a  seven-year  program  to  develop  a  viable 
future  C&C  system.  The  preliminary  results  that  have  been  achieved, 
as  well  as  projected  results,  are  to  be  presented . 
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SECTION  VIII.  SIG-D  EVALUATION 


On  19  January  1979,  Ur.  J.  B.  Huff,  Director  of  the 
Guidance  and  Control  Directorate  appointed  a  team  to  conduct  an  in- 
depth  review  of  the  digitally-controlled  SIG-D  program  (Appendix  G). 
The  team  was  placed  under  the  leadership  of  Hr.  Russell  T.  Gambill. 
Because  of  his  expertise  in  the  digital  control  field.  Dr.  Seltzer 
was  designated  to  evaluate  the  theory  and  design  portion  of  the  SIG-D 
program.  The  formal  review  took  place  during  8-9  February  1979. 
T.'Mtten  reports  were  submitted  by  each  of  the  ten  team  members  to 
Hr.  Gambill  who  combined  them  into  an  overall  report  entitled 
"SIG-D  Red  Team  Report?  dated  21  February  1979.  Section  II  of  the 
report  entitled  "Theory  and  Design"  was  prepared  and  submitted  by 
Dr.  Seltzer.  A  copy  is  included  in  Appendix  G. 


SECTION  n,  HELLF3RE  T.flIITE  PAPER 


A  "white  paper"  was  prepared  by  Dr.  Seltzer  describing 
how  the  Guidance  and  Control  (G&C)  Analysis  Group's  Hardware-in-the- 
loop  (HWXL)  simulation  facility  is  u-  *d  to  meet  that  Group's  HELLE3RE 
mission.  In  particular,  the  maimer  in  tdiich  the  G&C  Analysis  Group 
met  the  objective  outlined  in  the  "Lew  Cost  Laser  Seeker  Evaluation 
Test  Plan"  was  addressed.  The  paper  was  generated  to  meet  the  sharp 
criticism  leveled  against  the  HWH  facility  by  Mr.  Art  Kendell, 
Rockwell  International  -  Columbus.  A  copy  of  the  paper  is  included 
in  Appendix  H. 


SECTION  X.  PRESENTATIONS 

Several  presentations  were  made  during  the  course  of 
the  contract  period  to  demonstrate  progress  on  tasks  associated 
with  the  contract, 

A,  Presentation  to  Dr,  J,  B,  Huff  and  Dr,  Richard  Hartman, 

On  10  May  1979,  a  presentation  was  made  to  Dr.  J.  B, 

Huff,  Director  of  the  Guidance  and  Control  Directorate,  and  Dr,  Richard 
Hartman,  Director  of  the  Research  Directorate,  and  members  of  their 
organizations.  The  purpose  of  the  briefing,  held  at  Dr.  Huff’s  request, 
was  to  assess  the  program  status  of  the  research  on  the  Future  U,  S, 

-Army  Modular  Missile  "Advanced  Guidance  and  Control"  and  to  describe 
the  program  plans.  To  meet  this  request.  Dr.  Seltzer  organized  the 
presentation  to  cover  the  following  topics. 

a.  Program  objectives 

b.  Program  plan 

c.  Initial  achievements  and  status 

d.  Projected  accomplishments 

e.  Program  implementation, 

A  copy  of  the  vu-graphs  used  in  the  presentation  is  included  in  Appendix  I, 

B.  Presentation  to  Advanced  Sensors  Directorate, 

As  a  result  of  the  above  presentation,  IT.  Hartman  set 
up  a  presentation  for  members  of  the  Advanced  Sensors  Directorate  on 
13  June  1979,  Hie  purpose  of  the  presentation  was  to  describe  the 
Advanced  Guidance  and  Control  program,  to  describe  the  threats  that 
the  program  addresses,  and  to  enlist  the  assistance  of  Advanced  Sensors 
personnel  in  the  program.  Bie  description  of  the  program  was  presented 
by  Dr.  Seltzer  -  and  the  description  of  the  threats  and  their  analysis 
was  made  by  Mr,  Ray  'aples.  Copies  of  the  vu-graphs  used  in  the  pre¬ 
sentation  are  included  in  Appendix  J. 
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SECTION  XI.  TRIPS  IN  SUPPORT  OF  TEE  CONTRACT 

Several  trips  were  mode  in  support  of  the  contract 
during  the  implementation  of  Task  II. 

A.  Trip  to  Eglin  Air  Force  Base. 

Qa  12  October  1978,  Dr.  Seltzer  joined  Mr,  Russell 
T.  Gambill  and  Dr.  Harold  L.  Pas  trick  at  Egiin  Air  Force  Base. 

Their  major  contact  and  source  of  technical  and  program  information 
was  2nd  Lt.  Tom  Riggs,  U.  S.  Air  Force,  Air-to-Air,  Systems  Analysis 
and  Simulation  Branch,  Air  Farce  Armament  Laboratory,  Air  Farce  Systems 
Coctnand,  Eglin  Air  Force  Base,  Florida.  Information  was  exchanged 
between  the  principals  mentioned  in  considerable  detail.  Continued 
exchanges  of  information  were  sought  by  the  Air  Force  personnel  present. 

3h  addition  to  considerable  in-house  effort,  Eglin 
Air  Force  Base  has  contracts  with  three  universities  (Alabama,  Florida, 
and  Texas)  and  three  members  of  industry  (OrinCon,  SCI,  and  TASC)  to 
study  advanced  digital  G&C  for  air-to-air  tactical  missiles. 

The  Eglin  personnel  expressed  a  desire  to  develop  a 
memorandum  of  understanding  with  the  U.  S.  Airy  IQRACOM  to  formalize 
a  joint  effort  in  the  development  of  advanced  digital  G&C  systems. 

Dr.  Seltzer  was  requested  to  write  a  first  draft  of  the  agreement. 

This  was  presented  to  Lt.  Riggs  during  a  visit  to  IHRADCOM  on 
21  March  1979.  It  presently  is  under  review  at  Eglin. 

B.  Trip  to  the  Pentagon. 

On  22-23  April  1979,  Dr.  Seltzer  visited  the  Director 
of  Land  '-Jar fare,  Mr.  Charles  Bernard,  of  the  Office  of  the  Under 
Secretary  of  Defense  for  Research  and  Development,  and  members  of  his 
staff.  The  purpose  of  the  visit  was  to  determine  future  target 
characteristics  so  that  performance  requirements  night  be  determined 
more  realistically.  While  the  major  portion  of  the  discussions  was 
classified,  an  abbreviated,  unclassified  set  of  performance  require¬ 
ments  has  been  developed  for  use  in  this  contract  (see  Section  17). 


SECTION  XU.  SUMMARY 


This  report  contains  the  documentation  performed  during 
the  period  20  October  1978  through  7  August  1979,  and  referred  to  herein 
as  Phase  H.  It  is  intended  to  meet  the  documentation  requirements 
specified  in  Farm  11*23  and  in  the  Documentation  Addendum  to  Fora  11*23 
found  on  pages  17  and  18  of  the  contract  (Reference  5)  and  the  require¬ 
ments  that  comprise  Task  II  of  the  Statement  of  Mark  of  Technical  Require¬ 
ment  No.  T-0107  (Reference  l).  The  latter  have  been  repeated  (far  the 
reader's  convenience)  in  Section  I,  Part  A. 

Task  I  documentation  requirements  were  met  by  the 
Technical  Report:  Phase  1.3  Some  of  the  Statement  of  Work  requirements 
have  been  augmented  subsequently  by  this  report. 

Requirement  2.1  was  met  by  Section  m,  Part  C 
of  Technical  Report:  Phase  I. 

Requirement  2.2  was  met  by  Section  17  of  Technical 
Report:  Phase  I  and  augmented  by  Section  17  of  this  report. 

Requirement  2.3  was  met  by  Section  7  of  Technical 

Report:  Phase  I. 

Task  II  documentation  requirements  and  Statement  of  Work 
requirements  are  to  be  met  by  this  report. 

Requirement  2.U  is  met  by  Section  7,  Parts  A  and 
B  and  by  Appendices  A,  B,  and  C. 

Requirement  2.5  is  met  by  Section  7,  Part  C. 

Requirement  2.5  was  deleted.2 

Requirement  2.7  is  met  by  Appendices  A  through  E. 

Requirement  2.8  is  met  by  Section  71. 

With  regard  to  meeting  the  Documentation  requirements  of 

Task  I: 

Paragraph  2. a.  A  stability  analysis,  summary,  and  imple¬ 
mentation  of  the  discrete  control  ...  via  parameter  space,  including  com¬ 
parisons  with  extant  methods,  is  documented  in  Section  7,  Part  B. 

Paragraph  2.b.  All  simulation  initializing  operations 
are  described  in  Appendices  D  and  E. 


Paragraph  2.c.  All  simulation  input  parameters 
are  described  in  Appendices  3  and  E. 

Paragraph  2.d.  Analytical  results  produced  are 
documented  in  detail  in  Appendices  A,  B,  and  C. 

Paragraph  2,e.  Operating  instructions  are  implicit 
in  the  formats  documented  in  Appendices  D  and  E. 

Paragraph  2.f.  Dele ted. ^ 

Paragraph  2.g.  Surmaries  of  computer  simulation 
runs,  results,  and  analyses  are  documented  in  Section  V,  Part  C,  and 
in  Appendices  D  and  E, 

In  summary,  this  report  describes  the  performance  results 
of  an  effort  of  forty-one  weeks  duration  (Phase  H).  The  specific  State¬ 
ment  of  Hark  requirements  were  met  in  detail,  i-*e. ,  the  parameter  space 
technique  has  been  fully  developed  for  stability  determination  of  a  digital 
control  system,  and  the  dynamics  of  a  missile  system  investigated  using 
that  technique.  Also,  an  in-depth  technical  seminar  describing  the  use 
of  digital  control  techniques  was  presented.  In  addition,  the  Principal 
Investigator  developed  and  documented  an  analytical  design  tool  (termed 
"SAM")  to  enable  the  control  system  engineer  to  derive  the  closed -loop 
transfer  function  (in  particular)  and  any  selected  state  (in  general)  of 
a  digitally-controlled  system.  Also,  he  developed  and  documented  a 
means  (termed  "cross-multiplication";  of  finding  the  response  of  a  digitally- 
controlled  system.  He  co-authored  two  papers  describing  the  work  of 
which  this  contract  is  a  part  being  performed  by  the  U.  S.  Army  MHIADC01I 
G&C  Directorate  to  develop  an  Advanced  G&C  System,  both  of  which  were 
accepted  for  presentations  at  AIAA  national  conferences.  He  made  several 
presentations  to  personnel  of  the  G&C,  Advanced  Sensors,  Aeroballistics, 
and  Research  Directorates  describing  this  work  in  order  to  gain  their 
cooperation.  To  gain  important  inputs  for  this  work.  Dr,  Seltzer  made 
trips  to  the  Pentagon  and  to  Eglin  Air  Farce  Base.  Also,  he  performed 
an  evaluation  of  the  theory  and  design  of  the  digitally-controlled  SIG-D 
as  a  portion  of  the  formal  review  of  that  system.  Finally,  he  wrote  a 
"white  paper"  describing  how  the  G&C  Analysis  Group's  Hardware -in-the - 
Loop  simulation  is  used  by  that  Group  to  meet  their  HELLFIRE  mission. 

All  of  the  foregoing  was  accomplished  with  the  expenditure  of  approxi¬ 
mately  lliOO  direct  scientific  man-hours  (the  equivalent  of  17?  eight- 
hour  days)  or  an  average  of  approximately  3U  hours  per  week  for  Phase  II). 
Included  within  this  time  was  the  preparation  and  conduct  of  I23  contact 
hours  of  an  Intensive  technical  seminar. 


SECTION  nn.  COHCHJSIOMS 


The  requirements  set  forth  in  the  contract  have  been  net. 

Three  major  analytical  tools  have  been  developed  to  aid 
the  G&C  system  engineer  to  design  a  digital  system.  They  are: 

1.  A  means,  termed  "SAW ,  of  determining  the  states 
of  a  digital  system  in  terms  of  that  system's  transfer  functions  and  the 
inputs  to  the  system. 

2.  A  means,  termed  " Cross-Multiplication" ,  far 
obtaining  the  response  of  a  digital  system  from  its  closed -loop  transfer 
function. 

3.  A  means,  termed  "Parameter  Space",  for  determining 
the  stability  and  dynamic  characteristics  of  a  digital  system  in  terms 

of  several  selected  system  parameters. 

The  three  design  tools  have  been  compared  to  extant  tech¬ 
niques  and  applied  to  the  design  of  a  typical  missile  system.  They  appear 
at  first  blush  to  be  superior  to  other  classical  methods.  3h  particular, 
they  have  been  developed  with  the  aim  of  being  easily  implemented  an  digital 
calculators  or  computers. 
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1.  INTRODUCTION 


This  report  describes  an  alternative  to  the  use  of 
Signal  Plow  Graphs  (SFG)  and  Mason's  Gain  Rule  (Formula)  for 
analysis  of  complicated  sampled-data  control  systems. 

Usually  in  such  systems,  block  diagram  algebraic 
manipulation  may  become  unwieldly,  particularly  when 
such  systems  include  multiple  loops  and  samplers.  The 
Systematic  Analysis  Method  (SAM)  may  be  applied  to 
such  systems,  as  well  as  to  simple  single-loop  feedback 
systems.  This  is  shown  in  Section  2.  Also  shown  is 
how  to  apply  SAM  to  make  use  of  modified  z-transforms 
(Section  5) . 

The  advantages  of  using  SAM  are  that  the  cumbersome 
application  of  Mason's  Gain  Formula  can  be  avoided. 

Further,  the  entire  method  of  drawing  Signal  Flow 
Graphs  may  be  circumvented.  Since  only  the  equations 
describing  the  system  are  needed  for  SAM,  even  the 
customary  block  diagram  is  not  needed. 

If  the  analyst  prefers  to  use  one  of  the  Signal 
Flow  Graph  methods,  a  modified  SFG  technique  is  also 
described  (Section  3) .  It  is  simpler  and  less 
cumbersome  to  apply  than  the  conventional  SFG  method, 
which  for  purposes  of  comparison  is  described  in 
Section  4. 

All  three  techniques  are  applied  to  two  examples. 

This  is  done  to  better  describe  the  application  of  SAM 
and  the  modified  SFG  and  to  help  provide  a  basis  for 
comparison  (Section  6)  of  the  three  methods. 

To  obviate  searching  for  such  a  description.  Mason's 
Gain  Rule  is  described  in  the  appendix. 

2.  SYSTEMATIC  ANALYTICAL  METHOD  (SAM) 

SAM  is  implemented  by  performing  the  following  four 
steps.  If  the  equations  resulting  from  the  first  three 
steps  are  placed  in  a  table  of  three  columns  (one  for 
each  step) ,  they  are  easily  manipulated  to  perform  the 
fourth  and  final  step. 

STEP  NO.  1.  OBTAIN  "SYSTEM  EQUATIONS" 

The  equations  describing  the  system  are  written  in 
the  Laplace  domain.  If  the  system  is  described  by  block 
diagram,  the  "system  equations"  are  written  upon 
inspection. 
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STEP  NO.  2.  OBTAIN  "MODIFIED  SYSTEM  EQUATIONS" 


If  any  of  the  "system  equations"  contain  terms  that 
in  themselves  contain  the  product (s)  of  an  unsampled 


system  variable  and  an  unsampled  transfer 


e  unsampled  variable  must  be  replaced  by  an 
expression  containing  no  unsampled  variable (s) .  In 
complex  systems,  this  may  require  a  chain  of  several 
substitutions . 


An  "unsampled  variable"  is  recognized  as  one  upon 
which  the  pulse  transform  operation  has  not  taken  place. 

For  example,  when  the  pulse  transform  of  a  Laplace 
transform  function  and/or  variable  is  taken,  that  operation 
is  denoted  symbolically  by  placing  an  asterisk  immediately 
following  the  expression,  yielding  a  so-called  "starred 
quantity."  For  example,  the  pulse  transform  of  the  Laplace 
function  F(s)  is  denoted  as  F*(s).  One  manner  of  expressing 
F*(s)  in  terms  of  F(s)  is 


F*  (s) 


00 


(s  +  i  2 n  n/T) 


n»  —oo 


where  T  represents  the  sampling  period. 

STEP  NO.  3.  OBTAIN  "PULSED  SYSTEM  EQUATIONS" 

Pulse  transforms  are  now  taken  off  each  side  of  the 
"modified  system  equations,"  yielding  "pulsed  system 
equations."  Now  all  system  variables,  either  unsampled 
or  sampled  (pulsed) ,  may  be  solved  for,  either  in  the 
"system  equations"  or  the  "pulsed  system  equations." 

i.  _  — 

STEP  NO.  4.  OBTAIN  DESIRED  INPUT/OUTPUT  RELATIONSHIPS 


All  system  variables,  both  starred  and  unstarred, 
may  be  found  in  either  the  "system  equations"  (Step  No.  1) 
or  the  "pulsed  system  equations"  (Step  No.  3) .  The 
desired  output (s)  may  be  solved  for  by  selecting  the 
appropriate  "system"  or  "pulsed  system"  equations, 
substituting  as  necessary.  This  will  be  brought  forth  in 
the  examples. 
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EXAMPLE  NO.  1 


Given:  The  digital  system  of  Example  No.  1  is 
described  by  the  block  diagram  of  Figure  1. 

To  Find:  The  continuous -data  and  pulsed  (sampled) 
outputs,  C(s)  and  C*(s),  respectively,  in  terms  of  the 
system  input  R(s)  and  the  system  transfer  functions 
G(s)  and  H(s) . 

STEP  1.  SYSTEM  EQUATIONS1 

Obtain  these  equations  directly  upon  inspection  of 
Figure  1. 


E  -  R  -  HC 
C  *  G  E* 


(2) 

(3) 


R  ( s) - ~Q 


■ 

■ 

■ 

■ 

C*  (s) 


Figure  1.  Block  Diagram  for  Example  No.  1. 


1.  Note:  In  the  sequel,  the  following  shorthand  notation 
will  be  used: 

F  =*  F(s)  =^£{f(t)]  ,  i.e.  the  Laplace 
transform  of  F(t), 

F*  F* (S) , 

GH*  [G (s) H (s) ]  ‘ 
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STEP  4.  OBTAIN,  C*,  C 

First  solve  Equation  (8)  for  E 


Substituting  this  expression  for  E*  into  Equation  (9) , 
one  obtains 

C*  -  [ - G-_-  -1r«.  (11 

I  (1  +  HG* ) J 

"To  obtain  C,  one  substitutes. . .yielding" 


EXAMPLE  NO.  2. 

Given:  The  digital  system  described  by  the  block 
diagram  of  Figure  2. 

To  Find:  The  continuous-data  and  pulsed  (sampled) 
outputs,  C(s)  and  C*(s),  respectively,  in  terms  of  the 
system  input  R(s)  and  the  system  transfer  functions. 


D*  (s) 


Figure  2.  Block  Diagram  for  Example  No.  2. 


In  the  interest  of  being  systematic,  the  analyst 
may  wish  to  assign  "states"  (X j )  to  the  various  system 
variables.  The  "system  equations"  (Step  No.  1)  are 
written  from  inspection  of  Figure  2.  In  order  to  im¬ 
plement  Step  No.  2,  the  "system  equations"  of  Step  No.  1 
must  be  checked  for  possible  products  of  unsampled  sys¬ 
tem  variables  and  unsampled  transfer  functions.  Two  such 
products  exist:  G1X2  in  Equation  (15)  and  G2X3  in  Equa¬ 
tion  (16) .  The  first  product  is  readily  manipulated  into 


the  approved  form  by  making  use  of  the  relation  between 
X2  and  X*i  of  Equation  (14) ,  which  is  substituted  into 
Equation  (15)  to  yield  Equation  (17) .  The  product,  G2X3, 
is  easily  handled  by  substituting  Equation  (17)  into 
Equation  (16) ,  then  substituting  Equation  (13)  for  Xi,  and 
solving  the  resulting  expression  for  X4,  yielding  Equation 
(18) .  This  completes  Step  No.  2  (the  second  column  of  the 
array  of  Table  2) .  The  third  column  (Step  No.  3)  is 
obtained  merely  by  "starring"  each  side  of  each  of  the 
equations  in  the  second  column. 


TABLE  2.  SYSTEM  EQUATIONS  FOR  EXAMPLE  NO.  2 
Sys.  Eqs .  Mod.  Sys.  Eqs .  Pulsed  Sys.  Eqs. 


E 

E*5  X 


Xj  -  R  -  X4 


Y  5  X, 


D»GjX2  +  Xx 

g2x3 


(13) 

(14) 

(15) 

(16) 


R  -  X4 

(13) 

X  • 

*1 

(14) 

D*G1X1*  +  X1 

(17) 

G2(D*G,X1*  +  R) 

(18) 

1  +  G, 

E*  »  X*x 
E*  -  X2* 
Y*  -  X3» 
C*  -  x4* 


R*  -  X4*  (19) 
Xx*  (20) 
D'Gj^*  +  X3*  (21) 
G3*D*Xx*  +  R3*  (22) 


where 


d/GlG2 


i  +  G, 


(23) 


and 


A  G5R 

*,  —2— 

1  1  +  G. 


(24) 


To  find  C*,  one  merely  substitutes  Equation  (19) 
into  Equation  (22) : 


C*  =  X  *  *  G  *D*X  *  +  R  * 
4  3  11 


=  G3*D*(R*  -  X4*)  +  Rx* 
G3*D*R*  +  Rx* 

1  g3*d* 


(25) 


In  a  similar  manner,  one  finds  C  by  substituting 
Equation  (15)  into  Equation  (16) ,  substituting  Equations 


(14)  and  (13)  for  X2  and  X^,  respectively,  and  finally 
substituting  Equation  (19)  for  the  resulting  Xi*,  i.e. 

c  -  x4  =  g2x3 

*  G2(D*G1X2  +  Xx) 


g2(d*g1X1*  +  R  -  X4) 


g1g2D*X1*  +  G2R 
1  +  G. 


»  G^Xj*  +  Rx 

»  G3D*(R*  -  X4*)  +  Rv  (2 

Since  X4*  =  C*  was  just  determined  in  Equation  (25) ,  that 
value  is  substituted  into  Equation  (26)  to  yield: 


(*D*R*  -  G3*D*R*  -  Rx* 


g3*d*. 


[/g,*d*r*  +  R,A 

r*  -  ^  i  +  g3*d*  -)  +  R1 

(R*  +  G  *D*R*  -  G,*D*R*  -  Rx*\ 

- - - - -  +  R, 

1  +  g3*d*.  )  1 

/G  D*  ^ 

=  \l  +  G-*D*)  <R*  -Rl*>  +  Ri  •  C 

An  alternate  form  for  (R*  -  R^)  may  be  found  if  desired: 

**  -  V  *  R*  -(rraj 

im-w 

■W  • 

/  R  \  *  d.  _  * 

(i  +  g2  )  r2  •  (i 
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This  expression  for  (R*  -  R  *)  may  be  substituted  into 
Equation  (27)  to  yield 

G.D* 

C  “  1  +  G3*D*  R2*  +  Rl' 

which  may  be  slightly  simpler  in  form  than  Equation  (27) . 
3.  MODIFIED  SIGNAL  FLOW  GRAPH  TECHNIQUE 

If  the  analyst  prefers  to  use  a  Signal  Flow  Graph 
(SFG)  technique,  the  following  modified  SFG  is  proposed. 
It  incorporates  many  of  the  features  developed  in  SAM. 

As  such,  it  appears  to  be  simpler  to  implement, 
requiring  the  application  of  Mason's  Gain  Rule  at  only 
one  stage  of  the  analysis. 

The  first  three  steps  are  identical  to  those  of  SAM. 

STEP  NO.  4.  CONSTRUCT  EQUIVALENT  SFG 

The  equivalent  SFG  is  drawn  directly  from  the 
information  contained  in  the  "system  equations"  (Step  No. 
1)  or  from  the  block  diagram. 

STEP  NO.  5.  CONSTRUCT  SAMPLED  SFG 

The  sampled  SFG  is  drawn  from  the  "pulsed  system 
equations"  (Step  No.  3)  . 

STEP  NO.  6.  CONSTRUCT  COMPOSITE  SFG 

This  is  achieved  by  connecting  the  output  nodes  of 
the  samplers  in  the  equivalent  SFG  to  the  nodes 
representing  those  same  quantities  on  the  sampled  SFG. 

STEP  NO.  7.  OBTAIN  DESIRED  INPUT/OUPUT  RELATIONSHIPS 

Mason's  Gain  Rule  (Appendix  A)  is  applied  to  the 
composite  SFG  to  obtain  desired  outputs  in  terms  of 
system  transfer  functions  and  inputs  to  the  system. 

Examples  of  the  use  of  the  modified  SFG  method 
follow.  For  comparative  purposes,  the  same  two  examples 
to  which  SAM  was  applied  will  be  used. 
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EXAMPLE  NO.  1. 


STEP  NOS.  1-3.  SYSTEM,  MODIFIED  SYSTEM,  AND  PULSED 
SYSTEM  EQUATIONS 

Repeat  these  steps  as  shown  in  Example  No.  1  of 
Section  2  (SAM) .  They  are  summarized  in  Table  1  and 
are  Equations  (2),  (3);  (4),  (3);  and  (8),  (9), 
respectively. 

STEP  NO.  4.  CONSTRUCT  EQUIVALENT  SFG 

The  equivalent  SFG  is  constructed  directly  from  the 
"system  equations"  of  Step  No.  1,  Equations  (2)  and 
13) .  The  resulting  SFG  is  shown  as  Figure  3. 


-H 


Figure  3.  Equivalent  SFG  for  Example  No.  1. 

STEP  NO.  5.  CONSTRUCT  SAMPLED  SFG 

The  sampled  SFG  is  constructed  directly  from  the 
"pulsed  system  equations"  of  Step  No.  3,  Equations  (8) 
and  (9) .  The  resulting  SFG  is  shown  as  Figure  4. 


-HG* 


Note:  Encircled  number  refers  to  identification  of  loop. 
Figure  4.  Sampled  SFG  for  Example  No.  1. 

STEP  NO.  6.  CONSTRUCT  COMPOSITE  SFG 

The  composite  SFG  is  constructed  by  joining  the  two 
SFG's  of  Figures  4  and  5^  in  the  prescribed  manner.  In 
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this  example  a  single  line  connecting  the  E*'s  of  the 
two  SFG's  is  required  (Figure  5) . 


-H 


Note:  Encircled  number  refers  to . identification  of  loop. 
Figure  5.  Composite  SFG  for  Example  No.  1. 


STEP  NO.  7.  OBTAIN  C*,  C 

Looking  at  Figure  5,  one  sees  that  there  are  two 
inputs  to  the  system,  R  and  R*.  Applying  Mason's  Gain 
Rule  (appendix) ,  one  finds  only  one  possible  forward 
path  from  R*  to  C*  and  none  from  R;  hence,  k  =  1.  The 
gain  along  that  forward  path,  M^,  is  seen  to  be 

Mk  *  Mx  =  G*  .  (30) 

There  is  a  single  loop  whose  gain  is 

Kx  =  -  HG*.  (31) 

The  value  of  A  is  found  to  be 


A  =  1  -  Kx  =  1  +  HG*.  (32) 

Since  the  single  forward  path  touches  the  single  loop  of 
this  system, 

4ks4l=L 

The  gain,  M,  between  C*  and  the  input  R*  is  then 


(33) 


Solving  Equation  (34)  for  C*  we  obtain  the  same  result  as 
Equation  (11) . 

In  solving  for  C,  we  find  that  only  input  R*  has  a 
forward  path  to  C.  In  this  case  the  gain  along  that 
path  is  seen  to  be  G  (Figure  5) ,  i.e. 

Mx  =  G.  (35) 

The  value  of  A  remains  the  same  as  that  shown  in  Equation 
(32)  ,  as  does  the  value  of  K.  remain  as  shown  in 
Equation  (31) .  The  single  forward  path  touches  the  single 
loop  of  the  system,  so  Equation  (33)  still  applies. 

The  gain,  M,  between  C  and  the  input  R*  is  then 


C_  _  M1A1 
R* 


G  __ 
1  +  HG 


(36) 


Solving  Equation  (36)  for  C  one  obtains  the  same  result 
as  Equation  (12) . 


EXAMPLE  NO.  2 


STEP  NOS.  1_-  3 

Repeat  these  steps  as  shown  in  Example  No.  2  of 
Section  2  (SAM) .  They  are  summarized  in  Table  2. 

STEP  NO.  4.  EQUIVALENT  SFG 

Construct  the  equation  SFG  directly  from 
"system  equations" (13)  through  (16)  (Figure  6) . 


-1 


Figure  6.  Equivalent  SFG  for  Example  No.  2. 


STEP  NO.  7.  OBTAIN  C*,  C 

From  Figure  8  one  sees  that  there  are  three  inputs 
to  the  systiml  r7  R*,  and  R.*.  Applying  Mason's  Gain 
Rule,  one  finds  two  forward  paths  to  C*, one  from  R*  and 
one  from  R^*,*  hence,  k  =  2. 

There  are  two  loop  gains,  denoted  herein  as  and  K2« 
The  loops  are  designated  by  encircled  numbers  on 
Figure  8,  and  their  gains  are 

Kx  =  -G2,  (37! 

K2  «  -G3*  D*.  (38] 

It  is  observed  that  the  loops  are  nontouching  (necessary 
information  for  formulating  A)  .  A  is  thus  found  to  be 

A  *  1  -  (  Kj_  +  K2)  +  KlK2 
*  1  -  (-g2  -g3*  D*)  +  g2  g3*  D* 

-  (1  +  G2)  (  1  +  G3*  D*)  .  (39 

The  forward  path  from  R*  to  C*  may  be  designated  as 
k  =  1.  Since  it  is  touched  by  Loop  2  but  not  by  Loop  1, 
the  value  of  A1  is 

*1  "  1  -  K1 

=  1  +  G2.  (40 

The  forward  path  from  R.*  to  C*  may  be  designated  as 
k  *  ?.  Since  it  is  touched  by  Loop  2  but  not  by  Loop  1, 


The  gain  along  the  first  forward  path  (k  =1)  is, 
from  Figure  8, 

m  G3*  D*. 

The  gain  along  the  second  forward  path  (k  =  2)  is 


■ 


I 


The  gain  between  R*  and  C*  is 

1  d.  C1*  Ml*l 

M  =  R5-  *  ~A“ 

(,g3*d*)  (1  +  g2) 

(1  +  G2) (1  +  G3*D*) 

G3*D* 

1  +  G3*D*  * 

2 

The  gain  M  between  R^*  and  C*  is 

2  d.  C2*  _  ^2^2 

M  *  r]T  =  — - 

(1)  (1  +  g2) 

=  (1  +  g2)  (1  +  g3*d* 

1 

*  1  +  G3*D*  * 

1  7 

Solving  Equations  (44)  and  (45)  each  for  C  *  and  C  *, 
respectively, 

where 

C*  =  C1*  +  C2*, 

one  finally  obtains 

C*  =  M1  R*  +  M2  Rx* 


G,*D*R*  +  R  * 

-  1-T-g3^d*  1  ' 


(44) 


(45) 


(46) 


(47) 


which  is  seen  to  be  identical  with  the  earlier  SAM  result 
of  Equation  (25) . 


To  obtain  C,  one  must  first  observe  from  Figure  8 
that  all  three  inputs  to  the  system  can  find  their  way 
to  the  node  representing  C.  The  two  loop  gains,  K.  and 
K_,  are  the  same  for  finding  C* ,  as  is  A.  The  forward 
pith  from  input  R*  to  C,  designated  as  k  =  1,  has  a 
gain  of 

»  D*  G2.  (48) 

The  gain  along  the  path  between  R^*  and  C,  designated  as 
k  *  2,  has  a  gain  M2  of 

M2  *  -D*  G1  G2  .  (49) 

Finally,  the  gain  along  the  path  between  R  and  C, 
designated  as  k  =  3,  has  a  gain  of 

M3  -  G2  .  (50) 

The  k  =  1  and  k  *  2  paths  touch  both  loops,  so  the  values 
of  Ai  and  A2  are  both  unity.  The  k  =  3  path  only  touches 
Loop  1,  so  the  values  of  A3  is 

A3  =  1  -  K2  ■  1  +  G3*  D*.  (51) 

The  gain  between  R*  and  C  is 

1  d.  C1  M1A1 
M  =  R7  =  ~L~ 


(52) 


(53) 
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The  gain  M  between  R  and  C  is 


M3  = 


C3  _  M3A3 


G2(l  +  G3*D*) 

a  +  g2>  a  +  g3*d*) 


r+"®2' 

12  3 

Solving  Equations  (52)  through  (54)  for  C  ,  C  ,  and  C  , 
respectively, 

where 


(54) 


C1  +  C2  +  C3  , 


(55) 


one  obtains 
C  i 


M1  R*  +  M2  R^  +  M3  R 


D^Gj^Gj  R* 


d*gxg2  rx* 


g2  R 

(1  +  g2) (1  +  g3*d*) ■  (1  +  g2) (1  +  g3*d*)  +  1  +  G 


D*G? 

ih’g,»d«)  (R*  -  V>  +  R1  • 


2 

'56) 


Equation  (28)  may  be  used  in  an  attempt  to  simplify  the 
above  result,  yielding 


D*G, 

J  R  *  +  R 

(1+G,*  D*)  K2  K1  * 


(57) 


This  is  equivalent  to  Equation  (29)  obtained  using  SAM. 
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4.  THE  SAMPLED  SIGNAL  FLOW  GRAPH  METHOD 


The  standard  Signal  Flow  Graph  method  in  use  is  the 
"Sampled  Signal  Flow  Graph  Method."2  So  that  the  SAM 
and  modified  SFG  methods  exposed  in  Section  2  and  3, 
respectively,  of  this  report  may  be  compared  to  this 
standard  method,  it  is  described  briefly  in  this 
section.  The  same  two  examples  that  have  been  used 
previously  in  this  report  are  used  in  this  section  to 
better  permit  comparison  of  the  various  methods.  The 
other  popularly  used  SFG  method,  "The  Direct  Signal 
Flow  Graph  Method,"  will  not  be  described  herein. 3 

STEP  NO.  1.  CONSTRUCT  EQUIVALENT  SFG 

This  is  equivalent  to  Step  No.  4  of  the  modified  SFG 
procedure  of  Section  3. 

STEP  NO.  2.  CONSTRUCT  SAMPLED  SFG 

Write  system  equations  for  all  noninput  nodes  of  SFG, 
applying  Mason's  Gain  Rule.  A  "noninput  node"  is  defined 
as  a  node  that  is  not  an  input  node,  where  an  "input"  is 
defined  as  a  system  input  or  the  output  of  a  sampler. 

Take  the  pulse  transform  of  each  side  of  each  of 
the  system  equations. 

Using  equations  noted  in  the  paragraph  above,  draw  the 
sampled  SFG  for  the  system. 

STEP  NO.  3.  OBTAIN  RELATION  BETWEEN  SAMPLED  INPUTS/OUTPUTS 

This  is  achieved  by  applying  Mason's  Gain  Rule  to 
the  SFG. 

STEP  NO.  4.  OBTAIN  RELATION  BETWEEN  INPUTS/ 

CONTINUOUS-DATA  OUTPUTS 

Connect  the  SFG's  of  Steps  No.  1  and  No.  2  to 
yield  a  composite  SFG. 

Apply  Mason's  Gain  Rule  to  SFG. 


2.  B.  C.  Kuo,  Digital  Control  Systems ,  SRL  Publishing 
Company,  Champaign,  Illinois,  1977,  pp.  100-105. 


3.  Ibid,  pp.  106-115 
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EXAMPLE  NO.  1. 

See  Section  2  and  Figure  1  for  a  description  of  the 
example . 

STEP  NO.  1.  EQUIVALENT  SFG 


From  the  block  diagram  describing  the  system 
(Figure  1) ,  construct  an  equivalent  SFG.  This  is  the 
same  as  Figure  3. 


STEP  NO.  2.  SAMPLED  SFG 

Write  system  equations,  applying  Mason's  Gain 
Rule  to  equivalent  SFG  (Figure  3) . 

E  *  R  -  G  H  E*  ( 

C  ■  G  E*  ( 

Take  the  pulse  transform  of  each  side  of  the  system 
equations . 

E*  =  R*  -  GH*  E*  ( 

C*  =  G*  E*  ( 

Draw  sampled  SFG  from  pulsed  system  equations 
(Figure  4) . 


STEP  NO.  3. 

Obtain  C*,  E*  from  sampled  SFG  (Figure  4) ,  applying 
Mason's  Gain  Rule.  There  is  one  path  from  R*  to  E*;  k  =  1. 

M1  =  1  (62) 


1  loop:  Kx  =  -GH* 


A»l-K1»l  +  GH* 


The  forward  path  touches  the  loop: 


There  is  one  path  from  R*  to  E  (k  =  2)  ,  and  one 
path  from  R  (k  =  3)  . 


M2  -  -GH 

(77) 

m3  *  1 

(78) 

a2  »  1 

(79) 

A3  =  1  -  Kj  -  1  +  GH* 

(80) 

w2  d.  E2  M2A2  -GH 

»*  ‘  4  "  i  +  at* 

(81) 

rl3  d.  E3  _  **3*3  _  1  +  OI*  _  1 

R  A  1  +  GH* 

(82) 

2  3 

Solving  Equations  (81)  and  (82)  for  E  and  E  , 
and  using  the  relationship, 

respectively, 

E  =  E2  +  E3, 

1  (83) 

one  obtains 

E-R-  GH  R*7"‘ 

1  +  GH* 

(84) 

In  this  elementary  example,  it  is  seen  that  the  SFG's  and 
algebraic  relationships  are  identical  to  those  obtained 
with  the  modified  SFG  method  of  Section  3.  Such  will  not 
be  the  case  with  Example  2. 

EXAMPLE  NO.  2 

See  Section '2  and  Figure  2  for  the  description  of  the 
digital  control  system.  It  is  desired  to  find  C  and  C* 
in  terms  of  system  input  R  and  the  system  transfer 
functions . 

STEP  NO.  1.  EQUIVALENT  SFG 

This  may  be  drawn  directly  from  the  system  equations 
summarized  in  column  1  of  Table  2 .  It  is  shown  as 
Figure  6. 
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STEP  NO.  2.  SAMPLED  SFG 

Write  system  equations,  applying  Mason's  Gain  Rule 
to  SFG  of  Figure  6.  There  is  one  path  from  R  to  E  (k  =  1) , 
one  path  from  R  to  Y  (k  =  2) ,  and  one  path  from  R  to  C 
(k  =  3) .  There  is  one  path  from  E*  to  Y  (k  =  4)  ,  one 
path  to  C  (k  =  5) ,  and  one  path  to  E. 

Using  the  techniques  that  are  by  now  well^  established 
in  this  report. 


M] 

M. 

4 

M- 

M4 

MS 

M< 

K1 

A 

A, 


M  S*  X_ 


M 


*  l' 

(85) 

;  -  lr 

(86) 

1  "  G2, 

(87) 

=  D*G. 

>  t 

(88) 

-  D*Gl  G2, 

- 

(89) 

|  *  -D*G1  G2, 

(90) 

=  ~G2 , 

(91) 

=  1  -  K1  =  1 

+  G2, 

(92) 

=  1, 

(93) 

d.  E1  _  M1A1 

1 

(94) 

R  A 

1  +  G2' 

d.  Y^  _  M2A2 

.  1 

(95) 

R  A 

1  +  G2 

d.  c1  M343 

_  G2 

(96) 

R  A 

i  +  g2  ' 

d.  Y2  _  M444 

D*G. 

m  • 

(97) 

*  E*  A 

1  +  G-, 
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M5  3.  91 

M  -  E* 


M5*5 


D*GiG2 

1  +  g2  ' 


Vt  -°*G2G2 

A  1  +  G2 


One  solves  Equations  (94)  and  (99)  for  E^-  and  E^ , 
respectively,  and  using  the  expression. 


E  *  E1  +  E2, 


(98) 

(99) 


(100) 


one  obtains  E: 


M1R  +  M6E* 


1  ♦  G, 


d*g1g2 

1  +  G- 


E* 


=  R2  -  D*G3E*. 


(101) 


1  2 

One  solves  Equations  (95)  and  (97)  for  Y  and  Y  ,  and 
using  the  expression, 


Y  =  Y1  +  Y2, 


(102) 


one  obtains  Y: 


2  4 

MR  +  ME* 


D*G, 

R2  +  1  +''G  E*' 


(103) 


26 


(104) 


jktkUBHUMWOOi 


Similarly,  one  solves  Equations  (96)  and  (99)  for 
C1  and  C2,  and  using  the  expression, 


C1  +  C2, 


one  obtains  C: 


c  »  m3r  +  m5e* 


G-  R  D*G, G~ 

1  +  g2  1  +  G2 


■  R1  +  D*  G3E*. 


Take 

Equations 


the  pulse  transforms  of  each  side  of  the 
(101) ,  (103)  ,  and  (105) . 


(105) 


(106) 

(107) 


(108) 


Draw  the  sampled  SFG  from  the  pulsed  Equations  (106) 
through  (108)  (Figure  9) .  Note  this  is  not  the  same  as 
the  sampled  SFG  that  resulted  from  using  the  modified  SFG 
method  of  Section  3.  However,  the  final  answers —  the 
desired  responses  —  will  be  the  same,  as  will  be  seen  in 
the  sequel. 
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-D*G3* 


Note:  Encircled  numbers  refer  to  identification  of  loops 

Figure  9.  Sampled  SFG  for  Example  No.  2 
(standard  SFG  method) . 


STEP  NO.  3. 


Obtain  C*  from  the  sampled  SFG  (Figure  9) ,  using 
Mason's  Gain  Rule.  There  are  two  inputs:  R.*  and  • 
They  have  two  forward  paths  to  C*;  k  =  1  and  k  =  2, 
respectively.  Again  using  techniques  that  have  been 
well  established  in  this  report: 

Mx  =  1, 


(109) 


M2  =  D*G3*f 


(110) 


K1  =-D*G3*, 


(111) 


A  *  1  -  K. 


1  +  D*G  * 


(112) 


1  2 

One  solves  Equations  (115)  and  (116)  for  C  *  and  C  * , 
respectively.  Using  the  expression, 

C*  =  C1*  +  C2*, 
one  obtains  C* : 

C*  =  M1R1*  +  M2R2*. 


Equation  (118)  may  be  manipulated  into  the  same  form  as 
Equations  (25)  and  (45) ,  if  desired,  by  substituting 
the  expression  [resulting  from  Equation  (28)]. 


R*  -  R,*  =  R2* 

into  Equation  (118) : 

D*G  * 

C*  =  V  +  T-T D*G3*  R2* 


Rx*(l  +  D*G3*)  +  D*G3*(R*  - 

1  +  d*g3* 


Rx*  +  D*G3*R* 

1  +  d*g3* 


STEP  NO.  4.  OBTAIN  C 

Connect  the  SFG's  of  Figures  6  and  9  to  obtain  a 
composite  SFG  (Figure  10) . 


(117) 


(118) 


(119) 


(120) 
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Using  the  expression, 

C  -  C1  +  C2,  (130) 

1  2 

and  solving  Equations  (128)  and  (129)  for  C  and  C  , 
respectively, one  may  obtain  _  _ _ 

c2R  d*Gig2r2* 

C  =  1  +  G2  +  (1  +  D*G3*) (1  +  G2) 

g.d*r2* 

85  R1  +  ITFg^  *  (131> 

It  is  noted  that  the  values  of  C*  and  C  just  obtained  are 
the  same  as  those  obtained  using  the  SAM  and  modified 
SFG  techniques. 

5.  APPLICATION  TO  2-  AND  MODIFIED  Z-TRANSFORMS 

In  the  SAM  and  SFG  techniques  for  obtaining  sampled- 
data  outputs  of  a  system,  that  form  has  been  indicated 
as  C*  or  C*(s).  The  z-transform  of  C*(s)  is  merely 
written  as  C(z).  Hence,  anywhere  an  expression  C*(s) 
is  found,  it  may  be  replaced  by  C(z)  if  it  is  desired 
to  work  in  the  z-domain  rather  than  the  s-domain. 


If  it  is  desired  to  find  an  output  expression  in 
modified  z-transform,  that  is  denoted  by  the  symbol  C(z,m). 
This  form  may  readily  be  obtained  from  the  expression  for 
a  sampled  output,  such  as  C*  or  C*(s),  by  noting  that  such 
outputs  appear  to  be  equal  to  the  product  of  an  unstarred 
quantity  and  a  starred  quantity.  Let  A(s)  represent  the 
unstarred  quantity,  and  let  B*(s)  represent  the  starred 
quantity.  Then  variable  C*(s)  may  be  written  as 

C*  «  C*  (s)  =  A (s)  B*  (s)  .  (132) 

If  one  recognizes  that  A(s)  or  B*(s)  may  be  equal  to 
unity ,  Equation  (132)  will  always  hold.  _ 

—  The  modified  z-transform  may  always  be  obtained 
from  Equation  (132)  by  performing  the  following 
transformation: 

C(z,m)  =  A(z,m)  B(z),  (133) 

where  A(z,m)  represents  the  modified  z-transform  of  the 
quantity  A(s) ,  and  B(z)  represents  the  ordinary  z-transform 
of  the  quantity  B(s).  This  technique  appears  to  be  an 
attractive  alternative  to  obtaining  modified  z-transforms 
through  SFG  techniques  (which  may  of  course  be  done) . 

6.  COMPARISON  OF  METHODS 

The  Systematic  Analysis  Method  (SAM)  can  be  used  to 
determine  the  states  of  a  digital  system  in  terms  of 
that  system's  transfer  functions  and  the  inputs  to  the 
system.  This  can  also  be  done  by  the  application  of  other 
methods,  such  as  SFG  techniques.  The  usual  advertised 
advantages  of  the  latter,  when  they  are  compared  to  block 
diagram  or  algebraic  manipulation,  is  that  they  are 
particularly  amenable  to  the  analysis  of  complicated 
systems. 

It  has  been  demonstrated  in  this  report  that  SAM 
can  handle  digital  system  analysis  as  capably  as  can 
SFG  methods.  It  has  the  advantage  of  not  requiring 
the  cumbersome  Mason's  Gain  Rule.  Hence,  it  avoids  the 
oft-committed  errors  associated  with  SFG  analysis,  such 
as  overlooked  closed  loops ,  nonobvious  forward  paths , 
etc.  As  with  SFG's,  a  block  diagram  is  not  needed;  the 
system  equations  are  sufficient.  Finally,  it  has  been 
shown  that  SAM  is  easy  to  implement.  It  appears  to  take 
less  lengthy  analytical  manipulation. 


If  the  analyst  prefers  using  SFG's  to  either  block 
diagrams  or  algebraic  manipulation,  a  modified  SFG 
technique  based  on  SAM  techniques  is  proposed.  As  such, 
it  is  systematic.  While  the  SFG's  produced  by  this 
technique  are  usually  different  from  those  produced  by 
standard  SFG  techniques,  they  yield  the  same  results. 
Mason's  Gain  Rule  is  only  applied  at  one  stage  of  the 
analysis  in  the  modified  SFG  technique,  as  opposed  to 
the  standard  SFG  technique  which  requires  several 
applications  of  Mason's  Gain  Rule. 

7.  CONCLUSIONS 

An  alternative  to  the  Signal  Flow  Graph  technique 
has  been  presented  and  compared  to  a  standard  and  a 
modified  SFG.  The  alternative,  termed  Systematic 
Analytical  Method  or  SAM,  is  claimed  herein  to  be  simpler 
and  more  straightforward  to  implement  than  the  SFG  methods 
Not  only  does  it  appear  to  be  quicker  to  use  in  system 
analysis,  but  it  obviates  the  cumbersome  use  of  Mason's 
Gain  Rule. 

For  the  analyst  who  desires  to  use  SFG  techniques, 
a  modified  SFG  technique  is  presented.  It  is  systematic 
and  reduces  the  number  of  required  applications  of 
Mason ' s  Gain  Rule . 


APPENDIX  A  -  REVIEW  OF  MASON'S  GAIN  RULE 


M:  The  gain  (transfer  function)  between  two  nodes 
on  a  Signal  Flow  Graph  (SFG) . 

k:  The  number  of  forward  paths  leading  from  all 
system  inputs  to  a  particular  selected  output. 

M^:  The  gain  along  the  k^  forward  path. 

A^  -  ^  (all  individual  loop  gains) 

+  (gain  products  of  all  possible  combinations 
of  two  nontouching  loops) 

” 2  (gain  products  of  all  possible  combinations 
of  three  nontouching  loops) 

+  ... 

Ak  *  Value  of  A  for  that  part  of  the  graph  not 
touching  the  kt"  forward  path. 
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1.  INTRODUCTION 


This  report  describes  a  technique  for  obtaining  the  response  of  a  digital  control  system.  It  is 
assumed  that  the  closed-loop  transfer  function  is  available  in  thez-  or.  modified  z-transform 
domain.  The  numerator  and  denominator  of  each  side  of  the  transfer  function  are 
cross-multiplied.  The  Real  Translation  Theorem  is  then  applied  to  the  result,  yielding  a 
difference  equation  in  the  time-domain.  ‘This  may  be  solved  for  the  system  response  in  terms 
of  the  reference  (or  other)  input(s)  to  the  system  as  well  as  in  terms  of  system  state  initial 
conditions. 


Two  different  modifications  to  the  basic  technique  are  described  one  using  the  submultiple 
method  and  one  using  the  modified  z-transform  technique.  These  arc  applied  when  it  is  desired 
to  determine  intra-sampling  responses  of  the  system.  All  three  techniques  are  applied  to  a 
single  example.  A  summary  of  the  techniques  and  their  application  is  provided  at  the 
conclusion  of  the  report. 

2.  RESPONSE  AT  SAMPLING  INSTANTS 


It  is  assumed  that  a  given  digital  or  samplcd-data  system  can  be  described  by  a  closed-loop 
transfer  function  that  relates  the  controlled  output  of  the  system  to  the  reference  input.  If  there 
is  more  than  one  input,  the  technique  can  also  be  applied  to  the  resulting  sum  of  closed-loop 
transfer  functions  relating  the  controlled  output  to  each  of  the  inputs.  Although  the  report 
refers  only  to  a  single  controlled  output,  the  technique  can  be  applied  to  find  any  system  state  if 
it  is  related  to  the  inputs  to  the  system  in  the  /.-domain.  These  relationships  may  be  derived  by 
using  any  of  the  standard  techniques  (such  as  signal  flow  graphs)  or  by  the  newly  developed 
SAM  (Systematic  Analytical  Method)  technique'. 


It  is  assumed  that  the  state  whose  response  is  desired  is  denoted  in  the  z-domain  as  C(z).  where 


(1) 


1.  B.C.  Kuo.  Analx sts  and  Synthesis  of  Sampled- Data  Control  Sx stems.  Prentice-Hall,  New  Jer>ey.  I96.V 

2.  S.M.  Sclt/cr.  .V  I  U:  An  Alternative  to  Sampled- Dm  a  .Signal  riou  (iraplts.  I  S  Army  MismIc  K^earch  and  Pe\cl»»pmont 
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The  script  z  denotes  the  operation  of  taking  the  z-transform  further  assume  that  there  is  only 
one  input  into  the  system:  R(z), 
w  here 


R(z)  =•  ^  |r  ( t) I  . 


The  relationship  between  C(z)  and  R(z)  usually  can  be  expressed  as  a  closed-loop  transfer 
function  (or  several  such  transfer  functions)  which  is  a  ratio  of  two  polynomials  in  z,  i.e.. 


C  (z) 

RTzl 


M  j 

^  a  .  zJ 

X  D 

±12 _ 

N  v 

X  bk  z* 

^ _ A  J‘ 


where  coefficients  a,  and  b,  represent  the  system  parameters.  The  procedure  for  finding  the 
response  at  sampling  instants  —  by  the  cross-multiplication  method  —  consists  of  three  steps: 
Step  1.  Cross-multiply  the  numerators  with  the  denominators  of  Equation  (3),  yielding  the 


expression. 


bk  z  c<z>  = 

k  =0 


E  z-^R(z) 

j=0 


bQC(z)  +  b1zC(z)  +  ...  +  bkzkC(v)  +  ...  +  bNzNC(z) 


apRfz)  +  a^zR(z)  +  ...  +  ajZJR(z)  +  ...  +  aMz  R(z) 
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Step  2.  Each  side  of  Equation  (4a)  or  (4b)  is  divided  by  bNz\  The  resulting  equation  is  then 
solved  for  the  C(z)  term  that  is  not  multiplied  by  a  non-zero  power  of  z.  Equation  i.e.. 


c(z)  .  £  +  £ 

°N  N 


+  ...  +  g~ 
N 


aj  z^-Nr(z) 


+  . 


k  *“  ."-"RU)  b°  Z-KC(Z)  Z1-NC(z) 

+  b; 


bk  zk~NC  ( z)  Vl  z-1C(z)  . 

Step  3.  Apply  the  Real  Translation  Theorem  to  Equation  (5),  recalling  that 

OO 

y-1  jc(z)J  =  c*  (t)  =•  ^Pc(nT)  5  (t-nT) 


(5) 


n  =  0 


and 


|c  ( t-kT) |  = 


z  (z) . 


( 6a) 


(6b) 


The  asterisk  is  used  to  indicate  a  variable  that  has  been  sampled,  and  <5(t-nT)  denotes  a  Dirac 
function  occuring  at  the  instant  t=nT.  The  resulting  value  of  c(nT)  at  each  sampling  instant. 
nT,  then  becomes 

a  a. 

c  (nT)  =  r  [  (n-N)  T]  +  r[(n-N+l)T]  +  ... 

dn  n  _ 

+  r-i  r[(n-N+j)T]  +  ...  +  r^«r  [  (n-N+M) T] 


N 


N 


JkM* 

I 

I 

T 

-*» 

i 

4» 


! 

t 

! 


I 

l 

1 

1 

1 

I 
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where  j,  k.  M,  N,  and  n  are  integers. 

The  advantages  of  the  form  of  Equation  (7)  are  three-fold: 

•  The  value  of  c(nT).  for  any  t=nT,  may  be  obtained  for  any  form  of  r(t),  whether  or  not 
it  is  “z-transformable.” 

•  The  expression  for  c(nT)  docs  not  have  to  be  recalculated  every  time  the  form  of  r(t) 
changes,  as  is  the  case  when  the  response  is  determined  by  the  partial  fraction,  power  series,  or 
inversion  formula  methods'. 

•  The  form  of  the  expression  for  c(nT)  permits  the  inclusion  of  initial  conditions,  such 
as  c(0),  if  they  exist. 


FICTITIOUS  SAMPLER 
, - >< — —  C  (z) 


Figure  1.  Closed- loop  sampled-data  system. 


Example  1.  Given  the  closed-loop  sampled-data  system  of  Figure  /,  'the  closed-loop  transfer 
function  easily  is  found  to  be: 


C (z)  G(z)  # 

R  ( z)  1  +  HG ( z) * 


(8) 


If  G(s)  and  H(s)  are  given  to  be 


G  (s)  »  1/  s(s  +  1) 


and 


H (s)  =  1, 


(9) 


6 


respectively,  their  z-transforms  are 


G(z) 


-T 

z  (1  -  e  A) 


(z  -  1) 


) 


(10) 


and 


HG (z)  =  G(z)  .  (ID 

If  one  substitutes  the  expressions  of  Equations  (10)  and  (11)  into  Equation  (8),  one  obtains  the 
closed-loop  transfer  function. 


C (z)  _  (1  ~  e~T)  z 

2  -T 

R ( z)  z  -  2e  z  +  e 


(12) 


One  may  now  apply  the  three  steps  prescribed  for  the  cross-multiplication  procedure, 
obtaining: 


Step  1: 


z2C(z)  -  2e'TzC(z)  +  e"TC(z)  =  (l-e_T)  zR(z)  .  (13) 


Step  2: 


C(z)  =  (1  - 


e~T) z_1R ( z) 


-T  - 1 

2e  z  C(z)  - 


e  Tz  2C(z) 


(14) 


Step  3: 


c(nT)=(l  -  e”T)  r  [(n-l)T]  +  2 e-Tc  [(n-l)T]  -  e_Tc  t(n-2)T]. 

(15) 


If,  as  in  pp.  145-147  of  Kuo’s  book,  r(t)  is  assumed  to  be  a  unit  step  input,  the  system  is  assumed 
to  be  initially  at  rest,  and  the  sampling  period  T  is  assumed  to  be  I  second,  application  of 


7 


Equation  (15)  readily  yields  the  following  values  for  c(nT): 
c(T)  =  0.6321,  the  value  of  c*(t)  at  1  sec. 
c(2T)  =  1.0972  the  value  of  c*(t)  at  2  sec. 
c(3T)  =  1.2067,  the  value  of  c*(t)  at  3  sec. 

c(nT)  =  0.632J  +  0.736  c[(n-l)T]  -  0.368[c  (n-2)T],  the  value  of  c*(t)  at  t  =  n  sec. 

These  values  correspond  to  those  obtained  by  more  tedious  means  in  pages  145-147  of  Kuo’s 
book. 

3.  RESPONSE  BETWEEN  SAMPLING  INSTANTS  USING  THE 
SUBMULTIPLE  METHOD 


If  it  is  desired  to  find  the  intra-sampling  response  of  the  same  type  digital  system  described 
in  Section  2,  it  may  be  accomplished  by  applying  the  submultiple  method  found  in  pages  83-86 
of  Kuo’s  book.  Briefly,  let  c  J^)  represent  the  value  of  the  response  c(t)  at  the  instant,  t  = 
,  where  n-1  represents  the  nurttber  of  intrasampling  responses  desired  (n  is  an  integer 
with  value  greater  than  unity).  If  m  is  an  integer,  the  sampling  period  within  which  the 
submultiples  are  to  be  determined  is  denoted  as  mT.  The  z-transform  of  c  [^J  may  be  found 
from  the  ordinary  z-transform  in  the  following  manner.  In  essence,  c  [^f]  is  the  output  of  a 
fictitious  sampler  which  samples  n  times  as  fast  as  the  real  sampler.  The  z-transform  of  that 
output  is  defined  as 


{=<£>} a' 


C  (z)  =  C  (z) 


n 


z-z 

n 

T-Tn 


where 


,1/n 


n 


and 


Tn  =  T/n. 


(16) 

(17a) 

(17b) 


Now  the  closed-loop  expression  of  Equation  (3)  may  be  altered  to  read 


C<2>n 

R(z) 


M 

t 

N 


(18) 
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9 


I 

T 


«  ft  ' 


4 


i: 

[i 

!.] 


Step  3.  Similar  to  the  procedure  of  Section  2,  the  value  of  eft)  at  the  n,h  submultiple  of  the 
sampling  instant  mT  is: 


+  ...  +  >  r  riElSlWTl-  '!i  c  fiEHllll*... 
bjj  L  n  J  bN  L  n  J 


bk  _  T (m-N+k) T 

bN  L  n 


c  ria±itfl  •  (2D 

d^j  L  n 


where  all  values  of  r[.]  are  equal  to  zero  except  for  those  values  of  r  at  integral  multiples  of  T, 

i.e.. 


r  (AT)  = 


OVZ 


not  integers 


r(t) 


t=£T 


for  l  =  0,  I,  2, . (integers) 


(22) 


and  j.  k,  M,  m,  N,  n  are  integers. 

Example  2.  The  same  sampled-data  system  is  used  in  this  example  as  in  Example  1 .  It  may  be 
shown,  for  the  example  at  hand,  that 


C(z)n  =  G(z)n  E  ( z)  ,  (23) 

where  G(z)„  is  found  in  the  same  manner  as  demonstrated  in  Equation  ( 1 6).  From  Figure  I  it  is 
seen  that 


E(z) 


R(z) 

1  +  HG(zT 


(24) 


Substitution  of  Equation  (24)  into  Equation  (23)  leads  to 

C(2)n  _  G(2)n  . 

R(z)  1  +  HG (z) 


(25) 
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Replacing  z  by  z„  and  T  by  T„  in  Equation  (10)  leads  to  the  expression, 


T 


G(z)  = 


zn(l  -  e~Tn) 


n  (zn  -  1)  (zn  -  e-Tn) 


(26) 


If  one  substitutes  Equations  (17)  into  Equation  ( 1 1)  and  substitutes  the  resulting  equation  and 
Equation  (26)  into  Equation  (25).  one  obtains  the  closed-loop  expression, 


C(z)n  zn(l  ”  e  Tn^  ^  (z  -  e~T) 

R(Z)  (z-1) (z  -e'Tn) ( z  2n-2e"Tz  n+e~T) 

n  n  n  n 


•  (27) 


The  denominator  of  this  expression  is  of  course  (when  equated  to  zero)  the  characteristic 
equation.  It  may  be  expanded  into  the  following  polynomial  in  z„: 


D(z  )  =  z 
n  n 


2n+2 


.  -T  ,  2n+l  ^  -T  2n 
-  (1  +  e  n)  z  +enz 

n  n 


-  2e‘Tz„n+2  +  2e  A  (1  +  e  ±n) 


-T 


-Tm  _  n+1  -  - (T+T  ) „  n 


n 


n 


-  2e  v  n' z 


n 


+  *"V  -  ®'T  <1  +  e_Tnfzh  +  e-(T+T"\. 


(28) 


Similarly,  the  numerator  of  Equation  (27)  may  be  expressed  as  the  polynomial. 


N(z  )  =  (l-e'Tn)  [z  2n+1-  (l  +  e"T)z  n+1  +  e"Tz 
n  L  n  n 


.]• 


(29) 


-  ■ 


II 


Step  I.  If  Equations  (28)  and  (29)  are  substituted  into  Equation  (27)  and  the  resulting 
numerators  and  denominators  cross-multiplied,  one  obtains  the  equivalent  of  Equation  (1 9b): 


2n+l 


2n+l 


Step  2.  If  one  applies  the  procedure  of  Step  2  to  Equation  (30)  one  obtains 


(n+1) 


(n+1) 


(T+Tn) 


(n+2) 


(2n+l) 


(2N+2) 


I 

1 

i 

1 

I 

I 

Ir 

I 

l 

1 

I 

I 


Step  3.  Finally,  one  may  apply  the  procedure  of  Step  3  to  obtain  the  time-domain  difference 

mT  • 

equation  that  yields  the  value  of  c*(t)  at  the  instant  of  time,  t—  —  ,  where  m  is  any  desired 


integer. 


c 


=  (1  -  e~Tn)  jr 

•Tc  j~(m-n)  T  j 


(1  +  e'T)r 


pm-n+1)  T  j 


-T 

+  e  r 


+  ( 1  +  e  ^n) c I 


e  Tnc 


+  2e 


-  2e~T(l+e~Tn)  c  j^iSinzilTj+  2e'(T+Tn)c  ^<nt-n-2)Tj 


(m-2n) T 
n 


+ 


e-T  (1 


-T 

e  n)  c 


pm-2N-l)Tj 


e-(T+Tn) 


(m-2n-2) T 
n 


(32) 


Again,  as  in  Kuo’s  book  and  in  Example  1 ,  it  is  assumed  the  system  is  at  rest  initially,  and  the 
sampling  period  T  is  one  second.  One  may  desire  to  know  two  intra-sampling  values  of  the 
output  c*(t).  In  that  case,  n  =  2  +  1  =  3.  If  one  substitutes  these  numerical  values  into  Equation 
(32),  one  obtains  the  equation. 


c^=o.2835r(^)-0.4866 
+  0.1043r^^)+  1.7165 

-0. 7165c  (53^)+0. 5272c  (^)-l.  26  30c  (5^-) 

+0.5272c(2ji) 

-0. 3679c(^)+0.6315c(^)  -0.2636c(^)  •  (33) 
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It  is  worthwhile  to  pause  a  moment  and  consider  the  power  of  the  difference  equation  just 
obtained.  First,  the  output  c  ||j  can  be  obtained  for  any  deterministic  reference  input  r(t). 
It  only  need  be  specified  at  the  sampling  instants.  Second,  initial  conditions  or  instantaneous 
changes  can  be  accommodated  readily  by  the  equation.  Third,  difference  equations  are 
particularly  amenable  to  programming  on  desktop  calculators  or  digital  computers.  These 
three  advantages  are  not  enjoyed  by  all  methods  found  in  the  standard  textbooks  for 
determining  the  state(s)  of  a  sampled-data  system.  The  submultiple  method  does  suffer  from 
the  drawback  of  having  to  specify  beforehand  the  number  n-1  of  the  intrasampling  instants  for 
which  it  is  desired  to  know  a  given  state,  such  as  c(t)  in  this  case. 

In  the  case  at  hand,  a  brief  demonstration  of  the  calculation  of  c(t)  at  instants  t=0,  T/3, 
2T/3.  T,...,  is  provided  below.  A  unit  step  input  is  assumed  so  that  the  results  may  be 
compared  to  those  of  Kuo  in  pages  195-198  of  his  book.  The  symbol  c  'nj  denotes  the  value 
of  c(t)  at  t=  f|j  sec. 

m=0:  c(0)  =  0 

m=  1 :  c(  1  /  3)  =  0.2835  r(0)  +1.7165  c(0)  =  0.2835 

m=2:  c(2/3)  =  0.2835  r(  1  /  3)  +  1.7165  c(  1  / 3)  =  0.4866 

m=3:  c(  1)  =  0.2835  r(2/3)  +  1.7165  c(2/3)  -  0.7165  c(  1  / 3)  =  0.6321 

rn=4:  c(4 / 3)  =  0.2835  r(  1)  -  0.4866  r(0)  +  1 .7 165  c(  1 )  -  0.7 1 65  c(2/ 3)  +  0.5272  c(  1  / 3)  =  0.8407 

m=5:  c(5/3)  =  0.2835  r(4/3)  -  0.4866  r( I / 3)  +  1.7165  c(4/3)  -  0.7165  c(l)  +  0.5272  c(2/3)- 

1.2630  c(  1  / 3)  =  0.9901 

m=6:  c(2)  =  1.0972 
m=7:  c(7/3)  =1.1464 
m=8:  c(8/3)  =  1.1816 
m=9:  c(3)  =  1.2067 


4.  RESPONSE  BETWEEN  SAMPLING  INSTANTS  USING  MODIFIED 
Z-TRANSFORM  METHOD 

If  it  is  desired  to  find  the  intra-sampling  response  of  the  same  type  of  digital  system 
described  in  Section  2.  it  also  may  be  accomplished  by  applying  the  modified  z-transform 
method.  In  this  method  it  is  necessary  first  to  determine  the  modified  z-transform  equivalent 
of  the  closed-loop  transfer  function  of  Equation  (3).  The  reader's  mind  will  be  refreshed  (as 
needed). 


Let  c(t)  denote  the  response  of  a  digital  system  and  c*(t)  its  sampled  output.  The  value  of  c(t) 
at  the  instant  of  time,  t=(n-A)T.  is  the  value  of  c(t)  that  has  been  delayed  by  an  increment  of 
lime  AT.  The  latter  is  represented  symbolically  as  eft- AT).  If  the  relation. 


m  =  1  -  A,  0  £m  £  1  ,  (34) 

is  used,  the  value  of  c(t)  delayed  by  an  amount  AT  after  the  sampling  instant,  t=nT.  may  be 
denoted  as  c[(n-l+m)T].  Thez-transform  of  c[(n-l+m)T]  is  termed  the  modified  z-transform 
(  y  m)  of  c(t),  i.e.. 


tf  (c(t-T)}  =  2r(c  [  (n-l+m)  t]} 

y  (c(t)}  ='  C(z,m)  =  C(z,A) 
m 

It  also  may  be  shown  that 


A  =  1-m 


C(z,m)  =  z  ^  y{c(t+mT)}  =  z  ^  £  c  [ (m+n) T]  z  n. 

k=0 


(35) 


(36). 


Alternately,  Cfz.m)  may  be  derived  through  its  complex  convolution  definition: 


C(z,m)  =  <^{c  (t-AT)  xdT  ( t)  } 


7-1 


Ts 


-  [^{c(t- AT)}*^{aT(t)  }]  |z=£Ts  (37) 

ioo  - 

I.  C(5)c-"5[l-e-T'=-5>]-V 


7- Z 


where  oi(t)  represents  a  train  of  impulses  each  with  a  unit  area,  the  star  symbol  (★)  denotes  the 
complex  convolution  operation.  ^  denotes  the  operation  of  taking  the  Laplace  transform,  c 
denotes  the  abscissa  of  convergence  (a  positive  real  number),  and  £  represents  a  dummy 
complex  variable.  . 

Equations  (35)  -  (37)  represent  two  methods  of  determining  analytically  the  modified  z- 
transform  (it  is  recommended  that  one  look  it  up  in  a  modified  z-transform  table,  if  one  is 
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available).  The  above  relationships  may  also  be  used  to  determine  the  modified  z-transform  of 
the  transfer  function.  G(z.  m). 

For  the  purpose  of  exposition,  assume  the  modified  z-transform  equivalent  of  the  closed- 
loop  transfer  function  is  known  and  is 


C  (z,  m  ) 
R(z) 


M 

z 

Ifr 

E 


k=o 


(38) 


where  a,  and  b,  may  be  functions  of  m.  The  cross-multiplication  technique  is  similar  to  those 
two  previously  described. 

Step  1.  Cross-multiply  the  numerators  and  denominators  of  Equation  (38),  yielding  the 
expression, 

N  M 

E  bk  zk  C (z,m)  =  E  a_.z^R(z)  (39a) 

k=0  j=0 

b0C(z,m)  +  b^zC^Ztm)  +  b2z2C(z,m)  +  ...  +  bkzkC(z,m) 
+  ...  +  bNzNC(z,m)  =  aQR  ( z)  +  a^Rfz)  +  ...  +  ajZ3R(z) 

+  . . .  +  aMzMR (z) .  (39b) 


Step  2.  Each  side  of  Equation  (39a)  or  (39b)  is  divided  by  b\ZN.  The  resulting  equation  is 
then  solved  for  the  C(z,  in)  term  that  is  not  multiplied  by  a  non-zero  of  z,  i.c.. 


C  (z,m)  =  z-NR(z)  +  ^  z"^15 
bN  bN 


Cl  • 

R  ( z )  +  ...  + 


-""i'Rlz)  +  ...  +  ?!  - 

N 


N 


0  _-N 


z  C  (z,m) 
N 
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+  Zi  + 

N  -1  1 
b  z  J'C(z  ,m) 

+  Vi 

bN 


.  +  r~  z"(N'k)C(z,m)  +  ... 
bN 


(40) 


Step  3.  The  inverse  z-transform  of  each  term  of  Equation  (40)  is  determined,  using  the 
inverse  of  Equation  (39b), 


(41) 


cl  ^ 

C  [(n-1)  T,m]  EC  [(n-l+m)  tJ  =  ^  r[(n-N)T]+  gi  r[(n-N+l)T]  + 

a .  an 

...  +  g^-  r  £(n-N+ j )  tJ  +  . . .  +  g— -  r  £(n—N+m)  tJ 

-|g^  c[(n-l+m-N)T]+  gi  C  [(n+m-N)  tJ+  ... 

+  g^-  C  ^(n-l+m-N+k)T]+  ...  +  g^i  c[(n-2+m)Tj| 


where  j,  K,  m,  M,  N  are  integers. 


Example  3.  Again,  the  same  sampled-data  system  is  used  in  this  example  as  in  Examples  1 
and  2.  The  modified  z-transform  method  will  be  applied.  A  fictitious  time  delay.  e-,|-mlT.  is 
placed  following  the  forward  loop  gain,  Ci(s).  and  a  fictitious  time  advance  (of  the  same 
magnitude  as  the  time  delay)  is  placed  just  before  the  feedback  loop  gain  H(s).  The  closed-loop 
transfer  function  that  results  from  these  two  addition  elements  being  added  is 


C  ( z,m)  _  G  [z  rTii) 

R(z)  "  1  +  HG(z)  * 


(42) 


Using  Equation  (9)  for  G(s).  one  may  obtain  G(z.m): 


G(z,m)  ,  UG(S)}“  U  { s  (S  l  1)} 


ri  ~mT. „  ,  ,  -mT  -T, 

( 1-e  ) z  +  (e  -e  ) 

(z-1) (z-e  T) 
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(43) 


The  value  for  HG(/)  is  given  in  Equation  (9)  -  (10)  which,  with  Equation  (43),  may  be 
substituted  into  Equation  (42)  to  yield 


During  the  n'h  sampling  period,  i.c..  (n-l)T <t<  nT.  where  n  is  an  integer,  one  lets  m  assume 
values  between  zero  and  one  to  find  the  intrasampling  values  of  c(t).  For  instance  if  it  is  desired 
to  verify  the  values  obtained  in  Example  2.  let  m=l  3  and  2/3  and  vary  the  integer  n 
incrementally,  starting  at  n=0.  Again,  for  simplicity,  let  T  =  1  second 


n  =  0 :  c (m-1) . =  0 

n  -  1:  c (m)  =  (l-e“m)  r (0)  =  l-e_m 

m  =  0  (check  case) :  c(0)  =  0 

m  =  1/3:  c (1/3)  =  1  -  0.7165  =  0.2835 
m  =  2/3:  c(2/3)  =  0.4866 
m-1:  c(l)  =  0.6321 


n  =  2:  c (1+m)  =  (l-e_m)  r(l)  +  (e“m-e~1)  r(0) 

+  2e  ^  c  (m) 

m  =  0  (check  case):  c(l)  =  0.6321 

m  =  1/3:  c (4/3)  =  (l-e"1/3)  +  (e~1/3-e-1) 
+  2e_1  c (1/3) 

=  0.8407 

m  =  2/3:  c(5/3)  =  (l-e"2/3)  +  (e"2/3-e-1) 

+  2e_1  c (2/3) 

=  0.9901 

m  =  1:  c (2)  =  (1-e"1)  +  (e"1-®'1)  +  2e_1 


c(l)  =  1.0972 


n  =  3:  c (2+m)  =  (l-e‘m)  r(2)  +  (e-m-e_1)  r(l) 
+  2e  ^  c(l+m)  -e  ^  c (m) 


m  =  0 (check  case):  c(2)  =  (1-e  ) 

+  (e~°-e-1) 


+  2e-1  c(l)  -e"1  c (0) 
=  1.0972 


m  =  1/3:  c (7/3)  =  l-e_1/3  +  e_1/3-e_1 

+  2e-1  c ( 4/3)  -e'1  c(l/3) 
=  1.1464 

m  =  2/3:  c(8/3)  =  l-e_2/3  +  e_2/3  -e_1 

+  2e_1  c (5/3)  -e"1  c(2/3) 


=  1.1816 


m  =  1:  c (3)  =  1-e-1  +  e_1  -e_1  +  2e_1  c(2) 


-  e-1  c(l) 


=  1.2067 


etc. 


5.  CONCLUSIONS 

Several  analytical  techniques  for  obtaining  the  response  of  a  digital  control  system  have 
been  described.  They  are  based  on  a  single  principle:  cross-multiplication  followed  by 
applications  of  the  real  translation  theorems;  Each  is  applied  to  a  single  example.  As  a  starting 
point  for  application  of  each  of  the  techniques,  it  is  required  that  the  dynamics  of  the  digital 
control  system  be  described  in  the  z-  or  modified  z-domain. 
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The  advantages  of  the  three  techniques  over  extant  classical  methods  are: 

•  The  response  may  be  obtained  for  any  deterministic  reference  input  into  the  system  as 
long  as  its  value  is  known  at  the  sampling  instants.  It  need  not  be  described  by  a  differential 
equation,  and  the  /-transform  for  a  specific  reference  input  need  not  be  determined  before 
obtaining  an  expression  for  the  response. 

•  Initial  conditions  and  instantaneous  changes  can  be  accommodated  readily  by  the 
equations  obtained  through  use  of  the  cross-multiplication  methods. 

•  The  difference  equations  obtained  are  particularly  amenable  to  programming  on  a 
desktop  calculator  or  digital  computer. 

•  A  detailed  knowledge  of  the  theory  underlying  digital  or  sampled-data  control  systems 
is  not  required  (although  it  certainly  is  helpful)  by  the  analyst  in  order  to  apply  the  recipes 
described  herein. 


DISTRIBUTION 

No.  Of  No.  of 

Copies  Copies 


Commander 

Defense  Documentation  Center 
ATTN:  DDC-TCA  12 

Cameron  Station 
Alexandria,  Virginia  22314 

Commander 

US  Army  Missile  Research  and 
Development  Command 
Redstone  Arsenal,  Alabama  35809 
ATTN:  DRDMI-X,  Technical  Director  1 

-T,  Director  1 

-E,  Director  1 

-TR,  Director  1 

-TBD  1 

-Tl  (Reference  copy)  1 

-Tl  (Record  Set)  1 

DRSMI-LP,  Mr.  Voigt  1 

Commander 

US  Army  Research  Office 
ATTN:  DRXRQ-PH,  Dr.  R.  Lontz  5 

P.O.  Box  12211 
Research  Triangle  Park, 

North  Carolina  27709 

US  Army  Research  and  Standardization 
Group  (Europe) 

'  ATTN:  DRXSN-E-RX,  1 

Dr.  Alfred  K.  Medoluha 
Box  65 

FPO  New  York  90510 
Commander 

US  Army  Material  Development  and 
Readiness  Command 

ATTN:  Dr.  JameS  Bi  nder  1 

Dr.  Gordon  Bushy  1 

5001  Eisenhower  Avenue 
Alexandria.  Virginia  22333 


Headquarter,  Department  of  the 
Army,  Office  of  the  DSC  for 
Research,  Development  and 
Acquisition 

ATTN:  DAMA-ARZ  1 

Room  3A474.  The  Pentagon 
Washington,  DC  20310 

OUSDR&E 

ATTN:  Mr.  Leonard  R.  Weisberg  1 

Room  3D1079,  The  Pentagon 
Washington,  DC  20301 

Director 

Defense  Advanced  Research 
Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209  1 

OUSDR&E 

ATTN:  Dr.  G.  Gamota  1 

Deputy  Associate  for  Research 
(Research  in  Advanced  Technology) 
Room  3D1057,  The  Pentagon 
Washington,  DC  20301 

US  Army  Materiel  Systems 
Analysis  Activity 

ATTN:  DRXSY-MP  2 

Aberdeen  Proving  Ground, 

Maryland  21005 

IIT  Research  Institute 

ATTN:  GACIAC  1 

10  West  35th  Street 

Chicago.  Illinois  60616 

Control  Dynamics  Company 
ATTN:  S.M.  Seltzer  30 

701  Corlett  Drive.  Suite  2 
Huntsville.  Alabama  35802 


23 


APPENDIX  C.  SAMPLED -DATA  ANALYSIS  IN  PARAMETER  SPACE 
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This  report  descries  a  technique  for  determining  the  stability 
and  dynamic  characteristics  of  a  digital  control  system  in  terns  of 
several  selected  oyster,  parameters •  -he  method  rscu-ir:s  that 
the  system  characteristic  equation  be  available  ir.  the/ z-a c..nn . 
example  is  provided  to  further  elucidate  the  technique. 


' 

) 

i 


SECTION  I.  INTRODUCTION 


I 

I 


T 

I 

I 


^  %he  parameter  spaeS/^HJ  p^y'vi *s  an  analytical  tool  developed 

for  use  in  control  system  analysis  and  synthesis.  Although  not  necessary, 

its  application  is  facilitated  by  augmenting  the  analytical  results 

with  graphical  portrayals  in  a  selected  mult iTparame ter  spo.ee. 

The  method  reauires  that  the  control  system  be  described  by  a  character^ 
istic  equation  which,  for  sampied-aata  or  digital  systems,  if 

essed  ±n"  the  z-donain.  The  technique  is  based  on  the  analysis  and  synthesis  methods 

for  linear  and  nonlinear  control  system  design  which  are 


excellent  monograph  on  the  subject.' 


^ amply  described  in  ^iljak’s 
Reference  2  describes*^) 

Vhcapplicavloii  the  technique  to  the  analysis  and  synthesis  of  linear  sampled  data  control 

_ _ _  systems.  *  "  ' 

" - ^  fcOncc  the  system  change,  istic  eeuatien  has  been  obtained,  the  pamncler  plane  method 

'ii.»Vn,r»j  a\ enables  tiie  designer  to  c.uiu.itc  graphically  o;  ta.e  :  •i,'n.  i le  -ce.  he  may 

design  the  control  system  :n  tarms  of  the  chosen  pericinar.'  C  criteria^  e  g.  absolute  stability^ 


I 


dan 

roots 


aping  ratiq,  and  scuiiitio  time.  He  is  abh- to  sec  the  ertcei  o.ii  the  charactcriitic  equation  /C.|  i  \  J  i.  ■ 

is  of  dtsngim’  two  aiiji'si.ibie  parameters,  bil.iak  further  simplmed  the  design  procedure  f  ey  X  >n,  (j 


£prvr»n\j) 

fcl 


by  introducing  Cltebyshcv  functions  into  the  equations  thereby  putting  them  in  a  lomi 
paiticulariy  amenable  to  their  solution  by  a  digital  computer.  _  ,  . 

The  method  has  been  extended  to  pourtiv  the  ctlect  of  varying  the  sampling’  r-crioa,  D ,  K 


I 

I 

I 

I 


or.f^eir^see  the  effect  of  the  choice  of  values  assigned  to  the 
sampling  period  on  absolute  ar.d  relative  stability.  /--Iso,  the  recursive 
formulas  shown  therein  are  simpler  in  form  than  the  C he by she v  functions 

of  Ref.  2.  Tne  rcsllUn..  formulation  is  dchberatelscasuo^  ^  a?.aiu 

*  '  *  i  V)»*  *\  di  -ital  *'  + 

it  parlieuh  rh  -'nn1'''^1  analysis  and  computing  n..ul..n-  -^,^cn  pcrtrayecj  graphically,  the 

emphasumu  in--  n.»v  i  -  - 


/ 


results  show  the  dynamic  relation  between  the  solocted  parameters  r-rd  the. 

Cr|»T>gl>n»'VH>>i 

characteristic  equation  roots,  as  a  function  o.  v:.o^i::aeper.der.t  argument, 
cU^T.  Honce  one  (oar  dc-0  .\C2/rZr-J  ily  \he  dynamic  effect  upon  the  system  cf  various 
combinations  of  values  of  the  selected  parameters  defining  the  parameter  space. 
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7  The  history  of  the  continuous-tine  domain  version  of  the 

>*> 

parameter  plane  technique  is  well- described  wi+h  ■suitable  references 
-*■  sLn  l37oj 

n  Ref.  1.  Briefly  summarizing  that  history. .  .  ishnerradskv-. of _ 

Vaeve-Lcoed  a r.zj 

.the  Leningrad  "School  of  Theoretical  and  Applied  Mechanic s/^usea  "  1  . 

I 

first  version  of  the  parameter  plane  technique  to  portray  system 
stability  and  transient  characteristics  of  a  third  order  system  on 

a  two-paraneter  plane.  In  1949  Professor  Yu  I.  Neinark  of  the  Russian 
School  of  Automatic  Control  generalized  Vis  hne  grads  kJ^'s  approach  to  permit 
the  decomposition  of  a  two-parameter  domain  (D)  describing  an  nth  order  system 

- 

into  stable  and  unstable  regions 

ChT  technique  was  called  D-deconposition.  During 

the  period  1959-19 66  Professor  D.  Mitrovic,  founder  of  a  Belgrade 

group  of  automatic  control,  extended  the  method  to  enable  the  analyst 

to  relate  the  system's  variable  parameters  to  the  system  response, 

using  the  last  two  coefficients  of  an  nth  order  characteristic  equation. 

Beginning  in  1964,  PrOfessor  D.D.  Siljak,  then  a  student  of  Mitrovic ' S: 

at  the  University  of  Belgrade,  generalized  the  method  and  called  it 

the  parameter  Plane  method.  Kin  j-ethod  mgrmitted  the  analyst  to  select 

Qih'a rv. cz'-:"iszic  e c n a x i c rT^ 

an  arbitrary  pair  o^coefric ier.xs  v or  paia.v.exers  appearing  within  the 
coefficients)  and  portray  both  graphically  and  analytically  the  depend¬ 
ence  of  -the  system  response  upon  the  selected  parameters.  The 

method  was Csubsecuently>Sxtcr.dGa~"bv  ->iljak  and  others  to  encompass 


a  host  of  related  problems. 


“J.  George  modified  the 


ed  the  D-deconposition  method  to 

GTfiS?  Professor  u.  —  a  nul*-oaranete 

enable  the  portrayal  of  tbs  absolute  stabarty  r=  ltrtlllty. 

,  j  v,..,.  +  n  nor4rav  contouis 
space  (George  also  showed  how  -o  por  „ 

as  did  Bil.i^k)*  All  of  the  forcroing  work  is  carefully  ana  completely  referenced 


within  Ref.  1.  In  i960  and  subsequently,  Gcltzer  has  applied  the  parameter 
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v  space  method  to:  the  design, of  missile,  aircraft,  and  satellite  controllers, 

including  systems  containing  one  or  two  nonlinearities;  the  analysis  of  the 
dynamic  effects  of  the  nonlinear  "Solid  Friction"  (Dahl)  model  for  systems 
with  bali  bearings,  such  as  control  moment  gyroscopes  and  reaction  wheels; 
and  the  specification  by  the  system  designer  of  the  dynamic  structural  flexibility 
constraints  to  the  structural  designer.  Most  of  this  work  has  appeared  in  the 
technical  journals  of  the  IEEE ,  the  AIAA ,  the  International  Journal  of  Control, 
and  the  journal  .ne^nnters  1-  Electrical  Fr?ine?r_ing.A  portion 

of  the  history  that  has  not  been  reported  uponjUwith  one  exception;  to 
be  noted)  is  the  control  system  v/ork  conducted  by  the  German  rocket 
scientists  in  the  early  19^0* s  at  Peenemunde.  There,  Dr.  Haeussermann 

and  others  applied  the  D-deconpositior.  technique  to  the  design  of  the 

^  -  _  ^  ""*  > 

V-2. Rocket,  following  ur.  Haeussermann ' s  .earlier^  _(pre- 

^orld  »ar  II)  application  to  the  control  of/|inderwater  torpedo  .  This 

v/ork  was  not  published  in  the  open  literature  because  of  national 

-j  security  constraints.  when  the  group  cane  to  the  United  states  to 

work  v/ith  the  "rrty  ballistic  fissile  Agency  (first  at  ft.  ^liss,  Texas, 

(  then  at  Redstone  arsenal,  Huntsville,  Alabama),  Dr.  Haeussermann  and 

L 

•his  associates  continued  to  apply  the  method  to  UH  ^rny  missiles 

1 

(and  later  to  NASA  space  vehicles),  ''gain,  national  (this  time,  another 
nation!)  security  precluded  publication  in  the  open  literature  until 
1957-  5  ’ 

$  In  19£H  Professor  Sil jak^«i3^the  first  application  of  the 

l 

parameter  plane  technique  to  sampled-data  systems.  '‘•s  mentioned 
above,  this  was  extended  in  Ref.  3”H.In  1971  Geltzer  presented  an 
algorithm  for  systematically  solving  the  Popov  Criterion  applied  to 
>  sampled-data  systems-  .  applications  of  these  sampled-data  parameter 
space  techniques  are  found  in  “of.  3»  H,  and  6, 


SLwTiC.t  Ixa  A  - 


j^he  technique  requires  that  the  control  system  be  described  by  a  characteristic  equation  i- 
which  is  i.i'ar.yffrnu'il  nrte  the  c-dom:un.  Two  adjustable  parameters  ik0,  /.q)  are  selected, 
and  the  characteristic  equation  (CF.)  is  recast  in  tetris  of  them;  i.e. 

CE  =  f  yrJ  =  0,  •  (1) 

j=o 

z=en  =  »«i9,  (3) 

r  =  _  (4) 

B  =  B(C,  T)  =  cos  0  =  cos  K7\  T~-T’ )  s  ^  W  ‘  (5) 

where  «„  And  T  represent  the  damping  ratio,  natural  frequency,  and  sampling  period.. 
respcctively!|To  transform  the  ch.tcacierislic  equation  from  “,ri  nth  order  pe.t-ynorcia 


into  an  algebraic  equation,  zJ  may  be  defined  m  toms  o 


‘presents  all  system  parameters  other  -nan  ,:Q  , 


r-;al  and  imaginary  parts,.  R^and  * 


z3  =  Rj  +  1  ' 

It  readily  follows  from  Bq.  (3)  that  the  values  of  R 


R  =  r-  cos 

J 


I.  =  r3  sin  4*5 

While  Bqs.  (?)  are  satisfactory  for  determining  one  v«lu 
I  a  set  of  recursive  formulas  can  be  derived  that  are 
amenable  to  implementation  on  a  desktop  calculator  or  di 


jtr  $OJsn  <jnb<*&  *  - 


l^yj/J^cj  /?•  }  T- 


v.-lues/fTTha  and  I^*J\  "••” “S  -rc  obtained  from 
«q.  (3)^75  "‘.-n  j  =0,  the  value  of  zJ’  in  uni 


y.  of ini t ion  of  z  ir 
;  and,  fr on  Zq.  (£), 


V.hen  j  = 


Hence , 


z^  =  R„  +  i  I, 


1  *  R0 
0  =  IQ  . 

1,  the  value  of  zJ  is,  in  Zuler  for.n, 
7}  -  z  =  r  cos  0  +ir  sin  8 

=  Rj  +  i  Ij  . 


Rj  =  rB 


( io*0 

(101) 


II  =  ^1-3  • 


(12R) 


(121) 


It  will 
that  th 
Solely 

1i’  30 


turn  out  to  be  useful  in  the  sequel  to  observe  from  ~q.  (3) 
rsdical?*/l-B%  appears  as  a  factor  in  each  value  of  I ^ . 
for  simplicity,  it  nay  be  factored  cut  by  deiiring  a  new  term 


I.'  =  Ii/*V  1-3“  .  ^13) 

3  J  ’ 

" t'tf  C'lg.r j «ri  iiU-h  substituted  into  Ecyiuti»»ii  (1?  and  the  real  and  iniaginan  parts  of  the 
resulting  equation  are  set'ti.a'.cd.  two  siniMiancous  algebraic  equations  are  obtained.  _ 

These  two  equations  contain  t’\p.  adjustable  parameters,  or  variables, 

k^  ,  k^  .  “ence^tne  "two  equations  may  be  solved  £xplicitely  for  kQ 

and  kj  as  functions  of  the  other  system  parameters  (d^)  ar.d,  in  parti 

cular,  as  functions  of  the _ .independent  argument,  WbT. 

It  is  this  latter  observation  that  forms  the  h^sis  of  th?  p..ra: 

C*  y*' ■  ^  ^  ’  tOd  “'  q  V»  .  v»  r  •  *  y'l  *  v*<  -•  .-7  j*.  3*  ^  T  *  «*» « p  “  v  c  v  *  **■' 


LV.J  ..olect- a  system 


jr  ^  U  gij'*vn  ^  funct.o^j  0f  ^  ^5  3  0y  6r  ;?  w, 


fc%.  dA,‘ 


y  1 


SZCTICr:  III.  STABILITY  uZ?Z?Z:ZL,ZliVi 


The  method  involves  the  definition  of  stability  bcund-ries  or. 
a  multi- parameter  (to  include  kQ  ,  k1  )  space.  Those  boundaries  &r*t 
found  from  the  pair  of  simultaneous  algebraic  equations  « :•<*»' 
from  the  follow  inf-  operations  on  the  system  characteristic  equation 
written  in  terns  of  the  complex  variable  ,  z.  equation  (6)  is  suss  —  .- 
uted  into  characteristic  equation  (l),  and  the  resulting  real  ana 
imaginary  parts  ar°  separated  and  equated,  to  ?gyOi  „  T.-|.  “"r‘ ■- 

algebraic  equations  nay  be  solved  for  any  a..'G  p..^  -  ---  -VI ••  ^  .  * 

vv  >  >  ,_Lf+..r?.v  ■  ' .w .1 . '  '?  '-f*  'jy. .V  7  /J  '~l  J'  ^  ’ ~ 

There  nay  be  as  many  as  four  stability  boundaries  for  any  system 
described  by  characteristic  equation  (l)»  althouy.u  it  is  '.o*  nece_.&. j 

for  all  four  to  exist.  .  y 

1,  Cne  stability  boundary  separates  the  stable  complex 

^  *»  '^sbls  onca-  ”  ocnsiE‘s  c; £ , 
map  of  the  unit  eircle^Jp  8-pUn.  onto  a  r.leeted  ^ac.-. 

</////•/■  ill'  _  /Xy,Tv-v\t  0r,»  TAWS 


•*  — - Consider  the  oft-occurrir.g  case  where  the  coefficients  of 

*  '  * 

|  7  the  powers  of  z  in  the  characteristic  equation  are  linear  combinations 
of  Icq  i  f  i.e# 

I  %*  *}  h  +  biK+c,  ,  •  6*») 

# 

I  where  a.  ,  b.  ,  c;  reDresent  all  system  parameters  other  than  kA  ,  k,  . 

J  J  J  0  1 

In  this  case  the  two  simultaneous  equations  resulting  from  the  real  and  imaginary, 
[  i  respectively,  parts  of  the  characteristic  equat  ion  assume  the  form, 

/  Re  [c-£-]  -  A,  o  ♦  8,  h-C,-0. 

I  .  I*  {c.C.J  -  A2  o  ♦  U2  b-C2;0,  US!) 


m 


2.  The  next  tv/o  stability  boundaries  are  those  sop-rating 

the  stable  real  roots  fro:',  the  unstable  ones.  These  comprise  a  mapping 
of  the  0=1  and  s=-l  points  fron  the  z-plane  onto  the  selected  parameter 
space .  Those  are  found  r  lily  by  substitute v  *  s=l  and  z=-l,  respective! 
into  the  characteristic  equation  (l)  -  ch  of  the  tv/o  resulting 

equations  results  in  a  definition  of  the  tv/o  real  root  boundaries 

i,  the  ..l.ot.4  p««.t.p  space.  @ 

3.  The-  fourth  boundary  is  a  nuppingTtr.e  conditions  that 

cause  the  tv/o  simultaneous  al~ebrnic  equations^used  define  the 

connlex  conjugate  root  boundary)  to  become  dependent.  In  the  case 
where  Rq  and  k^  appear  linearly  as  in  Eq.  (l4),  this  caso  is  found 
by  d-ternining  the  conditions  that  cause  the  Jaccbian  of  Eq .  (i3) 

to  become  identically  equal  to  aero. 

for  a  linear  sampled-data  control  system  to  be  stable,  it  is 


nec 

e  scary 

that  all 

roots 

of  the 

characteristic 

equation  lie 

within 

the 

unit  c 

ircle  on 

the  z 

-plane . 

If  the  system 

is  low- pass 

in  nature, 

the 

stable 

j'Ofti  on 

is  bou 

nded  by 

the  semi-circle 

defined  by 

the  upper 

half  of  the  unit  circle^  (part  1,  above)  and  the  singularities  associate 
with  z  =  t  1  (part  2,  above)  and  J=0  (part  3»  above).  The  mapping  of 
those  boundaries  onto  the  parameter  space  will  bound  the  stable  region, 

if  one  exists.  That  ropier,  is  determined  by  applying  a  shading  criterion 

| 

or  using  a  test  point.  If  the  Jacobian  is  greater  than  zero,  then 
the  stable  region  (if  it  exists)  lies  to  the  left  of  the  complex 
conjugate  root  boundary  as  (Ok^T increases 

^the  left  side  of  the  lie:  is  double  cross-hatched  to  indicate  a  boundin'  associated  with 
double,  or  complex  e  jniueate.  roots,  ill  the  Jacobian  is  !es>  than  rare,  to;  stable  rartion  las 
*  to  the  ii;.;ht.)  Single  c/oss-lu'ichins  is  used  on  the  two  contours  associated  with  the  real 
S'  roots.  The  side  of  the  contour  on  which  to  place  the  cross-hatching  is  determined  by  the 

requirement  tlt:’t  cross-liatchiriir  be  continuous.  or  on  the  same  side,  ol  the  contouts  •>'  the 
intersections  c«,rres;  or.dine  to  :  =  ■*  I  and  z  —  —  1  arc  approached  alot.y  either  ,i.c 

complex  root  or  tea1  root  stabtiity  boundary,  q;._  ;u;  1  (V.tt  p'.y.'.ic  l'.-r  ill  ) 


t;'.e  d 


sient  characteristic 


of  the  system  can  he  specified  in  terns  of  the  locations  o-  the  rcces 
of  the  characteristic  equation  (pole  placement).  For  the 
complex  conjugate  ro 'ts,  these  locations  are  defined  in  terns  of  danpi.* 
ratio  $)  --.nd  system  natural  frequency  (^r»0  .  Contours  of  constant  ^ 
are  determined  as  functions  of  and  T  in  precisely  the  sane  r.an-er 
that  the  complex  conjugate  stability  boundary  was  determined  except 
that  $  is  not  set  equal  to  zero. 


1  The  rod  root  locations  corresponding  to  values  of  -  when  0  equals  (f  and  ISCi’^r.iav^Q 
**— olfrfl  be  plotted  on  the  parameter  plane  by  setting  r  equal  to  a  positive  or  negative  real 

constant,  substituting  that  value  into  (I).  and  solving  for  k,  as  a  function  of  k0. 

Each  resulting  contour  corresponds  to  a  location  of  a  real  root  in  the  e-domain.  When 

2  =  +Ct 


/ 

and  when  r  = 


X  VJV  =  0, 

J- o 


X  V-M:y  -  X  V,2;>  *’  =  o,  ii  even 

Jr-  0  (e|J 

(»-  1VJ 

X  li'ifli1  -  "  =  o, II odd 

r  o  _ 


m  ; 


w)  •• 


where  «( is  a  positive  i  cal  mintber.f 


i tut ad  3 


r-~  soondinr 


1T\r  'j  0  d  3  S  1 2T  Q  u  T  9 

)  -  (22)  to  cbtJ-i 
ooatiens cv/  the 


V  The  analytical  technique  developed  pen-uts  the 
taneously  changing  three  control  parameters  an 
convention:. I  technique's  tvtr.iit  the  obicvritmi  oi 
naramelc'  and  don’t  show  the  etirct  »>l  *.m r;> m-  ».i; 


Liner  to  observe  the  effect  of  m:«ui 
.•  N'miditv.:  pet  toil.  Most  exist:.', 
eii’t. cl  of  ehatiCtr.  ’  os«l »  one  contu 


yu»»sj  L  iArn a  .1  ... _ ,m.  -h  L.\ 


sor.'.e tines  the  system  parameters  nay  be  manipulated  within  the  equations 
defining  the  stability  and  dynamic  contours  so  that  more  than  two 
parameters  (such  as  kg  ,  k^  nay  be  used  to  define  a  parameter  space. 

<*n  example  of  a  three-parameter  space  is  provided  tv*.*  ! r ^  ,V|  ^ 


Sometimes  it  is  .-ipcedicd  that  the  seliimg  time  of  the  system  be  less  than  a  prescribed 
value.  This  corresponds  to  requiring  that  the  real  part  of  the  toots  of  the  characteristic 
equation  be  less  than  a  prescribed  negative  real  constant.  A  boundary  corresponding  to  this 
requirement  can  be  drawn  on  the  parameter  plane  by  mapping  a  circle  of  constant  radius 
(for  a  chosen  values  of  <  7  )  from  the  e-plane  onto  the  parameter  plane?'  .  Relations 

exist  for  estimating  the  maximum  overshoot  arid  peak  time  of  transient  rct.non.--e  v-  hen  it  is 
valid  to  assume  a  ..ccottd  order  system?  However,  a  simple  cs  t  i  in  a  t  dLcati/omet  ■  mostly 
made^bv/Tiicrdv'  looking  at  the  dillercnce  equation  representing  the  system  response, 
estimating  when  me  overshoot  will  occur,  and  plotting  a  corresponding  line  on  the  para¬ 
meter  plane.  Tbij  I'lruti-Ll.i.  j  iu  btou  ;rt  out  urtliu  i‘...'!.»|'»^Sic:idy  state  rcsponjejrtay  be 
found  from  the  open  loop  transfer  function  and  tnc  assumed  forcing  functions  in  the 
conventional  manner. 


«► 

i 

•4* 


mm 


I  _  section  X  .  g  X  ^ 

*  Consider  the  planar  model  that  portrays^ 

I  rotational  dynamics,^  -v.r  ~  '  7  It  is  shown  in  block 

j  diagram  form  in  Fig.  l.  |he  pi4nt  to  be  con- 


v3- 


O^r 


c5 l>) 

«.i|  i 


fig.  1  Block  Oiagrari  of  jYlcW 

trolled  is  represented  by  the  transfer  functions. 

G4(s)  e  d(s)/Tc(s)  «  l/Ovs 

G,(s)  s  «(s)/«(s)  *  Vs  -V 

9 

whnrf.  s  represents  the  Laolace  operator  0  re?re‘ 

_  ij  [in  the  tire  Ceram, 

•  TA  c  Lsploce  domain,  ofs)],  Tc. "presents 
the  co  r  an-J.  d  tereve.  and  the  ovcrcot  represent 
the  derivative  with  respect  to  tire.  The  two 

Censors  that  measure  attitude  (v)  and  attitude  rates 
/  \T)  arc  assure-/ to  he  perfect  with  unity  tracer 

’  functions  The  on-board  diqital  controller  de- 

/  -Soo  the  cUandcd  control  torpee  from  the  tnput 
states  o  and  o  and  the  cosmanded  attitude  signal 
dc:  C|,o(s)  represents  a  zero-oroer  hold  in  the 
computer:  i 


Gh  (s)  *  (1-c  Tj)/s. 


fltf  . 


where  T  is  the  sample  period  of  the  on-board  digital 
computer.  The  ?lrJ  control  algorithm  is  reprinted 
by  G^(s)  and  C&(s)  , 

G,(s)  *  K?  4  Kj/S  ^ 

where  Kp  and  K;  arc  the  position  ar.d  integral  feed¬ 
back  gains  ar.d 

W-h  ’  ^  ; 

where  r.n  is  the  derivative  feedback  gain  (r,dre 
cov.icnl!/  called  the  rate  can./. 


The  third  order  characteristic  equation  is  Eq.  (1)  with  n=3»  where  the  coefficients 


Kodificd  gains  a,  b,  c  are  defined  as 


%  The  stability  boundaries  are  presented  in  terns 
Of  a,  b,  c,  and  unT  by  observing  tnat  the  roots 
(in  this  case  three)  cf  a  characteristic  elation 
representing  a  stable  linear  ss."  pled-data  control 
system  rest  all  lie  witnin  the  unit  circle  in  the 
2-plane.  '  It  is  assured  that  the  syster, 

possesses  low-pass  filter  characteristics  so  that 
only  the  primary  strip  (corresponding  to  0  <  0  ±  r) 

'  need  be  considered.  The  stable  region 

may  be  defined  by  capping  its  three  boundaries  from  the  z-plane  onto  the 
a,b,c  parameter  space  by  first  considering  parameters  a,b  to  be  Hq  ,  in  Eq. 

Hs.)  * 

The  boundary  at  z  =  ♦!  is  found  by  substituting 
that  value  for  z  into  the  C.E.,  yielding  the  sta¬ 
bility  boundary 

tii  )?cusrut"s 

W  ■' 

— -  The  complex  conjugate  root  stability  boundary 

may  bo  found  by  sotting  z  equal  to  cos  i*'„r  +  -  rnriv°nient 

.1  sin  ,„T  in  the  C.t.  however,  at  that  poin*  it  IS  convenient 

to  use  the  recursive  formulas  of  Eq .  (S')  and  transform  tne  C.E. 

/  into  two  algebraic  equations  by  separating  the  real  and  inagirary 
parts.  If  *h®y  are  solved  for  a  and  b  and  the  associated  Jacobian  J,  one  obtains  as 
the  complex  conjugate  root  stability  bc'tr.oary, 

a  *  CO  +  B)/0  -  c)j  C  4  o  -  «),  (27)  ^ 

(It)  \ 


b  “  c  ♦  (1  -  C), 


and 


17 


£xa»\»'.ation  of  Eq.  (3^ )  reveals  that  the  stability  boundary  associated 

•%* 

with  the  singular  case  J=0  only  occurs  when  B=i,  which  is  the  already 


considered  z=l  case. 

th«  shbility  boundin' cs  of  Eqs.  >3$),  l~j&) ,  C^i) 

jy  t--1  ('  '  1  can  be  plotted  on  a  three-dimensional 
plot  with  axes  a,  b,  c.  A  sketch  of  these  bound¬ 
aries  is  shown  in  Fig.  t.  The  stability  region  is 


Fig.  2.  Stability  Boundaries  in  3-Dimensions 

found  by  applying  +V,t  '•>««' ^tVixwy.  rule?  or  by 
using  one  or  (.ore  test  points  (known  to  be  stable, 
to  lie  on  a, stability  boundary,  or  to  be  un¬ 
stable).  ‘  If  J>0  the  stable  region  (if  it 


exist'.)  lies  to  the  left  of  the  r.vr»W  conjugate 
root  stability  boundary  of  (gs.  (g-,j— ■ {/?> ,  as  0 
(or  c.nT )  increases;  the  left  side  of  l!  c  line  is 
double  cross-hatched  to  indicate  a  bound..  17  asso¬ 
ciated  with  double,  or  co.".;>)e/  conjugate ,  roots. 

If  J<0,  as  in  this  case,  the  stable  region  lies 
to  the  right  and  the  right  side  of  the  i'oc  is 
double  cross-hatched,  bir.glc  crr»ss-h*tcnir.g  is 
used  on  the  contours  of  tgs.  (Is)r(j4)  associated 
with  single  root  boundaries,  ire  side  of  tr.e 
contour  on  which  to  place  tne  cross-ha tchir.g  is 
dcterriir.cd  by  the  requi recent  that  cropp-natching 
be  continuous,  or  on  the  s are  side,  of  the  con¬ 
tours  as  the  intersections  corresponding  to 
2  *  +1  {0  5  0)  ar.o  =  -1  ('.  ■  oj  are  approached 
along  either  the  cr  root  or  real  rest  sta¬ 
bility  t 'under/,  t’sir.g  this  cri tenor.,  t re  starlc 
region  is  found  to  be  bounded  by  t'.e  c  *  Ci  plane, 
the  a  *  ?.  plate,  and  the  curve  sp'dfuJ  ty 
Er,s.  The  latter  m y  be  solved  for  a 

as  a  function  of  b,  resulting  In  the  sirple 
expression 


1 


V 

/ 

9 

I 

I 


I? 


The  limiting  values  for  a  an_  _ _ 

letting  0  approach  0  ar.d  njtfor  ire  case 
0-0,  represent  C* by  the  series  expansion 

B  ■  cos  0  ■  l-(02/2!  )+(C'fy4!)-{0®/6!  )+••• 


Then,  using  Eqs.  (*t).  (*3),  and  (Z6),  one  obtains 

lim  b  *  c 
0-0 

lim  a  •  *  lim  c  [(</02)-(2/3)  +  (O2/36)----,')~» 

0-0  0-0  :  fS3") 

■For  c  =  0,  Eg.  (27)  still  holds  (recognize  Kj  -  0 
is  really  the  case  one  is  interested  in).  Tr.en, 

■using  Eqs.(lT'c)  and  (LI,),  one  obtains 

Tima  =  Tim  (Ki03/<Jvu3)[2-(G2/2!)+(S4/<!)---] 

0*0  0*0  — - — - ; - 

(02/2!)-(0^/4 !)♦..• 

;  »  Tim  (KiO/2Jvu3)[2-(252/6)+(e4/72)----]'0 

0-0  j 

For  the  limiting  case  of  0  -  ita 

Tim  a  =  Tim  ft(l*8)c/ll-B)]+(j-B}}  =  2  J- 

0-fl  0  —  it  *  ** 


and 

Tim  b  *  Tim  (c+l-B)  »  c+2 

0  -  it  O-it 


,  M 


The  results  of  the  foregoing  are  sketched  graph¬ 
ically  in  two  dimensions  on  Fig..Jt  for  the  cases 
where  c  *  0  and  wnen  c  i  0. 


<ff  The  design  technique  e-ploys  the  pole  place¬ 
ment  approach.  It  is  based  on  the  premise  that 
th<~  dvr..?, tie  t-ehavicr  of  tre  system  ivy c sese ty 
tSi‘a';,-2rto  t..e  i'-  cation  or  t-s  roots' of  its  asse- 
cfUco  characteristic  conation.  Tito  r-otnod  shews 
both  analytically  and  graphically  >...o  direct 


correlation  between  these  roots  and  the  control 
gains  and  sample  period  of  the  controller.  The 
design  technique  then  involves  the  specification 
of  the  characteristic  equation  root  locations  ard 
the  subsequent  determination  of  the  control  system 
gains  and  sample  period  needed  to  attain  these 
locations.  The  control  system  designer  tr.en  rust 
determine  the  system  response  resulting  from  using 
these  numerical  values  and  assess  its  adequacy. 

If  it  is  not  adequate,  he  usually  relies  on  nis 
experience  to  re- 1  oca  to  the  roots  to  iniprove  tne 
response  in  the  manner  desired  for  his  particular 
system  (i.e.,  faster  settling  tire,  lower  peak 
overshoot,  etc.).'  It  is  assured  that  one  wishes 
the  pair  of  complex  conjugate  poles  of  the  C.E.  to 
dominate  tne  dynamic  response  of  the  system.  This 
response  will  ts.n  be  modi  tied  o  y  locating  the 
third  (real)  root. 


First  consider  the  pair  of  complex  conjugate  roots.  Cne  nay  use  Eos. 07 
(1ST),  substituting  a.b  for  kg  .  k,  ,  respective  lypUU->[i.lro  using  So.s. 

X 


r 


i 


/ 


i 


vo 


b  -  [(l-«)/1+4)]a«C2/0*5)]c-[(1-4)V0+4)] 

\1i) 

One  sees  that  £q.  C*i)  yields  a  straight  line  con¬ 
tour  of  b  vs.  a  for  a  given  value  of  c  and  a  given 
value  of  4.  Thus,  for  ?  given  value  of  c,  where 
4-<t££tour  of  Eg.  (8?)  crosses  the  c-ccntour  of 
EqjvIjMigbne  ootains  a  value  of  3  for  the  givc-n 
value  oT^cana  those  specified  values  of  '  and  o. 
From  Eg.  (|),  c  r.:ay  be  determined  as  a  function  of 
the  characteristic  equation  root  locations: 


C  *  (a?  +  8^  -  2a  ♦  1 )  ( 1  -  6 )/4 , 


'I'D 


where  a  and  3  are  the  real  and  imaginary  parts, 
respectively,  of  the  pair  of  complex  conjugate 
roots  of  characteristic  equation  (i).  For  a 
specific-^  location  (a.:/  ct  the  pair  of  complex 
conjugate  rests,  a  value  for  c  Is  established  by  ' 
setting  a  value  for  S  (location  of  the  real  root). 

An  additional  constraint  Is  imposed  by  Shannon's 
Sampling  Theorem: 

hi  <  w</2,  m 


where  u  is  defined  in  Eq.  (S )  and  uis/2ii  is  the 
Sampling  frequency, 


•  ias  *  2;i/T. 

From  Eqs.  (S"o)  -C?i)j 
may  be  restated  as 

0  <  it  „ 


[SI) 

this  constraint 
IS*) 


4 

The  stable  region  pprtrayc-d  in  Figs.  2-3  may 
be  shown  for  several  selected  values  of  c  on 
Fig.  ?J  using  Eqs.  (8  )  and  ft8).  As  c  increases 

in  value,  the  stable  region  shrinks  in  size  ar.d 
the  c  contours  lie  further  to  the  right.  «.r.en 
C  =  0.6,  the  r,  =  1//2  contour  lies  entirely  to  tha 

right  (ar.d  outside)  of  the  stable  rrcicr..  Finally.  _j  ,  a  u.- 

when  c_>.2,  the  stable  region  disappears  er.tirsly 

The  application  of  the  foregoing  technique  will 
now  be  surnarized. 

1.  Cased  on  desired  system  response 
characteristics,  tentatively  select  locations  for 
the  three  roots  of  '.re  characteristic  c-cuatisn. 

Although  the  zeros  of  the  system  characteristic 
equation  also  affect  the  dynamics,  placement  of 
the  poles  will  dominantly  affect  the  dynamics. 

a.  A  root  location  crite-ia  might  be 
to  select  numerical  values  for  c  ana  un. 

b.  Then  qualitatively  determine  hew 
much  integral  gain  is  desired  for  th*»  types  of 
disturbance  inputs  expected.  This  gives  an  indi¬ 
cation  of  the  value  of  c  to  use. 

c.  Determine  computer  cesign  con¬ 
straints.  Crc  migat  be  to  select  as  large  a  value 
of  T  as  possible  to  r  .r.-.-iizc  on-board  roiory 
capacity  requirements. 


r 


,/ 


2.  Using  the  value  of  c  tentatively 
established  in  paragraph  lb  above,  plot  the  sta¬ 
bility  Lounda-'es  in  r '•*?  b  vs.  a  parc-cter  plane, 
using  fc.s.  (3  4)--(3S- X  Thor,  plot  t.-.c  desirca 
t-contcur  (as  a  funcv-on  of  argor-e''*  ml)  on  the 
b  vs.  a  parameter  plane,  using  tqs. ( 19), 


3.  Using  various  values  of  4.  olot  6- 
contours  on  the  same  plot,  using  £q.  (/igj. 


4.  Find  an  intersection  of  a  4-contnur 
which  gives  a  desired  value  for  T  (ar.d  ner.ee  T, 
since  un  has  airejey  teen  specified)  and  4.  In 


/ 

j  general ,  if  4  •»  0,  its  effect  on  the  system 
/  d>nci»ics  is  small  compared  Co  the  effect  uf  the 


pair  of  complex  conjugate  poles  (placed  b.v  u!nT,  t). 


5.  Chech  the  response  of  the  system  using 
the  values  of  a,  b,  c,  and  i..,,T  (and  iier.ee  dp,  Kj , 
K0,  and  T)  associated  with  the  selected  inter¬ 
section.  If  it  is  unsatisfactory,  reiterate  the 
above  procedure,  selecting  another  value  of  c. 


If  the  design  procedure  were  merely  to  specify 
the  three  root  locations  (such  as  in  terms  of  a, 

6,  4  or  u'„T,  4),  iteration  would  not  he  required. 

Unfortunately,  in  design  practice  one  must  exercise 
engineering  judgement  in  selecting  the  tnree  root 
locations,  giving  rise  to  iterative  design  pro¬ 
cedures  such  as  the  one  outlined  above.  However, 
another- pos  si  bi-Li  J.y  does  exist.  If  the  real  root 
is  placed  near  the  origin  of  the  z-plar.e  (6  -  0), 
and  if  the  pair  of  complex  conjugate  roots  are 
located  near  the  unit  circle  in  the  z-plane  (i.e., 
specifying  and  uTl7 ),  the  latter  tv.'o  roots  will 
dominate  the  system  response.  "her.  the 
already  developed  tools  for  specifying  the  system 
response  (such  as  determining  explicitly  the 

.  maximum  overshoot  and  the  settling  tin?)  of  second-  .  .  _ 

'  order  systems  may  be  applied  handily.  i'fie  foregoing  procedure 


may 


be  amplified  through  the  application  of  realistic  numerical  values 


for  the  system  parameters.  ;,SSu-e  -mrr 

^ — s.  A«TeloHf' r yi."'  '»3i|»j»ic  a  digital  on-boar.'.' 

lor,  t1-'-  .v-Cur:  h;.r  assure  that  it  is  desired  to 


ppocosraft  slirfsr  le  tr.t-O-es'e* 

ccntrol- 


-  Q/-t-u  ru.ar  assure  tr.at  it  is  ....  __  .. 
Favn  a  controller  natural  frequency  (un)  of  six 
rad,s,  and  that  toe  rpsujr.t  of  inertia  (J7)  about 
the  single  axis  -  in  ij  ,T."'  ..~- 

^elweops-^'C  -i.ncii  is  -iQx 


ao'iei  is  IkJ  x  I  S'  'r.q-r.c.  Assume 
th3t  integral  control  is  cesirtd  to  drive  to  zero 
the  effect  of  constant  input  disturbances.  This 
means  a  r.on-zcro  value  of  c  is  desired.  !f  a  value 
of  r,  of  1//2  is  selected,  and  if  the  value  of  4  is 
desired  to  be  lopt  as  cicse  to  tero  as  possible, 
it  is  implied  from  Fig.  4c  tnat,  for  c  >  0,  r.o 
intersections  of  the  c-contuur  occur  when  c  >  0.2. 
Hence,  a  value  of  c  =0.1  is  chosen.  Tne  smallest 
value  of  4  whose  corresponding  c-line  intersects 


the  c-conteur,  with  a  reasoraole  factir  of 


t  Sc'Cvji 

is  6  *  0.4.  The  value  of  .pT  *  1.2  is  on  the  ;- 


contour  at  the  intersect". an,  yielci: 
0.20s  for  T.  Correspondin'!  values  0? 
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An  analytical  method  for  portraying  stability  regions  in  a- selected 
parameter  space  has  been  shown  for  a  digital  system.  The  method  requires 
that  the  system  characteristic  equation  be  availabla  -^d_  -acres  os 
in  the  complex  z -domain*  It  also  is  possible^.  ocxasn'ua sires  sqnamic 
characteristics  using  this  modified  parameter  space  technique.  The 
advantage  of  the  technique  over  existing  classical  sampled-data  methods 
is  that  the  stability  and  dynamic  response  characteristics  are  expres-ed 
in  terns  of  several  (rather  than  merely  one)  selected  parameters. 

Also,  the  sampling  psi'iod,  T,  need  not  be  expressed  numerically  before 
the  design  technique  begins,  giving  the  system  designer  one  more 
degree  of  freedom. 
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iFI-STDIX.  DaRIVAT.ICN  C?  3ICUR.il VS  RRI^TICKa 


It  is  possible  to  obtain  algebraic  recursive  relations  for  -the 
r ,al  and  ins. binary  parts  of  the  conplex  variable,  z,  when  it  is  raised 
to  the  power  (j  is  a  positive  integer)#  Fron  aq.  (3)»  z  is  asxine 
as  the  vector  (or  conplex  variable), 

~~ - -  z  =  eTs  =  r  eis=  r  cos  8  -g  r  sin  9  .  (Al) 


*1  Vy 

The  real  and  inarinary  parts  of  fined  as  Rj  and  I,  , 

respectively,  i.e. 


zJ  =  Rj  +  i  Ij  »  W^])  Kr 

Tf  i  =  i  .  an.  (.i?)  beccnes  identical  with  ~Qj  («1  A  loadina  to 
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APPENDIX  D.  PARAMETER  SPACE  PROGRAM:  STAB  HITT  REGIONS 


A  program  has  been  written  for  the  Hewlett  Packard 
9100B  Calculator  to  map  stability  boundaries  for  the  missile  system 
described  in  Section  7,  Part  C.  Since  the  real  root  boundary  is 
known,  only  the  complex  conjugate  root  boundary  need  be  plotted. 
This  ia  the  map  of  the  unit  circle  ffcom  the  complex  z -domain  onto 
the  selected  parameter  plane.  In  the  design  example,  the  control 
gains  kp  and  ky  are  selected  for  various  selected  numerical  values 
of  k»  (in  this  example,  kj  ®  0,  0.1,  and  0.5>),  Inputs  required  to 
start  the  program  are  the  numerical  value  of  kj_  (which  is  entered 
in  the  d  storage  register)  and  the  starting  value  of  «T  (which  is 
nominally  zero  and  entered  in  the  a  storage  register).  The  program 
automatically  steps  the  value  of  uT  upward  positively  in  increments 
of  one  (although  this  value  can  be  altered  by  inserting  a  different 
number  into  program  steps  Oc  through  10 ).  The  program  prints  out 
values  of  k_  (in  the  z  display)  and  kd  (in  the  y  display)  for  each 
printed  (inP the  x  display)  increment  of  w T. 
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APPENDIX  E.  SXSTH!  RESPONSE  PROGRAM 


Hie  purpose  of  this  program  is  to  compute  values  of 
the  missile  attitude,  0^,  at  each  sampling  instant,  t*kT.  3h  this 
case,  the  response  to  a  unit  step  input  has  been  computed,  using 
Equations  (17-1:7)  through  (lV-5>l).  The  program  has  been  written 
for  the  Hewlett  Packard  9100B  Calculator.  Hie  required  inputs  to 
the  program  are  the  selected,  numerical  values  of  the  control  gains, 
kp,  kj,  and  k^,  and  the  initial  conditions  on  0  at  t-0,  t*T,  and 
t*2T.  These  latter  three  are  calculated  from  Equations  (I7-U7) 
through  (17-1:9),  respectively.  The  program  can  be  easily  modified 
to  accept  initial  conditions  other  than  9o=0  and  for  inputs  other 
than  a  unit  step  input.  Hie  program  prints  out  values  of  9(kT) 
for  t  >  37  in  the  y*  display.  Hie  numerical  value  for  ki  must  be 
inserted  in  program  steps  97  through  9b.  Figure  h  indicates  a 
plot  of  the  outputs  of  this  program  for  selected  values  of  kp, 
kj,  and  kj  (see  Figure  3  for  rationale  in  selecting  these  values). 
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Abstrac  t 

A  comparison  of  guidance  laws  applicable  to 
short  range  tactical  missiles  is  made.  These  laws 
are  segmented  into  several  classes  and  the  princi¬ 
ples  underlying  each  class  are  discussed.  Specific 
attention  is  given  to  the  structure  of  the  guidance 
technique  and  the  requirements  for  its  implementa¬ 
tion.  Evaluation  and  comparison  of  the  performance 
of  each  guidance  law  versus  the  cost  of  implementing 
it  are  considered.  An  extensive  bibliography  of 
relevant  literature  is  included. 

I.  Introduction 

The  US  Army  Missile  Research  and  Development 
Cotmand  (MIRADCOM)  recently  began  a  task  to  develop 
an  advanced  guidance  and  control  system  for  further 
Army  Modular  Missiles.  The  intent  is  to  "leapfrog" 
systems  currently  under  development.  The  purpose 
of  this  paper  is  to  describe  the  work  that  has  been 
done  within  this  new  task  and  to  provide  an  indi¬ 
cation  of  future  efforts  that  are  now  planned. 

The  reason  for  embarking  on  this  task  is  now 
summarized.  Present  weapon  systems  performance 
may  be  seriously  degraded  in  engagements  against 
targets  with  predicted  characteristics  of  the  1990s 
and  beyond  and  in  the  battlefield  environments  of 
that  time  frame  (see,  e.g.,  p.  11,  Aviation  '.leek 
and  Space  Technology.  March  20,  1978).  It  has 
been  established  that  the  guidance  laws  currently 
in  wide  use  may  not  be  adequate  to  combat  those 
threats.  Thus,  it  is  projected  that  fundamental 
advances  in  guidance  and  control  systems  theory 
are  required  to  enhance  the  effectiveness  of  future 
weapon  systems.  Additionally,  missile  airframe 
and  propulsion  systems  may  require  advances  commen¬ 
surate  with  the  predicted  target  scenarios.  In 
particular,  air  defense  weapons  currently  in 
Research  and  Development  (R&D)  may  be  seriously 
hampered  in  the  combat  scenarios  env.sioned.  From 
an  overall  systems  viewpoint,  this  task  shall 
address  the  issue  of  creating  new  theory  in  the 
guidance  and  control  area  to  meet  the  high  perform¬ 
ance  threat  of  the  future  as  a  leading  technology 
item.  Closely  associated  with  It  and  in  parallel 
with  the  guidance  and  control  effort,  weapon  system 
work  shall  be  undertaken  to  modify  airframe  and 
propulsion  modules  to  be  capable  of  engaging  the 
1990s  threat.  General  support  weapons  shall  be 
viewed  initially  as  a  subset  of  the  air  defense 


system.  Heretofore  these  two  classes  of  weapons 
each  were  developed  independently.  This  research 
shall  attempt  to  view  them  as  potentially  similar 
systems  that  utilize  different  modules  such  as 
propulsion,  guidance,  warhead,  etc. 

The  first  step  in  implementing  this  task  was 
to  conduct  a  literature  survey  to  establish  a  tech¬ 
nology  base  starting  point.  Constraints  on  length 
of  the  paper  result  in  summarizing  (alphabetically 
by  author)  the  list  to  include  only  those  refer¬ 
ences  considered  by  the  authors  to  be  most  relevant. 

Following  this  survey,  guidance  laws  were 
placed  in  five  categories  and  defined  mathematically 
The  implementation  and  predicted  performance  of  each 
category  was  then  investigated  and  compared  in  light 
of  current  and  predicted  hardware  and  software  capa¬ 
bilities.  This  paper  describes  these  results. 

II.  Guidance  Laws 

The  development  of  guidance  laws  for  short 
range  tactical  missiles  has  become  a  well-researched 
topic  over  the  past  25  years.  A  summary  of  a 
detailed  literature  survey,  how  each  guidance  law 
can  be  implemented,  and  guidance  law  predicted  per¬ 
formance  are  described  within  the  five  guidance  law 
categories  stated  in  the  Introduction. 

Line- of- Sight  Guidance  I’Commanc-to-Line-of-Sight 
and  Beam  Rider) 

Clemow  (1960),  in  his  book  Missile  Guidance, 
provides  a  detailed  discussion  of  beam  riding, 
while  Mahapatra  (1976)  discusses  morphological 
design  based  on  beam  riding.  Clemow  (1960)  dis¬ 
cusses  command- to- line- of- sight  (CL0S) ,  while  Rain 
and  Yost  (1976)  use  this  method  in  their  ship 
defense  scenario,  using  Kalman  filters  to  reduce 
beam  jitter. 

At  the  MIRADCOM,  CLOS  and  Beam  Rider  (3R)  con¬ 
cepts  are  each  considered  a  subset  of  the  line-of- 
sight  guidance  laws.  They  differ  primarily  in 
their  mechanization.  The  CLOS  uses  a  wire  for  the 
transponder  link,  e.g.,  TOW  or  DRAGON.  3R  may  use 
an  electro-optical  link,  e.g.,  SHILLELAGH,  and  fly 
in  a  directed  beam  aimed  at  the  target.  Generically 
they  are  similar  and  will  be  discussed  as  one. 
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Th<*  )  ine-ot-sighC  (LOS)  guiii.'ivt*  scheme  (CLOS 
«“;-u  »*.R)  Is  one  in  which  the  missile  is  v.uiJed  on  an 
LOJ  cc-Uf&c  so  as  to  remain  on  a  line  adjoining  the 
f  r-.\  t  .*rul  the  point  of  control.  To  lly  along  the 
T.C*^,  ♦hi;  missile  requires  a  velocity  component  (V^) 

>*•  •  *  icular  to  the  LOS  that  is  equal  to  the  LOS 
Vti !  *v  i  ty  described  by  the  relation 

VM1  =  RSM  7  ST  (1) 

where  represents  the  range  from  the  missile  to 
the  tracking  station  and  represents  the  LOS. 

lr.  foneral,  the  missile  flies  a  pursuit  guidance 
cc'jl  se  at  the  initiation  of  the  entry  into  the  beam 
at  launch  and  flies  an  approximately  constant  bear¬ 
ing  course  near  impact.  This  is  observed  from  the 
velocity  equation 
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V  *  — ^  V  (2) 

Ml  RgT  VT  W 

where  the  range  of  the  tracking  station  to  the  tar¬ 
get  is  given  by  R  and  V  represents  the  target 
bT  T 

velocity  relative  to  the  surface  of  the  earth. 

Figures  1  and  2  are  simplified  control  block 
diagrams  highlighting  features  of  the  scheme  for 
each  of  the  two  cases.  It  should  be  noted  that  the 
projector  mount  dynamics  in  Figure  1  and  the  tracker 
mount  dynamics  in  Figure  2  may  be  considerably  more 


Go. ''.stein  (19/2)  gave  a  comparison  of  the 
qualiLic;  an &  sensitivities  of  LOS,  pursuit,  and 
proportional  navigation  guidance  for  air-to-ground 
and  air-to-air  missiles.  Tn  another  paper  in  the 
Sana  report,  Goodstein  discusses  the  guidance  and 
control  system  tradeoffs  in  missile  design. 

Two  pursuit  guidance  laws  are  discussed  herein: 
attitude  pursuit  guidance  and  velocity  pursuit 
guidance.  Attitude  pursuit  guidance  tries  to  keep 
the  centerline  of  the  missile  pointed  at  the  target. 
In  a  missile  which  flies  an  angle  of  attack  when 
maneuvering,  the  velocity  vector  will  always  lag 
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♦Implementation  scheme  provided  by  R.  H.  Fanner,  Technology  Laboratory,  US  Army  Missile  Research 
and  Development  Conmand. 
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the  vehicle  pointing  direction.  Miss  distance  is  a 
strong  function  of  the  maneuver  capability  of  the 
missile  and  can  be  reduced  by  a  fast  responding 
high-g  vehicle. 

Velocity  pursuit  guidance  attempts  to  keep  the 
velocity  vector  of  Che  missile  pointed  at  the  tar¬ 
get.  It  is  mechanized  in  some  less  sophisticated 
missiles  by  mounting  a  target  sensor  on  an  air  vane 
which  indicates  relative  wind  direction.  The  dif¬ 
ference  between  this  velocity  vector  and  a  true 
velocity  vector  is  the  primary  error  in  the  scheme. 
The  atticude  pursuit  guidance  mechanization 
decouples  the  angle  of  attack  from  the  target 
seeker  and  improves  miss  distance  performance  by 
an  amount  proportional  to  the  vehicle  angle  of 
attack.  Figure  3  depicts  two-dimensional  geometry 
useful  for  describing  the  pursuit  guidance  laws. 


Figure  3.  Pursuit  Guidance  Geometry 


Neglecting  the  case  of  a  maneuvering  target, 
for  simplicity,  one  notes  that 


VT  * 


and  the  LOS  rate  is 


(3) 


(4) 


For  an  ideal  pursuit,  0  -  X..  Since  0  *  the 
missile  will  always  have  to  turn  during  the  attack 
except  for  the  case  of  a  perfect  head  or  tail  chase. 

Figures  4  and  5  present  examples  of  simplified 
control  system  diagrams  for  the  attitude  and  velo¬ 
city  pursuit  guidance  lav/s,  respectively.  In  the 
former,  a  wide  angle  target  sensor  is  required 
since  it  is  typically  mechanized  to  be  body  fixed. 

In  the  latter,  a  narrower  f icld-of-view  (FOV)  sen¬ 
sor  may  be  utilized  as  a  result  of  the  decoupling 
of  the  body  from  the  sensor  mount  as  previously 
described.  In  each  case,  is  the  forward  guidance 

gain,  and  it  will  differ  for  each  as  will  the  feed¬ 
back  damping  gain  KR.  The  guidance  filter  in 

Figure  5  Indicates  that  higher  guidance  gain  in  the 
velocity  pursuit  law  requires  some  smoothing  to 
Inhibit  noise  of  the  target  sensor  optics  and  its 
associated  electronics. 
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Tigure  5.  Velocity  Pursuit  Guidance 

The  performance  to  be  expected  from  pursuit 
guidance  is  indicated  in  Figures  6,  7,  and  8.  These 
were  obtained  for  a  tactical  weapon  of  the  class 
known  as  close  support  antitank  weapons.  They  are 
indicative  of  the  quality  of  performance  one  may 
expect  for  these  guidance  laws.  Although  these  are 
simulation  results,  recent  experience  with  flight 
hardware  has  validated  the  simulated  performance  to 
a  high  degree  of  confidence. 

The  proportional  navigation  guidance  (PNG)  law 
performance  is  also  indicated;  therefore,  it  shall 
be  described  next. 


READING  ERROR  <d«> 


Figure  6.  Performance  of  Classical  Guidance  Laws 
for  a  Given  Heading  Error  at  First  Target 
Acquis itlon 
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Figure  7.  Performance  of  Classical  Guidance  Laws 
for  a  Given  Target  Velocity  and  as  a  Function 
of  Control  Authority  (Units  Expressed  in  g's) 


TARGET  ACCELERATION  (g> 

Figure  8.  Performance  of  Classical  Guidance  Laws 
for  a  Given  Target  Acceleration  Commencing 
5  sec  3cfore  Impact 

Proportional  Navigation  Guidance 

Adler  (1956)  provided  one  of  the  seminal  papers 
in  the  area,  considering  PNG  in  three  dimensions. 
Adler  provides  a  readable  development  of  the  theory 
using  vector  calculus,  whereas  most  subsequent 
authors  have  considered  only  the  planar  case  for 
simplification.  Murtaugh  and  Criel (1966)  provided 
another  fundamental  paper  developing  three-dimen¬ 
sional  PNG  for  a  satellite  rendezvous  problem. 


Clcmow  (1960)  gives  another  basic  derivation 
of  TNG ,  while  Pitman  (1972)  compiled  sensitivity 
functions  and  projected  errors  for  PNG  with  gains 
of  2,  3,  and  4.  Calculations  of  terminal  homing 
parameters  were  made  by  Rawlings  (1970).  In  a  1971 
paper,  Rawlings  considered  the  effects  of  saturating 
aerodynamic  surfaces  on  trajectories  flown  with  PNG. 
Many  authors  have  considered  augmenting  PNG  to 
account  for  target  accelerations.  Arbcnz  (1970) 
considered  making  the  closing  velocity  heading  race 
proportional  to  LOS  rate  and  developed  a  closed 
form  expression  for  a  moo i tied  PNG  law.  Siouris 
(19741  added  an  estimate  of  target  acceleration  to 
the  missile  acceleration  command  to  yield  an  aug¬ 
mented  PNG  law.  Guelman  used  geometric  arguments 
to  give  the  structure  of  the  missile  trajectory; 
he  showed  in  1971  that  PNG  will  almost  always  result 
in  an  intercept  for  a  constant  velocity  target  in 
1972.  For  constant  target  accelerations,  qualita¬ 
tive  trajectories  were  determined  and  target  acqui¬ 
sition  boundaries  assessed.  In  1976  Guelman  con¬ 
sidered  the  structure  of  trajectories  under  true 
proportional  navigation  guidance  (TPNG),  where  the 
commanded  acceleration  is  normal  to  the  LOS  rather 
than  the  missile  velocity.  In  this  work  he  showed 
that  TPNG  results  in  intercept  only  i.f  the  initial 
conditions  lie  in  a  well  defined  subset  of  the 
parameter  space.  Shinar  (1976)  considered  PNG  for 
a  rolling  missile  where  he  considered  the  cross¬ 
coupling  between  roll  and  the  control  system. 

Slater  and  Wells  (1973)  studied  optimal  evasive 
tactics  against  a  PNG  missile,  incorporating  a  lag 
into  the  missile  dynamics,  and  generated  two  stra¬ 
tegies  based  upon  different  optimality  criteria. 

A  comprehensive  study  of  classical  (PNG)  and  opti¬ 
mal  control  techniques  for  terminal  homing  of  cru¬ 
ciform  and  bank-to-turn  steering  missiles  was  per¬ 
formed  by  Salbirnie,  Sheporaitis,  and  Merriam  (1975). 

PNG  is  a  guidance  law  in  which  the  angular  rate 
of  the  missile  flight  path  is  directly  porportional 
to  the  angular  LOS  rate  of  change.  This  is  shown 
simplistically  in  Figure  9  for  a  two-dimensional 
case . 


Figure  9.  Geometry  of  Flight  Path 
The  geometry  of  Figure  9  suggests 

7m  -  N\  (5) 

where  represents  flight  path  angle  relative  to 

a  fixed  reference,  \  represents  the  LOS  relative 
to  a  fixed  reference,  and  N  is  the  navigation  ratio. 
A  general  expression  for  missile  acceleration  may 
be  written 

1  -  vM>J  +  V  T  (&) 

where  r,  represents  tctal  missile  acceleration; 

V„,  the  missile  speed;  y  ,  the  missile  flight 
n  n 


path  angle;  VHJ,  a  unit  vector  lateral  to  the 

missile  flight  path;  and  T,  a  unit  vector  tangen¬ 
tial  to  missile  flight  path.  The  definitions  of 
R  and  T  yield  the  relationship,  K  •  T  =  0.  One 
may  implement  proportional  navigation  via  the 
relationship  (in  the  Laplace  domain) 


i^Cs)  »  ke(s) 


t  kxs\f  s) ] 
(1  +  ts) 


(7) 


where  represents  Lhe  lateral  acceleration  in 

units  of  g's  and  c  represents  the  LOS  error  which 
is  measured  by  a  seeker  having  a  Lime  lag  constant 
of  a  sec  and  an  error  of  t\.  The  symbol  k  repre¬ 
sents  a  guidance  gain  factor. 


By  noting  that  a  constant  bearing  course  is 
determined  if  missile  and  target  are  flying  con¬ 
stant  speed  and  neither  is  maneuvering,  it  may  be 
concluded  chat  the  LOS  at  each  instant  of  time 
would  be  parallel  to  the  LOS  at  a  previous  inscant. 
Laterally  perturbing  the  collision  course  noted  by 
target  and  missile  positions  X^  and  X^  respectively 

yields  Z T  and  Z^  which  may  be  integrated  to  yield 

the  missile  trajectory  using  proportional 
navigation: 


+ 


dt 


-gkt\  +  t 


where  E(.)  des ignites  the  expected  value  operation, 
Sf  and  k  are  weighting  matrices,  x  is  the  state, 

u  is  the  control,  and  w  is  a  white  noise  process. 

In  most  formulations  the  dynamics  are  assumed  con¬ 
stant  to  obtain  closed  form  solutions. 

Bryson,  Denham,  and  Dreyfus  (1963);  Denham 
and  Bryson  (1964);  and  Denham  (1964)  were  among 
the  first  to  consider  optimization  techniques 
applied  to  missile  guidance  problems.  In  a  series 
of  papers  they  formulated  and  solved  the  optimal 
control  problem  with  inequality  constraints,  and 
then  applied  the.  result  to  the  trajectory  shaping 
of  a  surface-to-surface  missile  for  range 
maximization. 

Stallard  (1968)  gave  a  good  tutorial  review  of 
classical  and  modern  methods  for  homing  interceptor 
missiles.  In  an  earlier  paper  (1966),  he  dealt 
with  an  autopilot  design.  Ir.  a  1972  paper,  he 
applied  discrete  optimal  control  to  a  missile  sys¬ 
tem  with  undesirable  stability  characteristics, 
and  then  modified  the  dynamics  to  yield  a  new 
problem  with  more  desirable  characteristics. 

An  excellent  review  of  deterministic  optimal 
control  and  its  applications,  with  an  extensive 
bibliography,  is  found  in  an  IEEE  paper  by  Athens 
(1966).  Another  excellent  review  with  extensive 
references  is  in  a  1965  AIAA  paper  i  Paiewonsky. 
Kokotovic  and  Rutman  (1965)  provided  survey  of 
sensitivity  methods  drawing  heavily  upon  Soviet 
literature. 


n^L*“k,<x>-»  •  <■> 

The  performance  and  implementation  of  the 
guidance  law  are  best  appreciated  in  their  relation¬ 
ship  to  two  other  guidance  laws  which  are  similar  in 
in  that  no  requirement  is  established  for  range  or 
range  rate  information.  Velocity  pursuit  and  atti¬ 
tude  pursuit  are  in  this  category  with  PNG.  The 
performance  is  indicated  in  Figures  6,  7,  and  8. 

Optimal  Linear  Guidance 

Since  the  mid-1960's  the  missile  guidance 
literature  has  become  increasingly  permeated  by 
techniques  based  upon  optimal  control.  The  great 
success  found  by  linear-quadratic  regulator  theory 
and  its  dual  analog,  Kalman  filtering,  plus  the 
attractive  and  easily  determined  form  of  a  feed¬ 
back  solution  has  led  to  almost  all  work  in  this 
area  being  based  upon  linear  model  dynamics,  with 
quadratic  costs  and  additive  Gaussian  noise  (LQG). 
Most  formulations  consider  terminal  miss  distance 
and  running  control  effort  only  in  the  cost  func¬ 
tional.  Unlike  the  standard  regulator  format,  a 
running  cost  on  the  state  is  generally  not  appro¬ 
priate  In  this  framework.  The  general  optimization 
problem  then  becomes 


Some  early  authors,  in  an  attempt  to  justify 
the  resulting  guidance  laws  achieved  via  linear 
optimal  control,  showed  that  their  results  were 
an  extension  of  PNG.  Axelband  and  Hardy  (1969, 

1970)  used  linear  optimal  control  to  develop  what 
they  called  qu3Si-optimum  PNG.  As  with  almost  all 
linear  optimal  schemes,  time-to-go  is  required  for 
implementation. 

Lee  (1969)  described  several  techniques  in 
nonlinear  control  which  extended  the  linear  regu¬ 
lator  theory  and  drew  heavily  upon  recent  disser¬ 
tations  at  the  University  of  Minnesota. 

Deyst  and  Price  (1973)  used  linear  dynamics 
with  a  planar  engagement  to  develop  optimal  control 
laws  where  the  target  acceleration  was  a  first- 
order  Markov  process,  and  found  that  the  saturation 
of  control  surfaces  was  an  important  factor  in 
modeling  a  maneuver  limited  missile. 

Nazaroff  (1976)  formulated  an  LQG  approach  in 
which  he  assumed  extremely  simplified  missile 
dynamics  but  included  target  acceleration  and  jerk 
terms.  Stockum  and  Wiener  (1976)  assumed  exponen¬ 
tially  correlated  target  accelerations  in  generat¬ 
ing  an  LQG  guidance  law.  The  feedback  they  obtained 
on  projected  miss  distance  and  rate  was  analogous 
to  time-varying  PNG.  They  compared  their  guidance 
scheme  against  an  augumented  PNG  law. 


min  E(J) 
u(.) 


*’  (tf)S£x(tf)  + 
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u'  (t)Ru  (t)dt 
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Asher  and  Matuszewski  (1974)  considered  the 
optimal  control  problem  where  target  acceleration 
was  accounted  for  is  an  external  disturbance, 
resulting  in  a  tracking  formulation  of  the  regula¬ 
tor  problem.  They  constrained  the  final  miss  dis¬ 
tance  to  be  zero  but,  to  achieve  this,  needed  to 
know  the  target  acceleration  history  precisely. 


subject  to  x ( t)  «  Ax  (t)  +  3u(t)  +  Gw(t) 


Balbirnie,  Sheporaitis,  and  Merries  (1975) 
gave  a  means  to  obtain  weighting  matrices  which 
result  in  "classical  type"  control  gains  in  an  I.QG 
formulation.  Slieporaitis ,  Balbirnie,  and  Liebner 
(1976)  considered  a  quadratic  cost  on  the  angle  of 
attack.  They  used  a  second  order  Nowton-Raphson 
scheme  to  solve  the  resulting  optimization  problem. 

Speyer  (1976)  applied  his  1 inear-exponential- 
Caussian  (LEG)  controller  to  a  terminal  guidance 
problem.  Rather  than  minimizing  the  expectation 
of  a  quadratic  form  (as  given  previously), 

Speyer's  LEG  formulation  minimizes  E j,iexy0 pJ^j 

where  a  is  a  scalar.  The  dynamical  model  is  again 
linear.  Speyer  used  a  Kalman  filter  to  obtain 
estimates  of  state  variables  for  feedback.  His 
controller  did  not  enjoy  a  separation  principal; 
control  gains  depended  upon  the  filter  state 
covariance.  Speyer's  LEG  cost  functional  had  the 
effect  of  very  heavily  weighing  large  excursions 
and  thus  reduced  the  tails  of  the  terminal  miss 
distribution. 


position  variable  of  the  missile  projected  on  the 
ground;  Yj  is  the  derivative  of  Y^,  the  missile 

to  the  target  velocity  projected  on  the  ground; 
is  the  lateral  acceleration  of  the  missile; 

0  is  the  body  attitude  angle  of  the  missile;  and 
0 (  is  the  angle  of  attack  of  the  missile  shown  in 
Figure  10. 
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Youngblood  (1977)  used  very  simplified  linear 
dynamics  to  compute  inner  launch  boundaries. 
Youngblood's  report  contains  a  description  of  the 
Fletcher-Powell  functional  optimization  method. 

Fiske  developed  a  number  of  guidance  and  esti¬ 
mation  schemes  based  upon  LQG  formulations  with 
varying  degrees  of  model  complexity.  The  closed 
form  solutions  for  the  controllers  are  given  in  the 
report,  allowing  one  to  examine  the  effect  of  model 
parameters  on  the  gains.  Fiske  also  presents  a 
stochastic  guidance  law,  wherein  the  target  accel¬ 
eration  is  assumed  to  be  a  first  order  system 
driven  by  white  noise.  In  addition,  a  nonlinear 
law,  based  upon  nulling  projected  miss  distance 
over  a  single  control  interval  was  given.  York 
and  Pastrick  (1977)  looked  at  the  problem  of  mini¬ 
mizing  the  terminal  miss  distance  and  the  deviation 
of  the  missile  from  a  desired  orientation  at  the 
final  time.  A  formulation  was  given  for  a  system 
that  had  finice  time  delay.  In  fact,  the  increase 
and  decrease  in  time  delay  had  interesting  ramifi¬ 
cations  on  the  solution.  The  angle  of  attack 
assumption  was  investigated  and,  although  not 
solved  analytically  in  closed  form,  the  system  was 
derived. 

For  completeness,  the  following  example  sum¬ 
marizes  a  typical  optimal  control  law  formulation. 
The  geometry  of  the  tactical  missile- target  posi¬ 
tion  is  given  in  Figure  10.  Assume  that  the  angle 
of  attack  is  small  and  thus  can  be  neglected  (this 
assumption  will  be  considered  later),  and  choose 
the  following  set  of  variables: 


Figure  10.  Geometry  of  Tactical  Missile 
Target  Positions 

This  optimal  control  problem  will  have  a  con¬ 
troller  of  the  form 

«  ’  Vd  +  cY*d  +  cee  +  \\  <u> 

where  C^,  C^.,  Cg,  and  are  time-varying  coeffi¬ 
cients  chosen  to  minimize  the  cost  functional 


J  -  fd(t£)  +  70(tf)  +  S  J  u2(t)dt  .  (12) 

o 

For  the  case  where  the  angle  of  attack  probably 
cannot  be  ignored  e.g.,  for  the  larger  tactical 
missile,  the  system  of  equations  should  include  the 
angle  of  attack  3.  In  addition,  because  it  is 
feasible  to  achieve  only  a  small  angle  of  attack  at 
impact,  a  reasonable  performance  index  to  be  mini¬ 
mized  would  seem  to  be 
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where  Y^  M  the  position  variable  from  the  missile 
to  the  target  projected  on  the  ground;  Y£  is  the 
position  variable  of  the  target;  Y  is  the 

IQ 


The  performance  obtainable  from  the  optimal 
guidance  law  formulated  and  simulated  for  Equation 
(13)  was  comparable  and  even  better  than  the  PNG 
lav  in  terms  of  miss  distance.  Additionally ,  it 
had  the  added  feature  of  meeting  a  constraint  or 
the  impact  angle  which  the  PNG  law  could  not 
achieve.  In  particular  the  impact  angle  was  shown 
via  simulation  on  an  all-digital  6-D0F  missile 
simulation  to  be  within  l  deg  of  the  desired  impact 
angle  and  within  1  ft  of  the  desired  miss  distance. 
Other  researchers  have  corroborated  these  results. 

The  performance  of  any  realistic  optimal  con¬ 
trol  law  in  a  missile  application  is  dependent  on 
the  estimation  of  final  time  or,  equivalently,  on 
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timc-to-go.  Typically,  an  estimate  of  the  range 
between  the  target  and  missile,  and  the  rate  o£ 
change  of  this  range  are  obtained  from  radar  or 
other  ranging  devices;  the  time-to-go  estimate  is 
then  calculated.  This  estimate  works  quite  well 
as  long  as  the  range  and  range-rate  information 
are  accurate.  In  many  instances,  however,  the  data 
are  contaminated  by  noise  either  covertly,  as  in 
the  case  of  radar  jamming  devices,  or  by  the  pro¬ 
cessing  electronics.  This  adversely  impacts  the 
estimate  of  time-to-go  and  the  optimal  control  law, 
and  missile  performance  suffers.  Pastrick  and  York 
(1977)  present  a  discussion  of  several  aspects  of 
the  problem  in  the  context  of  a  realistic  applica¬ 
tion  and  provide  analytic  computer  algorithms  for 
its  solution,  as  well  as  a  closed-form  result. 
Another  more  recent  attempt  to  estimate  time-to-go 
was  made  by  Fiske.  He  also  addresses  the  possibility 
of  obtaining  this  variable  by  an  Intensity  ranging 
technique. 

Another  difficulty  with  optimal  guidance  laws 
is  their  sensitivity  to  initial  conditions,  as 
shown  by  York  (1978).  The  primary  message  therein 
is  the  strong  need  for  accurate  modeling  of  the 
system  and  the  importance  of  the  selection  of 
numerical  quantities  for  the  elements  of  the  weight¬ 
ing  matrices  in  the  performance  index. 

Other  Guidance  Schemes 

Whiting  and  Jobe  (1972)  considered  using  a 
"virtual  target"  approach  to  guidance  for  a  short 
range  air-to-ground  missile. 

Poulter  and  Anderson  (1976)  considered  a  non¬ 
linear  differential  game  framework  to  derive  an 
optimal  steering  law  for  a  terminal  homing  missile. 
In  a  prior  paper,  Anderson  (1974)  gave  a  method  of 
updating  a  differential  game  solution  via  linearized 
two-point  boundary  value  problems. 

For  background  material,  Fronlng  and  Gieseking 
(1973)  gave  a  description  of  autopilot  steering 
mechanisms  for  bank-to-turn  missiles.  Gido,  Jaffe, 
and  Wilson  (1974)  provided  computer  programs  to 
produce  autopilot  designs  using  both  classical  (PNG) 
and  modern  (LQG)  formulations.  George  (1974)  dis¬ 
cussed  trends  in  IMG,  guidance  and  control  hardware. 


Mahmoud  (1977)  described  a  dual  level,  hier¬ 
archical  optimization  scheme  based  on  invariant 
imbedding  which  is  used  to  obtain  approximate  solu¬ 
tions  to  many  nonlinear  control  problems.  Connor 
and  Vlach  (1977)  presented  a  new  augmented  penalty 
function  approach  to  optimal  control  problems. 

Their  formulation  is  applicable  to  finite  dimen¬ 
sional  optimization  problems  with  terminal 
constraints . 

Lansing  and  3attin  (1965)  presented  much  back¬ 
ground  material  on  random  processes  applied  to 
automatic  control.  Radbill  and  McCue  (1970)  pro¬ 
vided  background  material  on  quasi-linearizaticn 
methods  in  solving  coupled  nonlinear  two-point 
boundary  value  problems. 

Chin's  book  on  missile  design  (1961)  has  a 
section  on  missile  transfer  functions  which  is  very 
useful  in  assessing  the  contribution  of  aerodynamic 
factors  and  control  surfaces  to  missile  motion. 
Goodstein  (1972)  provided  an  overview  of  missile 
control  system  development  and,  in  the  same  AGARil 
report,  Acus  (1972)  presented  a  description  and 
comparison  of  inertial  guidance  technology  for 
tactical  missiles. 

III.  Discussion 

A  survey  has  been  made  of  guidance  laws  applied 
to  short  range  tactical  missiles,  and  they  have  been 
organized  into  five  categories  for  convenience  of 
discussion. 

To  better  place  them  in  relative  perspective, 
Figure  11  presents  the  hardware  requirements  for 
mechanizing  them.  Though  the  concepts  are  reason¬ 
ably  self  explanatory,  it  should  be  noted  that  the 
BR  and  CLOS  concepts  are  configured  with  aft  sensor 
and  optics  while  the  others  are  forward.  Also,  the 
optimal  guidance  scheme  suggests  the  need  for  a 
microprocessor  to  handle  the  more  complex  guidance 
algorithms . 

Guidance  concept  configuration  comparisons 
depicted  in  the  Table  feature  the  relative  cost  and 
other  items  that  show  complexity  and  performance. 
With  the  exception  of  the  microprocessor  requirement 
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Table.  Guidance  concepts  configuration  comparison 
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Guidance  Concept 

Accuracy 

(ft) 

Gyroscope 

Girabal  Mechanism 
(Sensor) 

Electronics 
(On- Joard) 

Control  Effort 

Relative  Cost 
(On-Soa  rd) 

Attitude  Pursuit 

>  30 

None 

None 

Same  J 

High  Turning 
Acceleration 

1.0 

Velocity  Pursuit 

>  20 

None 

Air  Vane 

Same  ) 

Required 

1.4 

PNG 

<  5 

Single 

2  DOF 

Gyro 

Gyro 

Same 

Low  Turning 

Acceleration 

Required 

1.6 

Beam  Rider 

<  2 

Attitude 

None 

Same 

0.9 

CLOS 

<  2 

Attitude 

None 

Same 

0.8 

Optimal 

< 

2  DOF 

Gyro 

Gyro 

Micro¬ 

processor 

High  Turning 

1.6 

for  the  optimal  guidance  law 
are  relatively  comparable.  1 

all  onboard 
However,  it 

electronics 
is  important 

IV. 

Conclusion 

to  note  that  Che 

CLOS  and  3R 

concepts  have  vastly 

When  the  performance  comparisons 

and  antici- 

more  complex  ground  stations  inherent  in  their 
mechanisations.  For  this  reason  alone,  their  air¬ 
borne  complexity  is  relatively  low  and  :.ess  expen¬ 
sive  compared  to  the  others  implying  that  their 
ground  trackers  are  not  included  in  that  cost 
estimate. 

Barring  unforeseen  breakthroughs  in  the  fifth 
category  (Other  Guidance  Schemes) ,  comparison  of 
performance  characteristics  favors  PNG  for  defeating 
targets  that  are  either  stationary  or  moving  with  a 
constant  velocity.  PNG  is  simple,  requires  a 
state  measurement  that  is  easily  attainable,  and 
is  relatively  easy  to  implement.  However,  for 
highly  maneuverable  accelerating  (or  decelerating) 
targets,  optimal  guidance  laws  are  superior, 
performance-wise,  to  other  forms  of  guidance. 
Further,  if  miss  distance  is  the  only  criterion  of 
performance,  PNG  is  favorable.  If  other  performance 
criteria  are  important,  such  as  impact  angle  at  the 
target,  then  these  criteria  can  be  Incorporated 
into  the  performance  index  and  optimal  guidance 
becomes  superior. 


pated  implementation  difficulties  are  considered, 
the  following  conclusion  is  reached.  PNG  is  a 
strong  contender  but  cannot  cope  with  all  expected 
targets  of  the  future.  Therefore,  optimal  guidance 
is  the  chosen  candidate  guidance  law,  even  with  its 
added  complexity  and  current  implementation  diffi¬ 
culties.  However,  these  may  be  alleviated  with 
judicious  use  of  the  microprocessor  technology  now 
enjoying  a  high  state  of  development  activity.  Con¬ 
siderable  effort  will  have  to  be  devoted  to  over¬ 
coming  what  appears  to  be  its  strongest  short¬ 
coming:  measuring  or  accurately  estimating  the 
time- to- go. 
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II.  THEORY  AND  DESIGN 


By  Sherman  H.  Seltzer 

1.  oraafaattm* 

In  an  attempt  to  provide  the  reader  with  a  logical  approach  to  the  review 
of  the  theory  and  design  associated  with  the  SIG-D  Program,  the  material  comprising 
that  review  is  presented  in  the  following  sequence: 

a*  Program  objectives , 

b.  Design  approach, 

c.  Analysis  program, 

d.  Simulation  program, 

e.  Hardware  testing  program, 

f.  Proof  of  concept, 

g.  Schedule, 

h.  Summary  and  conclusions. 

2.  Program  objectives. 

As  stated  by  Jack  Clayton  during  the  February  8-9,  1979  review,  "The  SIG-D 
Program  is  an  exploratory  development  effort  to  investigate  the  technical  risk  areas 
associated  with  a  low  cost  surface  to  surface  missile  utilizing  a  strap-down  in¬ 
ertial  navigation  and  guidance  system.*  If  it  is  intended  to  use  SIG-D  on  missiles 
other  than  the  particular  modified  LANCE  upon  which  design  has  been  based,  certain 
system  alterations  probably  will  be  required.  As  described  in  the  presentation, 
SIG-D  is  tied  closely  to  the  characteristics  of  the  modified  Lance  system  on  which 
it  is  carried.  In  particular,  the  bending  characteristics  of  the  modified  LANCE 
have  strong  impact  upon  the  SIG-D  design.  Bending  filter  design,  computer  storage 
capacity,  the  isolation  system,  and  the  selected  sampling  period  of  the  computer 
probably  all  will  be  affected  if  SIG-D  is  used  to  guide  a  missile  other  than  the 
modified  LANCE. 

3.  Design  approach. 

The  design  approach  is  orderly,  logical,  and  being  implemented  by  appar¬ 
ently  superb  engineers.  The  innovative  analytical  techniques  that  have  been  brought 
to  bear  on  the  SIG-D  design  are  outstanding.  An  exception  to  the  conservative  ap¬ 
proach  demonstrated  throughout  concerns  computer  memory  capacity.  It  was  stated- 


that  a  margin  of  75  wards  of  memory  still  remains.  I  feel  this  is  too  small  a  margin 
to  be  entering  into  the  hybrid  and  hardware- in-the-loop  (HWIL)  stages  of  (potential) 
design  refinement.  If  additional  memory  should  become  required,  it  may  be  necessary 
to  physically  add  an  additional  memory  module.  The  possibility  of  decreasing  memory 
requirements  by  performing  a  sampled-data  analysis  with  the  aim  of  decreasing  sampling 
period  might  be  considered  as  an  alternative*  Bending  mode  constraints  then  would 
have  to  be  considered  also. 

At  first  I  planned  to  recommend  looking  deeper  into  the  potential  problem 
posed  by  heretofore  undetected  sources  of  vibration.  An  example  (the  only  one  that  I  can 
think  of,  but  there  may  be  others)  is  the  external  cable  conduit  (or  "tunnel").  If 
it  were  to  become  unbonded  from  the  skin  of  the  missile  during  the  guided  portion  of 
the  flight,  its  vibration  characteristics  might  have  an  undesirable  dynamic  effect  on 
the  G&C  system.  Upon  further  investigation  I  found  that  this  problem  has  been 
studied  thoroughly.  In  the  case  of  the  cable  conduit,  the  masses  concerned  are  re¬ 
latively  negligible. 

A  design  decision  made  at  the  outset  of  the  program  was  to  use  the  Ring 
Laser  Gyroscope ,  presumally  because  of  its  low  cost.  According  to  Mr.  Clayton,  it 
has  undergone  extensive  sled  and  missile  ( "UPSTAGE  Program”)  testing,  although  never 
in  the  strapdown  mode  of  operation.  This  should  be  adequate. 

4.  AatiJBis  sraassi- 

As  already  indicated,  the  analysis  used  has  been  thorough  and  excellent 
in  quality.  Although  I  personally  am  wary  of  using  continuous  time-domain  techniques 
to  approximate  sarapled-data  dynamic  phenomena,  I  believe  it  is  Justified  in  the 
case  at  hand  because  of  the  relatively  high  sampling  rate  (the  latter  being  set  by 
the  frequency  of  the  highest  dominant  bending  mode).  Even  the  use  of  a  continuous 
time-domain  describing  function  analysis  appears  to  have  been  adequate  to  predict 
the  existence  and  characteristics  of  limit  cycles  and  was  stated  to  have  been  veri¬ 
fied  by  simulation.  This  is  a  tribute  to  the  engineering  Judgement  of  Bd  Herbert 
and  Or.  Paul  Jacobs  in  particular. 

Apparently  the  design  team  has  searched  for  all  dynamically  significant 
nonlinearities  and  has  identified  only  the  pneumatic  actuator  nonlinearity.  Cert¬ 
ainly  a  complete  HWIL  simulation  will  show  if  the  search  results  were  correct  (and 
I  think  they  are). 


If  tine  were  available,  I  would  recommend  a  mare  conventional  sampled- 
data  analysis  of  a  simplified  (far  analytical  tractability)  model  of  the  system. 

At  this  late  date  I  think  a  careful  simulation  should  be  adequate. 

The  IMU  mounting  plate  and  vibration  isolation  system  are  critically 
important  to  the  performance  of  the  G&C  system.  This  fact  has  been  duly  recognized, 
and  supporting  analysis  has  been  conducted.  Extensive  analysis  (supported  by 
hardware  tests— see  section  6,  below)  has  been  performed  to  validate  the  system 
design,  using  the  NASA  Structural  Analysis  (NASTRAN)  program.  This  analysis  in¬ 
cluded  the  identification  of  a  number  of  the  relatively  lower  longitudinal  and 
torsional  normal  (free-free)  bending  modes.  The  analysis  was  performed  for  both 
loaded  (i.e.  "powered"  flight)  and  unloaded  ("post-cut-off")  motor  case  conditions. 

5.  fijinmlrtim  Br9prwi» 

Usually  simulation  and  analysis  are  held  to  a  minimum  by  program  man¬ 
agers.  In  this  instance,  there  appears  to  be  an  over-abundance  of  simulation  pro¬ 
grams,  although  this  situation  is  the  result  of  an  orderly  progression  from  all- 
digital  simulations  toward  an  HWTL  simulation.  There  was  sene  indication  that 
agreement  between  the  various  M3RADCCM  and  Vought  programs  does  not  exist.  The 
reasons  far  this  disagreement  were  explained  during  and  subsequent  to  the  formal 
presentation.  Of  course,  ultimate  agreement  must  be  sought.  The  hybrid  and  the 
butt,  simulations  are  extremely  valuable  design  tools.  It  is  my  opinion  that  not 
enough  time  is  allocated  for  their  use.  It  is  not  apparent  why  the  rush,  when 
much  valuable  inf  carnation  should  be  forthcoming  from  the  use  of  these  two  simulations. 

6.  BfflftflTft  .testing  gggaa* 

Apparently  the  pneumatic  actuator  characteristics  have  been  confirmed 
by  actual  test.  It  also  was  stated  that  the  isolator  characteristics  have  been 
confimed  by  test. 

In  support  of  the  vibration  analysis  (see  para. 4,  above),  vibration 
testing  was  performed  on  a  suspended  post  cut-off -configuration  (i.e.  empty  motor 
ease)  of  the  modified  LANCE  missile.  There  was  reasonably  close  agreement  between 
analytically-derived  results  and  these  test  results.  Apparently  monetary  and 
schedule  constraints  precluded  a  test  of  a  simulated  boost  configuration  (i.e. 
simulated  propellant  characteristics)  as  well  as  tests  to  confirm  predicted  to¬ 
rsional  characteristics.  In  view  of  the  previously  noted  agreement  between  an¬ 
alytical  and  test  longitudinal  vibrations,  this  is  probably  an  acceptable  risk. 


J 
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7.  Proof  of  concept. 

If  the  tine  allotted  to  HWIL  and  hybrid  simulations  is  lengthened ,  and 
if  additional  vibration  analysis  is  performed,  there  appears  no  doubt  that  the 
design  concept  can  be  proven  sufficiently  for  flight  testing  to  begin.  With  a 
few  minor  exceptions,  simulation  results  appear  to  agree:  these  exceptions 
should  be  explained.  It  was  not  shown  if  the  several  error  analysis  discussed 
were  in  agreement. 

8.  Schedule. 

As  already  indicated,  it  is  my  opinion  that  the  schedule  is  too  tight 
to  permit  adequate  hybrid  and  HWIL  simulation. 

9.  Summary  and  conclusions. 

The  review  was  extremely  well-conducted  and  presented  ample  evidence  of 
a  thorough,  well-ccncieved  and  implemented  design.  The  SIG-D  design  team  is  to  be 
commended  for  their  fine  work.  The  following  recommendations  summarize  the  few 
shortcomings  associated  with  the  SIG-D  theory  and  design. 

a.  Prepare  a  contingency  plan  to  cope  with  the  possibility  that 
insufficient  computer  storage  capacity  may  exist  (para.  3)* 

b.  Conduct  a  careful  analysis  of  the  dynamic  effects  of  the  digital 
phenomena  as  manifested  in  the  HWIL  and  hybrid  simulations*  Be  careful  to  discrim- 
mlnate  between  simulated  on-board  sampling  phenomena  and  those  sampling  phenomena 
that  arise  because  of  the  digital  computer  used  far  the  simulation  (para.  4). 

c.  Allot  mare  time  in  the  schedule  for  an  orderly,  indepth  use  of 
the  hybrid  and  HWIL  simulations  (para.  5). 
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1.  On  10  August  1978,  engineers  from  the  Guidance  and  Control  Analysis  Group  were 
Invited  by  personnel  of  the  HELLFIRE  Project  Office  to  discuss  the  status  of  the 
group's  KWIL  simulation  capability. 

2.  A  paper  entitled,  "Utilization  of  Hardware-in-the-Loop  Simulation  for  HELLFIRE 

Tasks,"  has  been  prepared.  It  describes  how  the  Guidance  and  Control  Analysis  Group's! 
simulation  facility  Is  used  to  meet  that  group's  HELLFIRE  mission,  including  the  , 
support  of  the  Low  Cost  Laser  Seeker  evaluation.  | 
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UTILIZATION  OF  HARDWARE- I N-THE-LOOP  SIMULATION  FOR  HELLFIRE  TASKS 

Guidance  and  Control  Analysis  Group 
Guidance  and  Control  Directorate 
Technology  Laboratory 


1 .  Purpose 

The  purpose  of  this  paper  is  to  describe  how  the  Guidance  and  Control 
(G&C)  Analysis  Group's  Hardware-in-the-Loop  (HUIL)  simulation  facility  is 
used  to  meet  that  Group's  HELLFIRE  mission.  In  particular,  the  manner 
In  which  the  G&C  Analysis  Group  meets  the  objective  outlined  In  the  "Low 
Cost  Laser  Seeker  Evaluation  Test  Plan"  will  be  addressed.^ 

It  is  the  intent  of  the  G&C  Analysis  Group  to  support  the  evaluation 
of  the  two  low  cost  laser  seekers  currently  under  consideration  for  use  in 
the  HELLFIRE  Modular  Missile  System  (HMMS).  This  evaluation  is  being  conducts- 
under  the  direction  of  the  Optical  Guidance  Technology  Group,  Advanced 
Sensors  Directorate,  DRDMI-TEO.  That  Group  is  responsible  for  determining 
how  well  each  candidate  seeker  meets  the  seeker  specifications.  The  G&C 
Analysis  Group  is  responsible  for  comparing  the  two  seekers  in  a  facility 
that  simulates  the  missile  system  and  its  dynamics  under  predicted  environ¬ 
mental  conditions.  Further,  the  G&C  Analysis  Group  must  maintain  the 
technical  capability  to  monitor  and  spot-check  the  contractor's  [Rockwell 
International,  Columbus  (R.I.C.)]  technical  development  efforts  and  to 
assess  his  technical  approach,  particularly  when  he  recommends  altering  it. 

It  is  the  responsibility  of  R.I.C.  to  determine  quantatively  how  well  each 
candidate  seeker  meets  the  seeker  specifications. 

2.  Background 

A  comprehensive  HWIL  simulation  was  developed  to  perform  guidance  and 

control  system  design  verification  and  optimization  of  both  autopilot  and 

2 

seeker  for  a  laser  semi -active  terminal  homing  weapon  system. 

The  facility  consists  of  the  following  major  components: 

a.  A  Xerox  Sigma  5  digital  computer  with  64K  words  (32  bits)  of 
core  memory.  The  major  tasks  performed  by  this  computer  are  to  simulate 


the  missile  aerodynamics,  and  rotational  and  translational  dynamics,  and 
kinematics,  and  to  provide  commands  to  the  gimballed  mirror  and  the  laser 
designator. 

b.  Two  Comcor  Ci  5000  analog  computers  interfaced  with  digital 
computer  provide  a  hybrid  computer  simulation.  The  analog  computers  simu¬ 
late  missile  rotational  dynamics  and  can  simulate  autopilot,  actuator,  and 
seeker  dynamics  when  some  or  all  of  these  items  are  not  Included  within  the 
overall  HWIL  simulation. 

c.  A  laser  target  simulator  consisting  of  a  pulsed  laser  beam 
which  Is  reflected  onto  a  translucent  screen  by  a  gimballed  mirror.  The 
changing  distance  between  the  missile  and  the  target  is  simulated  by  varying 
the  laser  spot  size  and  the  energy  density. 

d.  The  three-axis  electrohydraulic  gimballed  Carco  flight  table 
simulates  missile  angular  motion.  The  laser  seeking  tracker,  which  senses 
the  laser  spot  through  the  translucent  screen.  Is  mounted  on  this  table. 

This  component  provides  steering  signals  to  the  actual  or  simulated  auto¬ 
pilot  and  actuators. 

The  HWIL  simulation  was  first  put  together  in  1969  to  evaluate  the 
(then)  SAM-D  missile.  It  originally  consisted  of  an  SDS  930  digital  computer 

Interfacing  with  a  Comcor  Ci  5000.  The  Carco  flight  table  also  was 
Incorporated  into  the  simulation.  In  1970-71  the  facility  was  improved 
to  permit  the  evaluation  of  the  Terminal  Homing  Accuracy  Demonstration 
(THAD),  which  was  the  predecessor  of  HELLFIRE.  The  improvement  consisted 
of  Incorporating  a  laser  terminal  homing  capability.  By  early  1971,  the 
facility  consisted  of  the  Comcor  Ci  5000,  a  newly  incorporated  Sigma  5 
digital  computer,  the  Carco  flight  table,  and  the  laser  simulator.  In 
the  1971-72  time  frame,  the  facility  was  used  to  evaluate  two  competing 
(Texas  Instruments  and  Martin-Marietta,  Orlando)  proposed  COPPERHEAD 
systems.  In  May  1973  it  took  on  a  true  hardware-in-the-loop  character. 

Both  Texas  Instruments  and  Martin  provided  12  missiles  each  to  be  evaluated 
competitively  for  the  advanced  development  phase  of  the  COPPERHEAD  program. 
Also  when  these  two  companies  bid  for  the  engineering  development  phase  of 
the  COPPERHEAD  program,  the  HWIL  facility  was  used  to  help  evaluate  the 
two  potential  contractors'  technical  capabilities.  Finally,  it  was  used 
to  support  the  source  selection  evaluation  board  (SSEB)  in  evaluating  the 
two  (R.I.C.  and  Hughes)  competing  bidders. 


Presently,  In  addition  to  HELLFIRE  support,  the  COPPERHEAD  program 
also  is  supported  using  the  HWIL  facility.  Each  produced  COPPERHEAD 
guidance  and  control  unit  (including  the  seeker  and  the  autopilot)  on  the 
Carco  table  and  subjecting  It  to  a  standardized  test  at  each  of  three 
roll  angles,  from  which  it  Is  determined  whether  the  unit  is  functionally 
operable. 

3.  HELLFIRE  Program  Support 

The  HWIL  facility  is  being  used  to  support  (by  performing  evaluations 
and  selected  analyses)  the  engineering  development  phase  of  the  HELLFIRE 
program.  It  is  also  intended  to  use  the  facility  to  evaluate  and  compare 
the  two  (R.I.C.  and  Martin)  competing  seeker  systems.  Rockwell  will 
evaluate  the  two  competing  laser  seekers  by  measuring  their  characteristics 
and  simulating  their  in-flight  operation  on  their  own  HWIL  facility.  The 
G&C  Analysis  Group  will  evaluate  each  of  the  two  laser  seekers  by  using 
their  HWIL  facility  and  determining  which  system  performs  best  under 
simulated  flight  conditions. 

The  G&C  Analysis  Group's  HWIL  simulation  facility  has  the  capability 
of  permitting  that  group  to  meet  Its  assigned  HELLFIRE  mission  of 
evaluating  and  comparing  the  two  competing  laser  seekers.  As  such,  it 
can  provide  Information  which  will  enable  the  group  to  determine  which 
system  has  the  best  chance  of  meeting  the  total  system  requirements. 

The  flight  table  used  In  the  HWIL  simulation  facility  is  a  modified 
2-axis  (pitch  and  yaw)  Carco  model.  The  bandwidth  of  that  table  has  been 
found  to  be  adequate  for  the  class  of  seekers  being  tested.  For  example, 
a  30°  phase  shift  existed  at  2  Hz  in  the  original  uncompensated  flight 
table.  However,  for  several  years  the  flight  table  has  been  compensated 
to  provide  the  proper  (i.e.,  0°)  phase  shift  at  that  frequency.  Performance 
has  been  established  over  a  60  db  range,  which  meets  the  specification. 

4.  Proposed  Hardware-in-the-Loop  Modifications 

In  addition  to  the  existing  modifications,  laboratory  conditions 
surrounding  the  HWIL  are  constantly  evaluated  with  the  objective  of 
improving  facility  performance  in  a  cost  conscious  manner. 


Modest  improvements  have  been  incorporated  from  time  to  time, 
such  as  the  recent  addition  of  anti-reflective  material  (for  the  wave 
length  we  are  using)  for  use  on  the  floor.  A  list  of  other  modifications 
under  consideration  follows.  Several  of  these  improvements  could  be 
Implemented  in  time  for  the  impending  seeker  evaluation. 

(1)  The  room  containing  the  facility  could  be  renovated  to  remove 
contaminating  inputs  to  the  systems  being  evaluated.  Where  feasible, 
mobile  equipment  that  does  not  need  to  be  in  the  room  (such  as  work  benches 
and  consoles)  could  be  removed  and  placed  outside  the  room.  The  ceiling 
and  walls  could  be  covered  with  an  anechoic  material  to  absorb  unwanted 
inputs.  A  sliding  screen  could  be  installed  to  separate  the  laser  area 
from  the  seeker  -  Carco  table  area.  Personnel  not  absolutely  needed  in  the 
room  during  simulation  runs  should  remain  outside  the  room.  Those  who 
remain  in  the  room  should  remove,  cover,  or  mask  light-colored  reflective 
clothing.  This  entire  set  of  improvements  could  be  implemented  in  a  short 
time  period. 

(2)  A  "Warm-up"  period  could  be  used  to  heat  up  equipment,  build 

up  energy  levels,  and  in  general  permit  transients  to  die  out.  This  should 
be  done  before  using  laser  pulses  in  order  to  provide  better  repeatability 
of  results.  This  improvement  could  be  implemented  immediately. 

(3)  A  study  could  be  made  for  improving  the  dynamic  response  of  the 
Carco  table  when  the  roll  fixture  is  attached. 

(4)  A  similar  (to  the  above  item)  study  could  be  instigated  to 
determine  If  the  mirror  system  dynamics  should  be  improved. 

(5)  An  analytical  study  of  a  means  of  improving  the  dynamic  range 
of  the  Carco  table  is  being  conducted  recently  be  Dr.  S.  M.  Seltzer  (a 
consultant  to  the  G&C  Analysis  Group). 

5.  Conclusions 

The  hardware-in-the-loop  facility  will  continue  to  be  used  to  support  th~. 
HELLFIRE  program.  Selected  scenarios  and  trajectories  that  are  within  the 
capabilities,  geometry,  and  optics  of  existing  equipment  will  be  utilized. 

It  is  felt  that  by  so  doing,  the  stated  HELLFIRE  mission  can  be  accomplished 
by  the  Guidance  and  Control  Analysis  Group. 
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1.  As  a  result  of  target  predictions  and  meetings  among  the  G&C  Analysis  Group  and  I 
the  Aeroballistics  Directorate  personnel,  the  following  decisions  regarding  the 
Advanced  G&C  System  have  been  arrived  at:  j 

a.  As  a  design  goal,  the  shape  of  the  missile  that  is  to  carry  the  Advanced  G&C  ’ 

System  is  to  be  as  clean  aerodynamically  as  possible.  Three  possible  shapes  are  to  be; 
considered  by  the  Aeroballistics  Directorate  in  their  studies  supporting  the  develop-  j 
moit  of  the  Advanced  G&C  System.  They  are:  j 

1)  Lew  aspect  ratio  long-cord  delta  configuration  (Fig.  1) .  This  shape  is 

similar  in  appearance  to  a  dart.  Maneuvering  is  provided  by  hinged  aft-tabs.  j 

2)  Flared-skirt  extendable  wine  surfaces  configuration  (Fig.  2) .  This  con¬ 

figuration  derives  its  stability  from  the  flared  aft-body,  wing  surfaces  are  extended- 
when  control  (maneuver)  authority  is  desired  —  typically  at  the  terminal  phases  of 
the  trajectory.  j 

3)  Lifting  body  configuration  (Fig.  3) .  This  is  to  obtain  more  efficient 
coordinated  ("bank- to  turn")  control  maneuvers .  Body  lift  is  achieved  by  using  ellip¬ 
tical  or  triangular  half  cross-sections. 

b.  Design  Point. 

The  following  design  point  will  be  used  by  Messrs.  Dave  Washington  and  Jimmie 
Derrick  (Aeroballistics  Directorate)  to  generate  aerodynamic  coefficients  for  G&C 
Analysis  Group  simulations.  { 


1)  MACH  No:  4 

2)  Lateral  g's:  12 

3)  Altitude:  70,000  ft 

It  is  expected  that  these  values  will  be  exceeded  by  a  significant  margin,  e.g. 
perhaps  20  g's. 

« 

2.  When  the  Aeroballistics  studies  have  been  completed,  the  G&C  Analysis  Group's 
simulation  development  will  begin.  The  analytical  models  of  the  plant  developed  by 
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Ed  Herbert  and  Mike  Kalanqe  (to  be  reported  upon  in  the  first  annual  report  on  our 
progress  and  achievements,  due  1  October  1979)  will  be  used  in  conjunction  with  the 
aerodynamic  data  as  the  basis  for  the  simulation. 
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LOW-ASPECT-RATIO  LONG-CHORO  DELTA  MISSILE  CONFIGURATION 


Figure  1 


FLARED- SKIRT- STA0ILIZEO  MISSILE  CONFIGURATION  WITH 
ALL-MOVABLE  FORWARD  SURFACE 


Figure  2 


